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A method is described of producing any desired con- 
centration of the isotope of hydrogen, H?, by electrolysis. 
Starting with twenty liters of water from an old com- 
mercial electrolytic cell and electrolyzing in four stages 
until only one-half of a cubic centimeter of water remained, 
this water had the specific gravity 1.073. Having shown 
that no large accumulation of the heavy isotopes of 
oxygen occurs and assuming that the density varies 
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linearly with the fraction of H®, 2/3 of the hydrogen in this 
water is H®. It is shown that a further reduction of volume 
to one-quarter will give 99 percent H*. In our electrolyses 
the percentage loss of H! is to the percentage loss of H? in 
the ratio of five to one. By two methods we have found 1 
in 6500, as a provisional value for the concentration of H? in 
ordinary water. 


INCE the announcement! of our success in 
attaining a high concentration of hydrogen 
isotope has aroused some curiosity as to the 
methods that were employed, we shall describe 
here our procedure, although we should have 
liked first to have made a more scientific analysis 
of the various steps in the process, and although 
it seems probable that better methods of pre- 
paring the isotope will shortly be available. The 
fact is, the difference in properties between the 
two isotopes of hydrogen is so much greater than 
between any other pair of isotopes that in spite 
of the very small amount of H? present in ordi- 
nary hydrogen, several methods will lead to an 
almost complete separation of H? from H!, as we 
have convinced ourselves by a number of 
preliminary experiments. 

Being engaged some months ago in the 
separation of the isotopes of oxygen, we found it 
Possible to follow small changes in the ratio 
of the isotopes in water by careful measurements 
of the density. By finding the temperature at 


‘Lewis, J. Am. Chem. Soc. 55, 1297 (1933). 
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which a 10 cc float will neither sink nor rise in 
ordinary pure water, and then finding the 
temperature at which this buoyancy balance is 
again established in a new sample of water, the 
specific gravity of that sample may be obtained 
with an accuracy of one part per million. At 
about 16°C a change of one part per million in 
density corresponds to a change of 0.006° in the 
temperature of the buoyancy balance. Of course, 
all the water samples must be carefully distilled 
and freed from dissolved air. 

While we were engaged in the slow process of 
concentrating the oxygen isotopes we realized 
how much more interesting and important it 
would be to obtain the isotope of hydrogen which 
was discovered last year by Urey, Brickwedde 
and Murphy.’ Since the concentration of H? in 
ordinary hydrogen has been recently estimated 
to be only about one part in thirty thousand, or 
less, it seemed hopeless to attempt such a 
separation unless a fractionating process could be 


2 Urey, Brickwedde and Murphy, Phys. Rev. 39, 164 
(1932). 
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obtained with an efficiency of an _ entirely 
different order of magnitude from that of the 
methods that have hitherto been used in the 
separation of isotopes. The more we considered 
the matter the more likely it seemed that such 
processes could be found, and this belief seemed 
to be corroborated by the discovery of Washburn 
and Urey* that the water in an old electrolytic 
cell contains an appreciably larger amount of H? 
than does ordinary water, although they did not 
determine the degree of concentration. 

In this laboratory there is an electrolytic cell 
which has been operated for four years without 
any change except the occasional addition of 
distilled water to replace the water electrolyzed. 
We distilled some of the liquor from this cell and 
found it to have a specific gravity of 1.000034, as 
compared with ordinary water at the same 
_temperature. This very startling increase in 
density, if due to a concentration of H? isotope, 
would indicate the presence of H? to the extent of 
one part in 3000, or ten times the amount which 
had been estimated to be present in ordinary 
water. 

The theory of such a fractionation is simple. 
If the cell runs long enough to approach a 
stationary state, water being continually added 
to maintain a constant volume, then the compo- 
sition of the electrolyte must approach a constant 
value. If a represents the percentage loss of H? 
divided by the percentage loss of H', then the 
ratio of H*® to H! in the cell must ultimately 
approach a value 1/a times as great as the initial 
value, and the composition of the evolved gases 
must approach that of the water added. There- 
fore our first experiment would seem to indicate 
that even if our electrolytic cell has nearly 
reached the equilibrium state, the value of a 
cannot be greater than 0.1. All this assumes that 
the evolved gases contain no water vapor. The 
fact that water vapor is lost in appreciable 
amount still further reduces the possible value 
of a. 

Instead of adding water to keep the electrolyte 
at constant volume we may allow the total 
amount of electrolyte to change. If we call x; the 
amount of H! and x2 the amount of H?, then in 


5 Washburn and Urey, Proc. Nat. Acad. 18, 496 (1932). 
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each stage of the process 

d log x2=ad log x1, (1) 
or, representing by x,° and x2° the initial amounts, 
log x2/x2° =a log x1/x1° Or (2) 


Our next step was to carry out such an elec- 
trolysis. Some of the liquor from the large 
electrolytic cell was placed in a beaker with two 
plates of nickel as electrodes and was electrolyzed 
with a current of about 15 amperes until the 
volume had been reduced to about 2 of the 
initial volume. The electrolyte was then distilled, 
and the increase in density over ordinary water 
was found to be fifty percent greater than 
before. This seemed once more to indicate that 
practically all of the H? was left behind. 

These optimistic conclusions were suddenly 
overthrown by a simple calculation. Professor 
Giauque, who had been using the large elec- 
trolytic cell, was able to make a rough calculation 
of the total amount of water added to the cell 
during its four years of operation. Taking the 
original concentration of H?® to H!' as 1 to 
30,000, it became evident that not more than 
one-third or one-quarter as much H? had ever 
been added to the cell as we were finding there. 
This discrepancy might be explained in one of 
two ways. Either the concentration of H? in 
ordinary water is much greater than had been 
supposed or our increase of density was not 
caused by the concentration of hydrogen isotope 
alone, but by concentration of the heavier 
isotopes of oxygen, which are initially present in 
much larger amount. 

The possibility that any large accumulation of 
heavy oxygen had occurred was soon eliminated. 
Some of the water obtained in the electrolysis 
which we last described, and of which the 
temperature of buoyancy balance was 0.31° 
above that of ordinary water, was distilled in 
such a manner that the steam passed over hot 
iron wool where part of it was converted into 
hydrogen. The gases passed on into a condenser 
from which the water was returned to the still 
while the hydrogen continued into a_ tube 
containing hot copper oxide, where it was once 
more converted into water and was then con- 
densed. After this process was complete a stream 
of hydrogen was passed through the hot iron 
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oxide, producing a second sample of water which 
contained all of the oxygen, while the first sample 
contained all of the hydrogen, of the initial water. 
This second sample had the density of ordinary 
water while the first, which contained only the 
hydrogen of the water under investigation, 
showed again a difference of 0.31° in buoyancy 
balance. 

We were thus forced to the conclusion that the 
concentration of H? in ordinary water is much 
higher than has been usually estimated, ap- 
proaching more nearly the value of 1 to 4500 


first predicted by Birge and Menzel* from 


discrepancies in the atomic weight of hydrogen. 

Our electrolytic experiment can now be ex- 
plained by assuming first that the ratio of H? to 
H! in the hydrogen evolved from the cell is 
higher than we had supposed and second that in 
the calculation of concentration from density, 
correction must be made for an appreciable 
difference between the densities of ordinary 
water and of pure H'H'O. While these new 
conclusions show a poorer efficiency in the 
electrolytic separation, this is partly compen- 
sated by the higher initial concentration of H®. 

We therefore engaged at once in a process 
designed to reduce by electrolysis 10 liters of the 
water from the large electrolytic cell down to 1 cc 
or less. In a large glass tube were placed eight 
parallel nickel plates, each 55 sq. cm in area, 
connected alternately with the positive and 
negative leads, which were also of nickel. A 
copper coil was introduced through which cooling 
water flowed. A current of 250 amperes was used, 
although this was reduced in the later stages of 
the electrolysis, when the volume of the elec- 
trolyte became small. In order to avoid too high a 
concentration of electrolyte at the end of the 
run, we used one liter of 5M alkali from the 
large electrolytic cell, and 9 liters of water 
distilled from the same liquor. 

At the end of five or six days the electrolyte 
was reduced to one liter. Of this, ninety percent 
was placed in a bath at 0°C, while CO. was 
allowed to bubble through it until all the alkali 
was converted to normal NaeCQOs, indigo carmine 
being used as indicator. This solution was then 
distilled from a copper kettle which could be 


‘Birge and Menzel, Phys. Rev. 37, 1669 (1931). 
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heated finally to drive over all the water of 
crystallization. The 900 cc thus obtained, to- 
gether with the 100 cc of alkaline solution 
reserved from the preceding electrolysis, were 
placed in a smaller but similar cell and elec- 
trolyzed to a volume of about 100 cc. This in 
turn, employing the same procedure as before, 
was electrolyzed in a third still smaller cell until 
it was reduced to 10 cc. In the fourth very small 
cell, which, instead of having a cooling coil, was 
surrounded by ice water, the volume of elec- 
trolyte was reduced to a little over half a cubic 
centimeter. This liquid, treated with carbon 
dioxide, was distilled in a vacuum into a small 
container. The water thus obtained was found by 
means of a small picnometer to have the specific 
gravity 1.035; in other words, 31.5 percent of all 
the hydrogen in this water is H’. 

In the next series of experiments changes were 
made in the first and the last electrolytic 
concentrators. In the first, the 10 liter glass 
cylinder was replaced by a 20 liter cylinder of 
monel metal with external cooling, and the 
current was increased to 400 amperes. In the 
last concentrator provision was made for saving 
the evolved gases. These were recombined on 
platinized asbestos, thus giving water containing 
several percent of H®. 

In this second run, concentrating from 20 
liters to 0.5 cc, we obtained water of specific 
gravity 1.073. Accordingly, in this water 65.7 
percent of the hydrogen is H?. While our previous 
experiment showed conclusively that there is no 
large change in the ratio of oxygen isotopes in our 
electrolytic process, and while spectroscopic 
analysis has shown that the hydrogen obtained 
from this sample of water contains over fifty 
percent of H?, two things must be done before we 
can, with certainty, obtain the composition from 
the density. First, we must subject this water to 
the process that we have described, which would 
free it from oxygen in which the isotope ratio 
may have changed, and second, we must prove 
that the density varies linearly with the compo- 
sition. These two experiments have not yet been 
completed. 

In our whole process of concentration we have 
attempted to adhere as nearly as possible to the 
conditions of our large commercial electrolytic 
cell until we could, by theory or experiment, 
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ascertain the best conditions for electrolytic 
separation. Thus all of our electrodes and leads 
were made of pure nickel. However, during the 
course of the experiments stray currents caused 
the transfer of copper from the cooling coils, and 
we have suspected that this copper plating 
may have slightly diminished the efficiency of the 
electrolytic separation, but we have as yet no 
proof of this. We have attempted to keep the 
temperature during electrolysis below 35°C, 
chiefly to avoid loss of water by evaporation, but 
partly also because it seemed likely that the 
efficiency of separation would be higher at low 
temperatures. Here again we have made no 
systematic efforts to prove the correctness of this 
view. 

By determining the density of the water in the 
various stages of concentration we have at- 
tempted to determine the efficiency of the 
electrolytic separation. While there seemed to be 
some individual variations, all the results agree 
pretty closely with the assumption that the 
value of a in Eq. (1) is 0.20. In other words, 
under equivalent conditions, five times as much 
H! as H? is evolved. 

According to this figure, if the water containing 
65.7 percent H? were reduced by electrolysis to 
one-quarter of its volume, it would contain 99 
percent of H®. We have, however, postponed 
making this experiment until larger amounts 
of heavy water become available. 

There is not much that can be said at present 
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regarding the theory of the electrolytic fraction- 
ation. It is evident that under the conditions of 
our experiments there is, between the two 
isotopes, a difference in the electrical polarization 
at the cathode amounting to four centivolts. 
Assuming that this represents chiefly a difference 
in rates, rather than in the true electrode 
potentials, it is hard to guess how much of this 
difference in rates is due to the difference in 
mobility of the two hydrogen ions, to the 
difference in rates of diffusion of the two atomic 
species in the nickel surface, and to a difference in 


rates of recombination of the atoms. These . 


questions are now under investigation. In addi- 
tion Professor Heyrovsky is planning to make 
direct determinations of the overvoltage of the 
two isotopes. 

In conclusion, let us make an estimate of the 
amount of the heavy hydrogen isotope in 
ordinary water, realizing that samples from 
different sources may vary considerably. There 
are two simple ways of making this determi- 
nation. By successive electrolyses, collecting and 
recombining the evolved gases, we may approxi- 
mate to pure H'H'O and find the difference in 
density between it and ordinary water. On the 
other hand, we may concentrate ordinary water 
in an electrolytic cell of known efficiency and 
find the increase in density. Preliminary measure- 
ments by both methods indicate that in Berkeley 
city water there is one part of H® to about 6500 
parts of H'. 
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Theoretical Considerations Concerning the Separation of Isotopes 


VOLUME 1 


Henry EyrinG AND ALBERT SHERMAN, Frick Chemical Laboratory, Princeton University 


The relative importance of effusion processes, the half 
quanta due to van der Waals forces, and half quanta 
associated with homopolar bonds is considered with 
respect to the question of the separation of isotopes. 
Specific reaction rate constants for various isotopic 
molecules in effusion processes have been calculated. It is 
shown that only in the neighborhood of liquid hydrogen 
temperatures or lower are the half quanta due to van der 
Waals forces important. Specific reaction rate constants 
involving the zero point energy of homopolar molecules 


(Received April 8, 1933) 


have been calculated for various isotopic molecules, 
principally those in which a hydrogen atom is one con- 
stituent, at temperatures of 20°K, 90°K, and 273°K. The 
advantage of working at as low a temperature as possible 
is rendered apparent. Even at ordinary temperatures the 
differences in the constants are in general much greater 
than in the case of effusion or van der Waals processes. A 
general equation is derived relating the concentration of 
the light and heavy isotopes, valid for any one of the three 
processes considered. 


ECENTLY theoretical discussions of the 

effect of zero point energy on reaction 
rates, with particular reference to the question of 
separating the hydrogen isotopes have been 
given.! This general question now possesses 
additional interest as a result of recent experi- 
ments.” 

Lewis reports the obtaining of watér of specific 
gravity 1.035, presumably by electrolysis. Taylor, 
Gould and Bleakney obtained a 5-fold increase in 
concentration of H'—H? over that present in 
ordinary hydrogen by desorbing hydrogen from 
an active charcoal surface at liquid air tempera- 
tures. It is the purpose of this paper to consider 
the subject in greater detail and to present the 
results of some calculations which show maxi- 
mum differences in specific reaction rate con- 
stants to be expected, at various temperatures, in 


teactions involving various isotopic molecules. 


We shall first consider effusion processes. The 
rate at which a molecule of molecular weight 
effuses through an orifice is proportional to its 
velocity component in the direction normal to the 
orifice. This in turn is inversely proportional to 


‘Cremer and Polanyi, Zeits. f. physik. Chemie 19B, 443 
(1932); Eyring, Proc. Nat. Acad. Sci. 19, 78 (1933). 

*Washburn and Urey, Proc. Nat. Acad. Sciences 18, 
496 (1932); G. N. Lewis, J. Am. Chem. Soc. 55, 1297 
(1933); H. S. Taylor, Gould and Bleakney, Phys. Rev. 43, 
496 (1933); Washburn, Smith and Frandsen, J. Chem. 
Phys. 1, 288 (1933). 
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the square root of the molecular weight, as is 
shown by kinetic theory. The ratio of the specific 
rates of effusion of gaseous molecules of molecular 
weights M; and Mz, will, therefore, be pro- 
portional to (M2/M,)'. This ratio is of course 
independent of temperature. 

Consider a gas containing two isotopic mole- 
cules of concentration (X) and (Y) effusing 
through an orifice. One can think of the two 
processes as independent unimolecular reactions 
in the ordinary way and obtain for the relative 
concentration of the isotopes the expression 
(Yo) /(Y) = The subscript 0 de- 
notes the initial concentration. In column three 
of Table I is given the relative increase in 
concentration of the heavy isotope to be ex- 
pected in the remaining gas when the pressure is 
reduced to the fraction of the initial pressure 
indicated in column two. 


TABLE I. 

Gas Pressure drop Relative increase in H'—H? 
H'!—H! 1/100 2.28 
H!—H? 1/10 1.52 

1/5 1.34 
1/100 1.15 
1/10 1.07 
1/5 1.05 
N“—N*" 1/100 1.08 
N“—N® 1/10 1.04 
1/5 1.03 
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In the experiments of Taylor, Gould and 
Bleakney an initially adsorbed 5025 cc, desorbed 
to 35 cc, showed a 3-fold increase in the heavy 
isotope, while with 5715 cc adsorbed a reduction 
to 90 cc in a two stage process showed a 5-fold 
increase. It is seen from Table I that the results 
cannot be accounted for by an effusion process 
alone because the 5715 cc to 90 cc reduction 
should produce, at most, a 2-fold increase, in 
contrast to the 5-fold actually obtained. It is 
therefore clear that other factors must play an 
important réle. If the molecules are not initially 
free but are escaping from a solid or liquid 
surface the forces with which they are held will 
further modify the specific rate for the two 
isotopes. To estimate the effect of the van der 
Waals forces in a liquid we shall use the empirical 
potential functions of Lennard Jones.* 

For the potential energy between two hydrogen 
molecules as a function of the distance, 7, between 
them, he gives 


E=(48.9X 10-5) 
— (15.6 10-*) /4r4 ergs/molecule. 


One might object to the use of this function for 
this purpose on the grounds that it was designed 
to explain observed results on the basis of 
classical mechanics but a detailed consideration 
of the potential energy between hydrogen mole- 
cules using perturbation theory gives essentially 
the same results.‘ 

Consider the removal to infinity of a hydrogen 
‘molecule initially at a distance r from a liquid 
hydrogen surface. Assume the molecules of the 
liquid are, on the average, at the corners of a 
cube, the size of which, calculated from the 
observed density of the liquid, is 3.6A. The van 
der Waals attraction of the molecule will be the 
sum of the attractions of the molecules in the 
liquid. The 196 nearest neighbors give the 
energy E presented in Table II, for distances, 
along the normal to the surface, equal to r. The 
contribution to the energy of the remaining 
molecules is small. 

We fit the parabola y—a=}k(x—b)* to the 
values in Table II and use for the vibrational 


-8 Fowler, Statistical Mechanics, Chapter X, Cambridge 
Univ. Press, 1929. 
4 Eyring, J. Am. Chem. Soc. 54, 3191 (1932). 
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frequency w=(1/27)(k/u)', where we take for 
the reduced mass of the oscillator consisting of 
one molecule vibrating against the surface, the 
mass of the molecule. In this way we obtain for 
the oscillators H'H! and H'H? half quanta of 56.3 
and 46.0 calories/mol, respectively. This gives a 
difference of 10.3 calories/mol as the difference 
between the half quanta so that, at the boiling 
point of liquid hydrogen, 20°K, the relative 
specific rates of desorption of H'H! and H'H? due 
to this energy alone is e!?-*/2*29 = 1,29, 


TABLE II. 


r (A) 3.20 3.40 3.60 


—53.96 108 —51.77x 105 


E (ergs/mol) —51.66X10* 


In considering the actual rate of escape from 
the surface this factor should probably also be 
multiplied by the factor arising from the 
difference in velocities for the effusion process 
discussed above. 

Thus we see that the van der Waals forces in 
the evaporation of liquid hydrogen will materi- 
ally increase the concentration of the heavier 
isotope in the residue. 

We now make a rough estimate of the half 
quantum for van der Waals adsorption of a 
hydrogen molecule on charcoal at liquid air 
temperatures. For homopolar molecules the heat 
of dissociation, D, is approximately proportional 
to the force constant k for a wide variety of 
molecules as is shown by the fact that the 
constant a of the Morse equation varies but little 
from molecule to molecule. This same approxi- 
mate proportionality holds for the van der Waals 
potential holding two molecules together. This is 
readily seen by a consideration of the empirical 
curves given by Lennard-Jones. For the evapo- 
ration process of any molecule held to a liquid or 
solid surface by van der Waals forces the force 

- constant will remain about proportional to D 
for different surfaces and different molecules. A 
test of this assumption for the evaporation of 
liquid He, He, Ne, Ne and A indicates that the 
force constant k does not rise quite as rapidly as 
the boiling point. 

In the table & is the force constant for the 
oscillator molecule-surface. These values of - 
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were calculated in the manner described for He, 
using the constants given by Lennard-Jones for 
the equation 


AIT, is the heat of vaporization and 7; the boiling 
point. The value of » used was 143 and of m 
used was 5. 

For the evaporation of any molecule held to a 
liquid or solid surface we may expect Trouton’s 
rule to hold, namely that the heat of vaporization 
is proportional to the absolute temperature at 
which the vapor pressure is atmospheric. Conse- 
quently we expect 3hw=(c/27r)(T/M)? to hold 
roughly for all substances. If c is chosen to fit for 
liquid hydrogen the calculated half quanta for 
substances with higher boiling points will proba- 
bly be too high. As before, / is the reduced mass. 
Thus the difference in the half quanta for H'H! 
and H'H? at liquid hydrogen temperatures is 10.3 
calories while for the same molecules boiling 
from charcoal at liquid air temperatures a 
difference of about 10.3(90/20)!= 21.85 calories 
or less is to be expected for molecules held by van 
der Waals forces. This gives for the ratio of the 
specific reaction rates at 90°K the value 1.13, 
which is less than the difference due to rates of 
effusion. Because of the larger reduced masses 
van der Waals forces will play an even less 
important rdle for substances other than hy- 
drogen. 

It thus seems clear that Urey, Brickwedde and 
Murphy’s® separation of the hydrogen isotopes at 
the triple point (insofar as it departs from an 
equilibrium process) was due to van der Waals 
forces, helped by differences in the effusion rates 
of the isotopes, while the results of Taylor, Gould 
and Bleakney® at liquid air temperatures can 
only be explained on the basis of surface forces 
with greater associated half quanta than those of 
the van der Waals type. 

When the hydrogen molecules instead of being 
held by the long range polarization or van der 
Waals forces are dissociated and form homopolar 
bonds with atoms in the surface the correspond- 
ing half quanta are much larger. In Table III are 


a Brickwedde and Murphy, Phys. Rev. 39, 164 
1932). 
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OF ISOTOPES 


TABLE III. 
Gas He H2 Ne Ne A 
k 34 186 441 238 740 
k/AH, 1.4 0.87 1.1 0.18 0.49 
k/Tp 7.9 9.2 12.9 3.0 8.5 


tabulated values for the extreme case where the 
surface atom and the adsorbed atom behave as 
isolated oscillators. The actual half quanta for 
surfaces will no doubt in general be smaller but 
the table should help in the selection of a 
desirable surface for separation of isotopes by 
either desorption or electrolysis. 

The ratios of the factors e”/*” for light and 
heavy isotopes in columns 3, 4 and 5 give the 
relative specific rates at 20°, 90° and 273°K. The 
advantage of using low temperatures becomes 
strikingly obvious. Thus, for C!®?— H! and C® — H? 
these ratios are 1.310", 5.010?, and 7.7 at 
20, 90 and 273°K respectively. (Table IV.) 

The calculations for electrolysis or desorption 
presuppose reaction processes with identical 
potential barriers for both isotopes. The chief 
difference in relative rates arises from the 
greater half quanta of the lighter isotope. Con- 
sider the case where hydrogen atoms are being 
desorbed from a solid surface to make hydrogen 
molecules. This can occur by two light atoms 
coming off together, or a light and heavy one, or 
finally by two heavy ones combining. If we 
represent the concentration of heavy isotope by 
y, and light ones by x, we have the two simul- 
taneous equations, 
—dx/dt=kx?+ksxy; —dy/dt=ksxy+ks*y*. (1) 
The two equations become identical if we take 
y=cx* so that this relationship is maintained 
between light and heavy isotopes. Evaluating c 
by putting in the initial concentrations x» and yo 
we obtain 


¥/yo= (x/xX0)*. (2) 


Of course s=e~4#/”7 where AE is the difference 
between the zero point energy for an oscillator 
composed of a surface atom and H! and for a 
second oscillator composed of a surface atom and 
H?. Further, s is a proper fraction. 

It is interesting to observe that (2) is the same 
relationship obtained for the separation by 
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TABLE IV. 


Bond 


Bond 3hv(kg.cal.) 


7.2 10% 
.9X 10 
9x 10% 
10% 
10" 
4 K 10" 
5X 10% 


10" 
10% 


XXX 


XXX 


"7X 10% 


DH CO 


1.1 10"! 


Ni—H! .2X10° 
Ni—H? 4X 10° 
.7X 10? 
-1X 10° 
10° 
108 
-7X 104 


wn 
xX X 
3 8 


RO: 


x XX 


bd 
“INIO NN WOH HK KHON OO WWNHO 


effusion or by van der Waals forces if s is inter- 
preted, there also, as the ratio of the specific 
rates for the light and heavy isotopes. Thus the 
determination of the efficacy of separation by the 
three processes reduces to comparing the ratios 
of the specific reaction rates. 

It is of course possible to distinguish experi- 
mentally the first order desorption of molecules 
held by van der Waals forces and the desorption 
of atoms which obey the Eq. (1). Eqs. (1) become 
after the substitution of (2) 


f dx 
zo 


and for yo/xo9 very small this behaves for a 
considerable time as a typical second order 
reaction. 


From Table II we see that the separation of- 


other isotopes than those of hydrogen should at 
least be more readily attained by reactions at 
suitably low temperatures than by a straight 
diffusion process. The particular virtue of elec- 
trolysis is that, even though the activation energy 
(overvoltage) may be considerable, still, by 
applying a suitable potential, the reaction may 


be made to go at low temperatures. The same - 


object may be accomplished by carrying out 
reactions photochemically. With this in mind 
Taylor and Gould are carrying out the hydrogen 


and chlorine reaction at suitably low tempera- 
tures. Of course any other agent than light which 
will supply atoms at low temperatures would do 
equally well. 

Many reactions which in a laboratory would 
proceed too slowly at ordinary temperatures to 
secure the desired separation may have occurred 
geologically where an entirely different time scale 
is available. Low temperature mineralizing and 
erosion processes immediately suggest them- 
selves. 

Since all chemical reactions having a positive 
temperature coefficient will go about half as fast, 
or less, at ordinary temperatures with the heavy 
as with the light hydrogen isotope it will be 
possible to single out such reactions in biclogical 
processes. 

Separation of isotopes due to equilibrium 
differences have been frequently discussed and 
require no additional consideration here. In 
general such separations will be less efficient than 
the corresponding dynamic processes. The heavier 
isotope approaches equilibrium more slowly from 
both directions. 


CONCLUSIONS 


Detailed calculations of the rdle of zero point 
energy in the separation of isotopes by dynamic 
processes are given. Apparently only in the 
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|| 
H!—H! 1.2X10% 94x10 6.5 10! of 
H!—H? 1.2X10%  2.0X108 2.0 10! 

H?—H? 46X10 3.3108 3.6 X 102 tu 

5 1.4104 6.5 X10! 
3 1.5X10° 1.1 108 6.8 X10! 

Ov 3.9X102 1.4104 2.010! 
Ov —H? 3 1.4X10° 1.0K 108 4.3X10! 

—H! 5 3.8X102 1.410 || Hg—H? 1.4X 10! 

—H? 3 1.3X10° 1.0108 || Hgt—H! 2.0 10? 
4 13X10 2.1108 || Hg+—H? 4.210! 
3 2.6X107 2.8X102 || 6.3 X10! 
4 1.2X10% 2.1108 || 6.0X 10! 
3 24X10? 10? || 5.6 X 10! 

Be—H! 3 18X10" 2.5 x10? || 5.3 X 10! 

Be—H? 2 2.4X10° 6.010" || 5.0 10! 
Al—H! 2 6.3X10° 8.210! || 4.8X10° 
Al—H? 1 15X10! 2.410! || N4¥—N™ 4.3X10? 
Fe—H! 1 x10 3.7X10' || 3.9 102 
Fe—H? 1 x10" 1.810% 1.210" |} 44 
Pt —H! 1 1.610! 2.410" |} 4.6X 108 4.3 
Pt —H? 1 x10" 94x10 9.6 Ch — Cl" 3.4108 4.2 
Cu—H! 2 x10 15x10 

Cu—H? 1 3.5 X10! 


SEPARATION 


neighborhood of liquid hydrogen temperatures 
are the half quanta due to van der Waals forces 
important. Processes which involve the breaking 
of homopolar bonds at moderate or low tempera- 
tures are far more efficient than effusion processes. 
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Any agents which provide atoms at low tempera- 
tures which then enter into subsequent reactions 
are particularly effective. Familiar examples are 
catalysts, light, electric discharges, and electroly- 
sis. 
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Molecular Vibrations of Three Particle Systems with Special Applications to the 
Ethyl Halides and Ethyl Alcohol 


Pau C. Cross AnD J. H. VAN VLECK, University of Wisconsin 
(Received April 4, 1933) 


The vibrational potential of polyatomic molecules is 
discussed in the light of the theory of directed valence. 
Andrews’ assumptions of ‘‘valence forces’’ are seen to be 
rather special, but at the same time are fairly good 
approximations under certain conditions. The normal 


- vibrations of three body systems are calculated for the 


potential function 
= +ho(Are)? +r2")(Ay)?] 


with 7, 72 radial coordinates and y an angular one. The 


potential terms proportional to Ar;Ar2, etc., are treated as 
perturbations. These calculations are applied to the ethyl 
halide and ethyl alcohol molecules by considering the 
CH;, CH, and OH groups as dynamic units. Force constants 
were determined which gave approximately the experi- 
mental frequencies. The calculation of the normal modes of 
vibration illustrate the approximate validity of the 
designation of the modes of vibration of unsymmetrical 
three body systems as twe independent radial motions and 
an angular motion. 


I. INTRODUCTION: RELATION OF DIRECTED 
VALENCE TO VIBRATIONAL POTENTIALS. 


N a recent study of the infrared absorption 

spectra of the ethyl halides (CH;CH,CIl, 
CH;CH2Br, CH;CH2I) Cross and Daniels! sug- 
gested for each of the fundamental bands a type 
of motion which was believed to represent the 
chief characteristic of the normal modes of 
vibration involved. Three of these types of 
motion (types F;, F;, Fs of reference 1) may be 
considered as the characteristic modes of a three 
body system in which the units are the halogen 
atom, the CH; group, and the CHe group. These 
were given as illustrated in Fig. 1. The experi- 


7 


(1) (we) (ws) 
Fic. 1. 


mental evidence for each may be summarized 
briefly as: (1), the occurrence of a characteristic 
frequency whenever the group —CH2X is 
present; (2), the occurrence in all the ethyl 
halides of a frequency nearly equal to that of a 
strong Raman line of ethane; and (3), the low 


1 P. C. Cross and F. Daniels, J. Chem. Phys. 1, 48 (1933). 


value of the bending frequency in the ethyl 
halides and the fact that no such low frequency 
is found in the methyl] halides. 

Cross and Daniels also noted the apparent 
absence of the frequency characteristic of the 
H;C—CHg: bond in the case of ethyl alcohol 
(CH;CH2OH). They suggested that the sym- 
metry of the molecule resulted in a strong inter- 
action between the different types of motion in 
such a way that a given frequency could not be 
assigned to a type of motion involving only a 
part of the molecule. (The word symmetry, as 
used here, applies to the equality of the masses 
attached by the central atom. In ethyl alcohol 
the attached masses are 15 and 17 on the atomic 
weight scale.) The same should, incidentally, be 
true for ethyl fluoride. 

The purpose of this paper is to calculate the 
normal modes of vibration of three body systems 
analogous to each of the ethyl halides and to 
ethyl alcohol. The masses, dimensions and 
normal frequencies shall correspond to those of 
the molecules, considering the CH; group and 
the CH: group as units (Fig. 1). It is reasonable 
to assume that the hydrogen atoms, by virtue of 
their small mass and the strong forces holding 
them to the carbon atoms, do not produce any 
first order deviations in the dynamical behavior 
of the molecules from that of three body systems 
in the vibrations we are considering. Hence the — 
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test of the accuracy of the above suggestions. 
Most of the previous work on the small vibra- 
tions of triatomic molecules has dealt with 
symmetrical systems, i.e., systems of the type 
Radakovic® has, however, given a 
general solution for the potential function 


V =} c23(Ares)? +¢12(Arie)? ], (1) 


in which Ar;; is the displacement of the distance 
between the particles 7 and j from the equilibrium 
distance. Instead of (1) we shall use the potential 
function 


+ +12) (Ay)? ]. (2) 


The significance of the coordinates is clear from 
Fig. 2. The coordinate y is an angle rather than 


m,(x+r, cos($ ry), y + ty) 


m, (x,y) 


Fic. 2. 


a distance, while 7; and 7 are identical with 7; 
and ro; of (1). The equilibrium values of 1, re 
are denoted by 1°, 72°. The factor 7:° +72" is 
inserted in (2) to make k; dimensionally similar 
to ki, ke; otherwise no direct significance could 
be attached to comparison of the numerical 
values of the various force constants. 

A system obeying (2) is often called a ‘‘valence 
force system” in the literature. Eq. (2) has a more 
reasonable form of the potential function than 
(1) if the directional valence of the central atom 
is more important than the repulsion between 
the attached atoms. 


*N. Bjerrum, Verh. d. D. Phys. Ges. 16, 737 (1914). 
*D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 
‘R. C. Yates, Phys. Rev. 36, 555 (1930). 
*K. W. F. Kohlrausch, Der Smekal-Raman-Effekt (and 
teferences therein) Julius Springer, Berlin, 1931. 
*M. Radakovic, Wien. Ber. 139, 107 (1930). 
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results of this calculation should be a qualitative 
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The most general form of the potential func- 

tion for small, i.e., harmonic, vibrations would 
include terms proportional to A7,Are, ArAy, 
and Ar,Ay (or terms proportional to Ar;:Arj;, 
etc., in (1)). In writing (1) and (2) we have 
assumed only three atoms, but the extensions to 
molecules with an arbitrary number of atoms are 
obvious, and our discussion on cross-terms in the 
next few paragraphs is not restricted to triatomic 
molecules. 

Actually the terms proportional to Ar;Ay and 
Ar,Ay do not vanish completely. In the first 
place, the theory of directional valence requires 
that the angular equilibrium positions be some- 
what modified when one of the bonds is stretched 
or contracted. This has been demonstrated by 
one of the writers in an article in this journal.’ 
It was there shown, for instance, that in tetra- 
valent carbon compounds, quite irrespective of 
repulsions between the attached atoms, the 
tetrahedral angles 109.5° should be somewhat 
modified when the four bonds are not of equal 
intensity. The bonds must be considered unequal 
not only when the four attached atoms are not 
identical, but also in symmetrical molecules such 
as CX, when an atom X is moved radially from 
its equilibrium position, so that the symmetry 
is lost. Hence the terms in Ar,;Ay and ArsAy 
cannot disappear, for absence of such terms 
implies that the equilibrium angles are unmodi- 
fied by small radial displacements. However, it 
is probably true that the equilibrium angles do 
not vary greatly from 109.5° even when the 
bonds are of unequal intensity or even when the 
central atom is not tetravalent. (Witness, for 
instance, the equilibrium angle in the neighbor- 
hood of 100° in the case H2O. Of course COs is 
collinear and so has an angle of 180°, but here 
there is the complication of double bonds, to 
which our discussion does not apply.) This is a 
very happy circumstance, for otherwise any 
discussion of angular structure in molecules 
would needs be much more difficult than in 
reality. The reason why 109.5° is usually not far 
amiss is the following: for a structure C-X-Y 
involving a single valence angle X-C-—Y the best 
value of the angle was shown by one of us? to 


7J. H. Van Vleck, J. Chem. Phys. 1, 219 (1933); (cf. 
especially Sections 5 and 6). 
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vary from 90° to 180° depending on the relative 
importance of the bonds emanating from the 2s 
and 2p states of the C atom. If the s—? inter- 
play were the only important thing, the best 
angle would be 180°, while 90° would be required 
for pure p-bonding. Actually, the p-bonding is 
perhaps dominant, so that the best angle is 
nearer 90° than 180°. In a compound of the 
form CWXYZ, each angle W-C-X, W-C-Y, 
. . . Y-C~Z has an optimum value between 90° 
and 180° as far as its own particular directional 
valence is concerned. These values are in general 
not compatible, and so must be somewhat modi- 
fied, with the most powerful directional valences, 
of course, in general favored at the expense of 
the others. The important thing, however, is 
that the best value for any individual angle is 
between 90° and 180°, and so never extremely 
different from 110° even when the bonds are of 
unequal intensity. Hence the alteration in 
equilibrium angle upon stretching or contracting 
a bond is probably very small, and so the terms 
in Ar;Ay and Ar2Ay originating from this cause 
are probably subordinate. A term in Ar,Arz is, 
in principle, required because the configuration 
of the C atom is somewhat modified by any 
alterations in position of the attached atoms, so 
that when one bond is stretched, the equilibrium 
radial values for other bonds are altered. Such 
terms, however, are probably very small. Other- 
wise particular bonds would not have nearly 
invariant vibration frequencies in a wide range 
of compounds. 

In the preceding paragraph we have considered 
only the effect of the valences emanating from 
the central atom. A more important cause of the 
“cross-terms’” proportional to Ar;Are, Ar,Ay, 
Ar,Ay is perhaps the repulsions between the 
attached atoms. If ri2 be the distance between 
the attached atoms, this repulsive energy is a 
function of the single argument rj2, and will 
generate cross-terms since 
+cAy }. The cross-terms originating from this 
cause will be computed in the next paper for the 
H,0 molecule. It will be found that these terms 
are sufficiently small to let (2) be the dominant 
part of the potential energy, but are at the same 
time by no means negligible, and are sufficiently 
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great to make purely fortuitous any exact agree- 
ment with experiment achieved by using a 
potential of the form (2). For this reason, it has 
seemed desirable to give in Section V a theory for 
treating the cross-terms as a pe turbation super- 
posed on (2). 

When the specification of the molecule requires 
several angular coordinates yi, ye, °***, the 
vibrational potential will in general contain 
cross-terms proportional to Ay:472, etc., because 
the equilibrium position of any given valence 
angle is modified by angular displacements in 
other bonds. Such cross-terms do not arise in the 
cases treated in this or the following paper, 
because we restrict our attention to effectively 
triatomic molecules, where only one angular 
coordinate is necessary. The well-known assump- 
tions of Andrews*® concerning molecular vibra- 
tions appear to be tantamount to assuming that 
all the various types of cross-terms are neg- 
ligible. Hence his hypotheses are rather re- 
strictive. We have seen that the quantum 
mechanics of directed valence requires that the 
cross-terms, at least in ArAy, are not excessive, 
and this probably explains why Andrews’ theory 
has been attended with some qualitative suc- 
cess. 

After our calculations were completed an 
article by Lechner? appeared, likewise giving the 
derivation of the characteristic frequencies for 
the potential function (2). His paper does not 
give the details of his derivation or analyze the 
perturbing effect of cross-terms and so it seems 
advisable to include the mathematical analysis 
which we used, as well as the applications, which 
are to different molecules than those studied by 
Lechner. Also, we believe the Lagrangian func- 
tion (6) which we give may be useful for other 
three body problems besides the ones which we 
consider. In particular the kinetic energy terms 
of (6) would be useful if one desired to carry 
through a calculation in which the oscillations 
were not infinitesimal, so that the method of 
normal coordinates could not be used. 


8 D. H. Andrews, Phys. Rev. 36, 544 (1930). 
°F, Lechner, Monatshefte f. Chemie 61, 385 (1932). 
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B Il. THE PROBLEM OF THE THREE Bopy SYSTEMS 


The kinetic energy in the Cartesian coordinate system is 


T= 3m cos sin sin cos $}?] 
+ cos —12(d+7) sin }?+ sin cos (6+) }?] 


]. 


(3) 


Let us take the origin of these coordinates at the center of mass, which we may suppose to be 
stationary.’ Then x and y may be eliminated through the relations 


~mixi= 


and the resulting Lagrangian function is 


(4) 


duel + + +1277? | 
cos sin sin y+nre(¢?+¢y) cos 
— ho(Are)?+ ks +172) (Ay)? 


mM2+ Mm 
where 


MM2+ M2M3 


’ 


M2 


m+ M2+ 


Since ¢ does not occur explicitly in L, it is ignorable and ¢ may be eliminated by solving dL /d¢6=P 
for @ and substituting this value of ¢ in the function” R=L—4¢(dL/d¢). Since P is the angular 
momentum of the system, it will be taken as zero inasmuch as, for our purposes, the system may be 


considered devoid of rotation. This gives 


R= — pus?re? sin? y me — pus?rs? sin? y +3 — us” cos? 


+[—us cos y+ pus*rire sin® y sin [uri — cos y cos y 


— +12) (Ay)? ], 


in which 1/p=pury?+ pore? cos 


Since the vibrations are small, the equilibrium 
values 71°, 72°, y° of ri, 72, y, respectively, may 
be substituted in the coefficients of the velocity 
terms in (6). Then if we introduce new coor- 
dinates qa=An, g2= Aro, (ri Ay, we 
have 


(7) 


The values of the constants a;; are immediately 
obtained from Eq. (6). The three equations of 


motion obtained from the Lagrangian function 
(7) are 


Gigi ta 


(aij =4ji). 


(i=1, 2, 3). (8) 


™ Cf., for instance, Whittaker, Analytical Dynamics, 
Third edition, page 54, Cambridge, 1927. 


(6) 


The usual assumption of a periodic solution 
qi=Aie-™ yields the following relations for 
determination of the amplitudes. 


=(, (4 = 2, 3). (9) 


The existence of a non-trivial solution of the 
three Eqs. (9) requires that the determinant of 
the coefficients vanish, whence 


a nr” ky a 12 
ko 193d? =(), 


332 ks 


(10) 


From this equation the constants necessary 
to yield solutions of the Eqs. (9) may be deter- 
mined. 
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III. EVALUATION OF THE FORCE 
CONSTANTS 


The interatomic distances given in Table I 
were calculated from the values of the atomic 
radii given by Sigwick and Bowen.!! The carbon 
to carbon distance 72° was taken throughout as 
1.54A. No experimental values of the angles 


TABLE I. 
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between the valence directions were available, 
but as shown in Section I, they cannot vary 
greatly from the tetrahedral angle. 

Values of the force constants, ki, ke and ks, 
were determined which made the roots of (10) 
equal to the experimental values of \?, which 
have the significance \? = 47°c’w* in terms of the 


7r\9(A) k, ko 
CH;CH.CI 1.74 2.2 3.9 
CH;CH2Br 1.9 1.6 3.8 
“HAC 1.2 3.8 
H; 4.1 3.6 


wi (cm) w2(cm™) 

k; X10 Cale. Obs. Calc. Obs. Cale. Obs. 
0.22 660 655 980 966 6335 
0.22 570 970 960 285 290 
0.18 510 500 960 950 245 200 
0.26 910 885 1020 1045 430 9425 


wave number w.'? In case no reasonable real 
values of the force constants would give the 
required roots of (10), values were selected which 
gave roots nearly the same as those obtained 
from experimental data. Values of these con- 


CH;CH2Cl; os=2.5 X10", 
CH;CHeBr; = 2.25 & 10°, 
CH;CH2I; C13 = 1.96 X 10°, 
CH;CH2OH; c3=4.0 105, 


Insofar as the consistency of the values of the 
force constants is concerned, there is little basis 
for a choice between the two potential functions. 
In particular, the value of ke or C2; should not 
vary greatly from molecule to molecule since the 
the intensity of the CH;—CH2 bond should be 
but little affected by the nature of the other 
atom attached to the C atom. This condition is 
seen to be met with either the potential function 
(2) based on valence forces, or with the less 
likely alternative (1). 


1N, V. Sidgwick and E. J. Bowen, Ann. Rep. Chem. 
Soc. 28, 402 (1931). 

12 Eq. (10) is a cubic, \®*+ad*+6d?+7=0, in d2. In the 
evaluation of its roots A2”, it is convenient to 
utilize the well-known expressions 


stants, together with the calculated values of w 
and the wave numbers of the corresponding 
Raman lines, are given in Table I. 

As an alternative the force constants of the 
potential function (1) were calculated, and are 
as follows: 


X10°, co=1.2 10%. 
93 = 3.75 X 10°, o2=1.0 10%. 
623 =3.7 10°, c2=0.78 X 10°. 
3=3.5 K10°, ce=1.7 K10°.% 


IV. THE NorRMAL MopDEs OF VIBRATION 


After determination of the force constants, ki, 
the roots of Eq. (10) were substituted in (9) to 
obtain the ratio A; : A2 : A; for each frequency. 
The corresponding values of A, were obtained 
from the relation 0L/d0¢=P=0. The resulting 


normal modes of vibration are shown in Figs. 
3 to 6. 


Explicit formulas for a, 8 and y in terms of the masses, 
force constants and molecular dimensions have been given 
by Lechner® and so need not be repeated here. Lechner 
expresses the equations of motion in the form @;= aig 
rather than (8), but we have verified that his secular 
equation is the same as ours. ' 

13 There are no real values of ci, which give the exper!- 
mental values of w for ethyl alcohol. The values given here 
correspond roughly to values 1015, 935 and 425 cm for w. 
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\ (CH,)7 
w, = 660 w2=980 cm! = 335 cm7! 
Fic. 3. 
(CH,)77 
o,=510cm"! w2=960 w; = 245 
Fic. 5. 


(CH,.7 


@,=570 w,=970 cm w3= 285 


Fic. 4. 


Fic. 6. 


#,=910 cm 


V. Tue Errect or Cross TERMS IN THE 
POTENTIAL FUNCTION 


We shall now consider briefly the effect of 
adding the cross-terms 


2 


to the potential function of (7). The existence 
and origin of these terms was discussed in 
Section I. The addition of (11) merely introduces 
into the equations of motion (8) the extra term 
DL ikiq;. The resulting secular equation is still a 
cubic, but as it can be solved only numerically, 
the addition of (11) introduces so many arbitrary 
constants that the problem ceases to be tractable 
or significant. Instead, the extent to which (11) 
influences the frequencies can be more easily 
studied if one introduces it as a perturbation. 
The result of the perturbation is perhaps best 
examined by introducing normal coordinates, 
A, +++, for the unperturbed problem. The coor- 
dinates qi, ge, gs will be linear functions 
gi= > j;a:;A; of these coordinates. After addition 
of (11) these coordinates lose their ‘“‘normal’’ 
property and the Lagrangian function in terms 


(kij =k (11) 


= (0.6212 — 1.1713 — 1.8523) X 10% 


of the A’s may be thrown into the form: 
(12) 
with 


The secular equation is then 


Bie Bis 
Bor Boo + do? — Bog =0. (13) 
Bs1 B32 Baa — 


Here primes are attached to the \’s when we wish 
to designate their values as modified by the 
inclusion of (11). The approximate solution of 
(13) is 


if the off-diagonal elements of (13) are small in 
comparison with the differences \,?—,?. 

The effect of the introduction of the cross- 
terms may best be examined by writing down a 
specific numerical form of (14) for a particular 
example, ethyl alcohol. 


(0.0612? + 1. 1 5ki3? + 0.47 kos” 0.62ky2ki3 +0.35 2.34ki3k23) x 10", 


he!” — = (— 0.9412 X 10% 


(0.0612? + 1.6513? + 1.0623? — 0.54 +0.41 — 2.64 X 


= (0.0312 +0.31h1s+0.37 hes) X 10% 


— (0.01 kis? +0.47 + 0.6423” + 0.082i2k13 + i2k23 — 0.3kiskes) X10. (15) 
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A rough estimate of the k;;’s by a procedure 
similar to that used in evaluating g(ri2) in the 
following paper yields the values ki2=0.13 X 10°, 
ki3=0.039 X 10°, ko3 = 0.044 X 10°. Introduction of 
these terms gives the values \;”=2.9X 10°, 
he” = 3.58 108, =0.685 X for the per- 
turbed frequencies. The unperturbed values are 
= 2.95 1078, Ap? = 3.7 K 1078, A3? = 0.655 X 107°. 
These values could be altered by choosing other 
values for the constants ki, ke, k; in the unper- 
turbed system, but it appears unlikely that the 
separation of \,? and i,” could be made great 
enough to agree with the experimental data. 
The separation is probably partly due to certain 
factors that have been neglected, such as the 
effect of the hydrogen atoms on the potential 
function and on the effective masses of the units, 
or to errors in the dimensions of the molecules. 
Despite this discrepancy, however, the approx- 
imate validity of the calculations is retained. 


VI. CoNcLUSION 
It is apparent from the normal modes of 
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vibration of the systems studied that the desig- 
nation of a type of motion involving only a part 
of the molecule to account for certain frequencies 
is justified under certain conditions of dissym- 
metry. This may in a large measure be respon- 
sible for the well-established fact that certain 
chemical bonds are accompanied by frequencies 
which vary but little from a mean value as 
drastic changes are made in other parts of the 
molecule. However, as was suggested by Cross 
and Daniels, the presence of certain types of 
symmetry, such as the approximate mass sym- 
metry of ethyl alcohol, may cause an interaction 
of the different types of motion in such a way as 
to largely obliterate the characteristic features 
of the separate types of motion. 

The authors are very grateful to the Wisconsin 
Alumni Research Foundation for a financial 
grant to one of us (P.C.C.), and to Professor 
Farrington Daniels for his continued interest in 
the problem. 
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The calculation is based on the Slater-Pauling theory of 
directed valence, but supplemented by data from the 
band spectrum of the diatomic OH molecule. The heat of 
dissociation and valence angle, obtained incidentally, are 
in good agreement with experiment. The “‘cross-terms’”’ in 
the vibrational potential prove to be of subordinate 
importance. It is found that two of the fundamentals 
should fall near 2.84, and that they should almost coincide. 
This designation of the observed wide absorption band at 


2.6-2.8u as the superposition of two fundamentals agrees 
with Mecke’s interpretation of the H:O spectrum, but 
contradicts Johnston and Walker’s. The computed value 
1660 cm-! of the third fundamental frequency »; differs 
from experimental determinations by only 4 percent, but 
any closer agreement than 15 percent is accidental in view 
of the approximate nature of the theory. Since w; represents 
a primarily angular vibration, its calculation is a particular 
test of the concepts of directed valence. 


HE water molecule is of interest because it 

is possible to calculate its vibration fre- 
quencies and other constants by a method based 
on quantum mechanics supplemented by spectro- 
scopic data from the OH and Hg (but not the 
HO) molecules. This process should not be 
confused with the less exacting, more customary 
procedure, used in the preceding paper, of 


simply adjusting force constants to yield agree- 
ment with observed vibration bands. Our calcu- 
lation of the fundamental frequencies will, of 
course, include a theoretical evaluation of the 
coefficients in the vibrational potential, 


(Ay)? 


(1) 


where 71, 72 are the O—H separations, and y is 
the angle H—O—H (cf. Fig. 2 of reference 1). The 
factors 7,°, 7;° are inserted, as in reference 1, 
to secure dimensional similarity in the k’s. It 
turns out that the coefficients ki2 and ki; are 
relatively unimportant. This fact furnishes at 
least partial justification for treatment of the 
“cross-terms” in the potential as small per- 
turbations, as in the preceding paper on the 
ethyl halides,’ or even the complete omission of 
these terms, as in the hypothesis of ‘valence 
forces” often made in the literature. 


*P. C. Cross and J. H. Van Vleck, J. Chem. Phys. 1, 
350 (1933). See this reference for more general discussion of 
molecular vibrations. 


According to the quantum theory of directed 
valence, the potential energy of an H2O molecule 
should be approximately of the form 


V=f(ri) +f(r2) +2 

+3(Neo— sin? (3y—47). (2) 
Here N,, and N,, are, except for sign, exchange 
integrals between a 1s state of H and respectively 
the 20 and 2p7 states of O, space quantized 
relative to the corresponding O—H axis.” The 
terms f(r), f(re) in (2) are radial O— H potentials, 
and, to our approximation, are of the same func- 
tional form as the potential energy curve for the 
diatomic OH molecule. The term g(ri2) repre- 
sents the repulsion between the two H atoms, 
whose separation is denoted by 712. The remain- 
ing term gives the effect of directional valence, 
and proves to be considerably more important 
than g(ri2). That this term has the form of 
dependence on y given in (2) follows from the 
Slater-Pauling theory of directed valence pro- 
vided one makes the hypothesis of electron 
pairing.* 


2 The notation is chosen in conformity with that used by 
one of us in J. Chem. Phys. 1, 219 (1933). The reader is 
referred especially to Sections 5 and 6 of this reference if 
he desires further information concerning the concepts of 
the quantum mechanics of directed valence, electron 
pairing, etc. which we are using. Note especially that our 
N’s are the negative of exchange integrals. 

3 The angular term of (2) is essentially the same as that 
given by Slater in Phys. Rev. 38, 1139 (1931), except that 
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Details of the derivation of the last term of (2) are as 
follows: Let us take the O nucleus as the origin and suppose 
that the H nuclei lie in the xy plane. Let us further take the 
bisector of the angle H—O—H to be inclined at 45° to the 
x axis. Then 


cos x) =cos (r2, y)=cos ir). (3) 


The configuration of the O atom may, with this choice of 
axes, be taken‘ to be 15°2s*2p072p0,2po,. Here 2p0., for 
instance, denotes a 2p wave function having m;=0 when 
space quantized relative to the z axis. From the rotational 
transformation properties of 2p wave functions, it follows® 
that the exchange integral J,, connecting the 1s wave 
function of H atom 1 and the 20, wave function of O is 


Jiz= — Noo 008? (11, x) — Nex sin? (1, x)=Jey. (4) 


Similarly 
— Nyx, (5) 


Jez = Jiy= — Noo (11, x) — Nyx cos? x). (6) 


We shall let J denote the sum of the exchange integrals 
for the pairs (1s of H; 1s of O) and (1s of H; 2s of O). 
Clearly J has the same value for either H atom and is inde- 
pendent of y. If we make the hypothesis of electron pairing 
the exchange energy is® 


+ Joy — S22 — tSiz— J). (7) 


The terms with coefficients +1 in (7) are those which 
represent pairs. When one expresses the exchange integrals 
involved in (7) in terms of Noo, Nex, y with the aid of (3-6), 
one obtains the last term of (2) except for an additive 
constant —2J—2Nge+3Nrzx. This constant must not be 
inserted in (2), since it is included in f(7:)+/(r2). Namely, 
the binding energy for OH is 


= —J— Noo and Coulomb terms, (8) 


so that f(r:)+/(r2) absorbs all of the exchange energy at 
y=}. Eq. (8) is proved from (7) by dropping all terms 
which involve atom 2, and keeping atom 1 on the x axis. 
The polar and Coulomb terms are to be added in (8) since 
(7) aims to give only the exchange energy. We have tacitly 
assumed that the polar and Coulomb effects do not depend 
on direction. This appears to be an allowable approxi- 
mation. 

It must be cautioned that in the foregoing, we have 


he has introduced explicit numerical values of the N’s. In 
fact our calculation of the equilibrium position and energy 
is very similar to that given by Slater except that we use 
different, and we believe better, numerical constants. Our 
equilibrium calculations are intended primarily as pre- 
requisites to the vibrational analysis, which was not 
considered by Slater. 

* By the methods of reference 2 one can show that any 
other, less symmetrical, choice of spatial quantization for 
the 2p wave functions of oxygen will yield a higher value of 
the energy if the hypothesis of electron pairing is made. 

5 See Eq. (6) of reference 2 and related discussion. 

§ Cf. Eq. (35) of reference 2. 
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assumed that the O atom always remains in the con- 
figuration s*p*, or in other words that there is no s—p 
hybridization for this atom. Any amount of this hybridiza- 
tion tends to intensify the bonds, but the s— p separation 
of the O atom amounts to over 15 volts, and so the O 
atom is probably unwilling to make a partial change from 
s°p* towards sp* in order to hold the H atoms tighter. On 
the other hand, hybridization apparently must be con- 
sidered in connection with the molecular vibrations of 
NHs, as otherwise it seems impossible to explain the low 
potential hill’ (4 volt) involved in converting one isomer of 
NH; into the other. This difference between H2O and NH; 
does not appear unreasonable, since the s— separation is 
smaller (11 volts) in N than in O and since hybridization is 
a greater advantage with 3 than with 2 bonds. 


We shall assume that f(r:), f(r2) and g(ri2) 
can be represented by Morse functions. This 
procedure has been at least partially justified by 
the success of Eyring’s® well-known investiga- 
tions of activation. The values of f(7:) and 
g(ri2) can then be determined from the potential 
energy curves yielded by the band spectra of OH 
and H: and are as follows? 


= 120(1 — 120, 
=— 35(1 — 4 35, 


Here distances are supposed measured in Ang- 
stroms, and the amplitudes have been in terms 


of a kg-cal./g-mol scale of energy. In evaluating 


7D. M. Dennison and G. E. Uhlenbeck, Phys. Rev. 41 
313 (1932); N. Rosen and P. Morse, ibid. 42, 211 (1933): 

8 H. Eyring and M. Polanyi, Zeits. f. physik. Chemie 12, 
279 (1931); H. Eyring, J. Am. Chem. Soc. 53, 2537 
(1931); Sherman and Eyring, ibid. 54, 2661 (1932). 

® Our choice of Morse functions has been determined to 
fit the vibration frequency and the heat of dissociation 
rather than the former and the anharmonic correction. For 
the band spectrum data see Weizel, Bandspektren, or 
Mulliken, Rev. Mod. Phys. 4, 1 (1933). In our formula for 
g(ri2) we have not attempted to make allowance for the 
fact that the polar character of HxO makes the H atoms 
each have an excess of positive over negative charge and so 
exert ionic repulsive forces upon each other. This repulsion 
is perhaps roughly counterbalanced by the van der Waals 
force between the H atoms, not included at large separa- 
tions in the Morse function. Both the ionic and van der 
Waals forces are probably subordinate to those covered 
by the Morse function. The van der Waals term is hard to 
estimate accurately, as Eyring notes (J. Am. Chem. Soc. 
54, 3201 (1932)) that the inverse sixth power polarization 
potential gives excessive values at distances as small as 
ria. In this reference he compares qualitatively the forces 
resulting from directional valence and H—H repulsions in 
H.O, and notes that the latter are, as we find, not entirely 
negligible compared with the former. 
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g(ri2), we have, like Eyring, supposed that the 
binding energy of an He molecule is 10 percent 
Coulomb, and 90 percent exchange. This makes 
g(ri2) 35 percent of a Morse function for He, 
inasmuch as exchange terms connecting electrons 
involved in different pairs involve the coefficient 
—} rather than +1, making* the proportionality 
factor 0.1— (3 X0.9) = —0.35. On the other hand 
f(r;) and f(r2) are complete Morse functions for 
OH. 

The values of NV,, and N,, will be taken as 2 
and —0.6 volt, respectively. These estimates 
have been obtained largely from the potential 
curves of the OH molecule, and are perhaps 
accurate to 20 percent or so. The details of the 
method of estimation will be published in a 
future paper by J. R. Stehn and J. H. Van Vleck. 
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Incidentally, the agreement with the values 2.3 
and —0.6 volts obtained by Coolidge’® by actual 
computation of the integrals with Zener’s wave 
functions is gratifying. We believe the spectro- 
scopic method to be perhaps more reliable than 
the computational one. As a matter of fact, the 
spectroscopic method gives a value of N,, 
somewhat below 2 volts, but we have used the 
round value 2 in order to give some weight to 
Coolidge’s calculations. With our estimate of 
the N’s, one has 


3(Noe— Nex) = 180 kg: cal./g- mol. 

One must now make a Taylor’s development 
of (2) about the equilibrium position r,°, 72°, y°. 
Inasmuch as = 7;°>+ 72? —2rrecos y and 7;°=r2° 
this development is 


V=V(n°, 72°, =) (Ar; +Are) — Nex) cos Ay 


Orie 


+3(Noo— Nex) sin 7: 


sin? 
Orie }y° 


The values of 7;°, 72° and y° are respectively 
determined by the requirement that the coef- 
ficients of Ar;, Arz and Ay vanish in (9). This 
gives 


(An — 


Angstrom, y°®=100°. (10) 


This value of r° agrees with the fact that the 
experimental determinations are practically the 
same as the equilibrium distance 0.98 Angstrom 


"A. S. Coolidge, Phys. Rev. 42, 189 (1932). 

"For instance, Eucken finds y°=110° 56’ while Mecke 
obtains a value between 102° and 111°. 

The good agreement of the calculated and observed 
heats of dissociation should not be considered any par- 
ticular triumph, for any theory based on electron pairing 
will predict approximately twice as great a heat of dis- 
Sociation for HsO as for OH. (Cf. J. H. Van Vleck, J. 
Chem. Phys. 1, 177 (1933). If we assume electron pair- 
ing, it is not necessary to neglect interaction between the 
H atom and the 2s or 2px states of O, as there stated, in 


=) [sin (Ary+Are) +171" cos Ay 


Ars)?-+} cos sin 34?- (6) 


in the diatomic OH molecule. The value of y° 
is in sufficient accord with the various estimates 
of from 102° to 111° which have been obtained 
from band spectra.!! The heat of dissociation 
into 2H +O is V(r, 2°, 7, 112°), or 205 
kg-cal./mol. Expressed in volts, this estimate is 
8.9 volts. Thermochemical data lead to a value 
of about 10 volts.'? Our theoretical estimate is 


order to obtain the factor 2.) For instance, Slater* secures 
equally good agreement with entirely different assumed 
values of the exchange integrals. On the other hand, the 
calculation of the transverse vibration frequency w; is 
sensitive to the evaluation of these integrals. Our semi- 
quantum mechanical, semi-spectroscopic calculation of the 
heat of dissociation should not be confused with Coolidge’s 
purely quantum mechanical treatment, which led to the 
excessively low value —5.7 volts. The writers believe that 
two reasons why Coolidge obtained too low a result may be 
the following: (a) He assumed the same exponential factor 
for the carbon 2s and 2 wave functions, whereas actually 
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raised to 9.4 when correction™ is made for the 
L, S structure of the O configuration s*p*. 

By substitution of (10) in (9), we find that the 
coefficients in (1) have the following numerical 
values: 


ki =6.96 X 10°, 
ki2=0.165 X 10°, 


ks =0.379 X 105, 
his = 0.04 X 108, 


The values of the constants involved in (7) of 
our preceding paper are 


ay 15.6 10-*, 
dy. = 0.185 10-*, 


33 = 3.86 X 10-5, 
143 = A23 = 0.32 x 10-*, 


Because of the addition of the cross-terms, the 
simultaneous equations (8) of the preceding 
article for determining the molecular vibrations 
are to be enhanced by the term > jk:;g;, with 
corresponding modification of the cubic secular 
equation (10) of the foregoing paper. In virtue 
of the symmetry of the molecule, this cubic 
factors into a linear and a quadratic equation in 
47°c*?w?; as follows: 


(di — — ki =0, 
ki — Riz — | 


oF 2[aisd? — kis P =(, (1 1) 


Substitution of numerical values yields 
= 3520 cm, we = 3560 cm™"!, w3 = 1660 (12) 


Experimentally regions of intense absorption 


the 2s states are more firmly held and so have a larger 
negative exponent. This implies that Coolidge over- 
estimated the exchange repulsion from the 2s states. (b) 
Examination of potential curves for OH suggests that the 
ionic attraction per atom is somewhat larger than his 
estimate of 1.1 volts. Coolidge’s difficult and noteworthy 
calculation goes to the refinement of considering more than 
simple permutations. The results of his paper suggest that 
possibly the first order Heitler-London calculation is 
somehow a better approximation without than with 
allowance for the higher order permutations. 

13 Under the hypothesis of electron pairing the valence 
state of H,O may be shown tobe }*P+}'D. The fractional 
coefficients occur because the pairing prevents L, S of O 
from being good quantum numbers. When we double the 
Morse function for OH we allow the O atom twice the 
activation energy necessary to raise it from *P to the 
valence state. Actually the activation energy of O is the 
same for OH as for H,O, so that our estimate of the heat of 
dissociation should be increased by one quarter of the 
interval 'D—§4P, or by 0.5 volt. 
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occurs at 1600 and 3700 cm." It is not at 
all surprising that the calculated radial vibrations 
come near 3600 cm for this is merely because 
the primarily radial vibrations have nearly the 
same frequency as that characteristic of OH. 
Conclusions of more interest are the following: 

(1) The coefficients ki2, ki3 in (1) are com- 
paratively small, so that inclusion of the ‘‘cross- 
terms’ does not greatly modify the calculated 
frequencies. In fact, if one sets kiz=k13;=0, the 
computed frequencies become 3560, 3520, 1660 
cm instead of (12). The frequencies w; and w» 
have the significance of being associated with 
unsymmetrical and symmetrical primarily radial 
vibrations: i.e., w; has Ar;= —Are, Ay=0 while 
like W3, has Ar,=Aro. 

(2) The difference \;?—),’ is given approx- 
imately by the equation 


2(—anki 


ay" — 12" 


2[ais*(ki t+ Riz) — 


a33(du 


A? d2? = 


(13) 


The second term is quite small in the case of 
water, contributing only about 10 percent of the 
total separation. It is interesting to note that 
the contribution of the Riz term is about twice 
as large as that of the aj. term, and is of the 
opposite sign. The value of ai2 is the following 
function of y°: 


sin? 


2(u—s cos y°) 


(14) 


The second term of (14) is about 3 as large as 
the first one. Attempts have been made in the 
literature to calculate y® for H2O from the exper- 
imental values of \,? and Az? by means of the 
relation 


Mut+Mo 
Mu 


cos 7° = — (15) 


obtained by setting cos y° and 
=kie=kis=0 in Eq. (11). It is obvious that 
such a method of obtaining the equilibrium angle 
is not reliable since small values of the cross- 


144 W. W. Sleator, Astrophys. J. 48, 125 (1918). 
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terms ki. will completely alter the character of 
the right side of (13)." 

(3) The difference w;—w2 is so small that it is 
comparable with the rotational structure. Hence 
it is not surprising that the frequencies w1, we 
have not yet been resolved experimentally. Ac- 
tually, the band at 3700 cm™ is extremely 
diffuse. Johnston and Walker" claim to have 
located we at 984 cm~ from the Raman effect, 
but such a value of we is entirely out of the ques- 
tion unless the theory is devoid of all meaning. 
Our interpretation of the band as a superposition 
of w1, w2 is similar to that proposed by Mecke.'” '* 

(4) The HOH? molecule, involving the isotope 
H?, may be assumed to have the same potential 
function as the HOH molecule. If the force 
constants ki, k; are adjusted to give absorption 


'*In molecules where the central atom is not of large 
mass compared to the attached ones, the difference 
\2—2? is not small, and must be calculated from the exact 
formulas (11) rather than from the approximate relation 
(13). When the central atom and attached masses are 
comparable, as in the examples studied by Lechner, 
(Monats. f. Chem. 61, 395 (1932)) the terms involving a2 
overshadow those proportional to ki: and 3. Hence 
Lechner’s calculation of equilibrium angles would not be 
materially modified by inclusion of cross-terms in the 
potential. 

6H. L. Johnston and M. K. Walker, Phys. Rev. 39, 535 
(1932). Our «1, we, w; correspond to Johnston and Walker's 
Dennison’s) v3, “1, v2, respectively. 

' R. Mecke, Zeits. f. physik. Chemie 16, 431 (1932). 

'S Note added in proof: Papers by Mecke, Baumann and 
Freudenberg (Zeits. f. Physik 81, 313, 445, 465 (1933)) have 
just appeared which give a very careful analysis of the rota- 
tion-vibration bands of water. They find w; and a» at 3756 
and 3600 cm~, respectively. Thus w2—w;, according to their 
analysis, has a different sign and larger numerical magni- 
tude than in our theoretical estimate w2—w; =40 cm. This, 
however, need not cause concern, since w:—.« is still small, 
in qualitative accord with theory, and quantitative accuracy 
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regions at approximately 3750 cm for w; and we 
and at 1600 cm™ for w; of the HOH molecule, 
the fundamental bands of HOH? should occur 
at 3750, 2720, and 1400 cm—. 

(5) The agreement between the calculated and 
observed values of w; is particularly gratifying, 
for ws; represents an angular mode of vibration 
not found in OH, and originates in the resistance 
of directional valence to any bending of the 
molecule from the equilibrium angle. It must be 
remembered that the calculations, from their 
very nature, do not claim any great precision, 
and agreement only to within 15 percent'® or so 
would not have been disconcerting. The calcu- 
lations on w; thus confirm the essential correct- 
ness of the Slater-Pauling concept of directed 
valence. 


in the theoretical calculation of a difference as small as 
w:—w) is not to be expected in view of the many approxima- 
tions which it is necessary to make, notably the assumption 
of electron pairing. Also, the experimental determination of 
w2 is not free from difficulty, and Mecke mentions that the 
region from 3300 to 3600 cm has so far defied exact 
analysis. An interesting problem is presented by the fact 
that Mecke finds w». to be an inactive frequency except when 
involved in combinations with wo. If one assumes that the 
motions of the water molecule can be represented by the 
vibrations of a system of three particles, each with a con- 
stant charge, then one must expect w: to be an active fre- 
quency. However, it is probable that the water molecule 
becomes increasingly ionic in character when 73 and 723 are 
simultaneously diminished. If so, migrations of charge 
back and forth between the O and H atoms accompany the 
vibrations. It is thus conceivable that equal displacements 
in f3 and 723 at nearly constant y may give rise to little 
change in moment, due to a cancellation of the change due 
to migration and the change calculated for a model in 
which the atoms have constant charges. In this fashion w2 
might be inactive. 

19The error in w; can be less than in the exchange 
integrals because w; is roughly proportional to (Nos— Nez). 
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The Nature of the Chemical Bond. V. The Quantum-Mechanical Calculation of the 
Resonance Energy of Benzene and Naphthalene and the Hydrocarbon 
Free Radicals’ 


Linus PAULING AND G. W. WHELAND, Gates Chemical Laboratory, California Institute of Technology 
(Received March 21, 1933) 


The secular equations corresponding to the five canonical 
structures for benzene and the forty-two for naphthalene, 
considered as six and ten-electron systems, respectively, are 
set up and solved with certain simplifying assumptions, 
leading to energy values differing by 1.1055a and 2.0153a, 
respectively, from those corresponding to unexcited 
(Kekulé-type) structures, a being a single exchange 
integral involving neighboring carbon atoms. Equating 
these values to the empirical values of the resonance 
energy, « is found to be about —1.5 v.e. 

It is pointed out that the dissociation of certain substi- 
tuted ethanes into free radicals is due not to weakness of the 
carbon-carbon bond in the ethane but to the stabilization of 
the free radicals resulting from resonance among the 
structures in which the unpaired electron is located on the 


methyl carbon and those in which it is on other atoms 
(ortho, para, etc., to the methyl carbon). The secular 
equations for a number of such radicals have been solved, 
neglecting excited structures. The experimentally de- 
termined heat of formation of hexaphenylethane from 
triphenylmethyl, 0.5 v.e., when equated to the calculated 
value C—C+2.2156a, with C—C=3.65 v.e., leads to 
a=-—1.4 v.e. The calculated tendencies towards dissoci- 
ation are in satisfactory agreement with observation, such 
features as the smaller dissociating power of @-naphthyl 
than of a-naphthyl and of biphenylene than of diphenyl 
being accounted for, so that resonance among the structures 
considered may be accepted as the principal effect causing 
the stability of the hydrocarbon free radicals. 


THE STRUCTURE OF BENZENE AND NAPHTHALENE 


NUMBER of structural formulae have been 
proposed for benzene, but none of them is 
free from very serious objections. 


Claus 


The oldest and best known structure is that 
proposed by Kekulé (I). The objections to it are 
twofold. First, it suggests that ortho-disubstituted 
derivatives should exist in two isomeric forms—a 
phenomenon which has never been observed. 
Kekulé avoided this difficulty, however, by 
assuming that the double bonds were in a state of 
constant oscillation such that any two adjacent 
carbon atoms were connected part of the time by 
a single bond and part of the time by a double 


bond. For example: 


1 The fourth paper of this series appeared in the J. Am. 
Chem. Soc. 54, 3570 (1932). 


-X 


Second, the Kekulé structure suggests that ben- 
zene should be a highly unsaturated and com- 
paratively unstable compound, in complete 
contradiction to the observed facts. Since this 
point has been discussed in considerable detail in 
various places, we shall not go into it here. 
The two centric structures (II and III) were 
proposed to avoid both of the above difficulties. 
They are indeed successful in eliminating the 
necessity for the oscillation hypothesis, but they 
can hardly account for the stability of the 
molecule. In the Claus structure the diagonal 
bonds would be very weak, as a result of the 
large distance between the atoms in the para 
positions, and the structure would probably 
represent a molecule less stable even than the 
Kekulé structure. In the Armstrong-Baeyer 
structure the meaning of the six lines pointing 
toward the center is not at all clear. If these 
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“central bonds” are left undefined, the structure 
is essentially meaningless; if they are considered 
to represent free valences, the molecule should be 
highly unsaturated. 

The Dewar structure (IV) has never been 
seriously considered since it suffers from the 
same defects, in an accentuated form, as the 
Kekulé structure. Its advantage over the other 
structures is largely that it explains the intimate 
relationship that usually obtains between the 
para positions in the benzene ring. 

The Ladenburg prism structure (V) can be 
definitely ruled out, since we know that the 
benzene ring lies in a plane, or very nearly in a 
plane. The chemical evidence is also unfavorable 
in this case. 

This by no means exhausts the list of structures 
proposed for benzene. The remainder, however, 
such as the structure of Thiele based upon his 
theory of partial valence and that of Collie based 
upon a dynamic model, are usually complicated 
and cannot be represented in any simple valence 
scheme. 

In the case of the aromatic hydrocarbons with 
condensed ring systems the state of affairs is 


quite similar. Thus the following structures have _ 


been proposed for naphthalene: 


for anthracene: 


and so on for the other aromatic hydrocarbons. 
As in the case of benzene, however, the proposed 
structures are all unsatisfactory, on account both 
of the non-occurrence of predicted types of 
isomerism and of the unexpected stability of the 
molecules. 

Apparently here we have a case where the 
classical ideas of structural organic chemistry are 
inadequate to account for the observed properties 
of a considerable group of compounds. With the 


development of the quantum mechanics and its 
application to problems of valence and molecular 
structure, it became evident to workers? in this 
field that the resonance of benzene between 
the two equivalent Kekulé structures was an 
essential feature of the structure of this molecule, 
accounting for the hexagonal symmetry of the 
ring and for its remarkable stability; and it 
seemed probable that the quantum mechanical 
treatment of aromatic molecules would lead to a 
completely satisfactory explanation of their 
existence and characteristic properties. E. Hiickel* 
has made a valuable start in this direction in a 
series of papers on the quantum mechanics of 
benzene. His method of attack, however, is very 
cumbersome. In this paper we present a treat- 
ment of the problem which is rather closely 
similar to that of Hiickel’s and which leads to the 
same result in the case of benzene, but in which 
the calculations are simplified to such an extent 
that the method can be extended to the naph- 
thalene molecule without undue labor. Further- 
more, with the aid of additional simplifications 
and approximations we have been able to treat 
the problem of free radicals and to obtain results 
in surprisingly good qualitative agreement with 
experiment. We shall leave the discussion of these 
latter calculations to the second part of the 
paper. 

We shall set up the problem in essentially the 
same way as Hiickel. We assume that each of the 
six carbon atoms possesses two K electrons, and 
four L electrons, one for each of four orbital 
functions formed by linear combination of the 
2s and the three 2p orbits. Three of these L 
orbits, each a combination of the 2s orbit and the 
two 2? orbits in the plane of the ring, form single 
bonds to the attached hydrogen atom and the 
two adjacent carbon atoms. The fourth orbit for 
each atom remains a pure p-orbit, projecting at 
right angles to the ring. We neglect the energy of 
the electrons forming the system of single bonds 
in the plane and of their interaction with the 
electrons occupying the pure p-orbits (these 
energy quantities occurring in the same way for 
all the structures considered, and hence leading 
to only a change in the arbitrarily-chosen zero of 

? As, for example, J. C. Slater, Phys. Rev. 37, 489 (1931). 


3 E. Hiickel, Zeits. f. Physik 70, 204; 72, 310 (1931); 76, 
628 (1932). 
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energy), and consider only the interaction energy 
of the latter electrons, which may interact with 
one another in different ways. That is, we treat 
benzene simply as a six-electron system with spin 
degeneracy only, and naphthalene as a ten- 
electron system with spin degeneracy only. 

Considerable justification for this choice of 
orbital functions and of bond distribution is 
provided by the fact that each of the single 
exchange integrals between a pure p-orbit and an 
orbit in the plane of the ring is positive (arising 
from the e?/r;; term in the Hamiltonian only), 
and the chosen distribution of bonds causes these 
integrals to occur with the negative sign in the 
expression for the energy of the molecule, while 
the bonding exchange integrals, which are 
negative, occur with the positive sign. 

In carrying out the calculations we make 
certain further simplifying assumptions. We 
neglect all exchange integrals of unity, and all 
exchange integrals of the energy H except single 
exchange integrals involving two adjacent atoms. 
The single exchange integrals involving adjacent 
atoms, (abcdef/ H/bacdef), (abcdef /H/achdef), etc., 
are represented by the symbol a. These integrals 
are seen to be equal in benzene; in naphthalene 
and other aromatic molecules, in which they are 
not all required to be equal by the symmetry of 
the molecule, the reasonable assumption is made 
that no serious error is introduced by giving them 
all the same value. The Coulomb integral 
(abcdef/H/abcdef) is represented by the symbol 


Q. 


i(0—W)+3a/2 

(Q—W)+3a/2 
2(Q—W)+3a/2 
2(Q—W)+3a/2 
2(Q—W)+3a/2 


(Q—W)+3a/2 
i(Q— W)+3a/2 
3(Q—W)+3a/2 
3(Q—W)+3a/2 
3(Q—W)+3a/2 


This quintic equation is easily reduced to three 
linear factors and one quadratic factor, the roots 
being —2a, —2a, 0, (— (13)!—1)a, and ((13)!—1)a. 
Since @ is negative, the last of these roots, 
((13)'—1)a=2.6055a, represents the normal 
state of the molecule. The eigenfunction corre- 
sponding to this is (before normalizing) ~y=Wa 


}(Q-W) +3a/2 
(Q-W)+3a/2 
(Q-W) 

4(O—W)+3a/2 
1(Q—W) +3a/2 
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Benzene 

The benzene molecule can now be treated very 
simply by the Slater method, with the help of the 
rules formulated by one of us‘ for finding the 
matrix elements occurring in the secular equa- 
tion. The bonds between the six eigenfunctions 
can be drawn so as to give the independent 
canonical structures shown in Fig. 1. Any other 


Fic. 1. The five canonical structures contributing to the 
normal state of the benzene molecule. 


structures that can be drawn can be represented 
by eigenfunctions which are linear combinations 
of the five above. It will be seen that $8tructures A 
and B are simply the two Kekulé structures, and 
C, D and E are three different forms of the 
Dewar structure. The Claus centric structure 
does not belong to the canonical set, but it can be 
represented as a linear combination of canonical 
structures; and the same is true of the Ladenburg 
structure if it be considered spread out in a plane, 
as: 


(The Claus structure =A+B—C—D-—E; the 
Ladenburg structure =A+B—D-—E.) 
As can be easily verified, the secular equation is 


2(Q—W)+3a/2 
2(Q—W)+3a/2 
1(Q—W)+3a/2 
(Q—W) i(Q—W)+3a/2 
(Q—W) 
yz). If we had neglected 
the resonance phenomenon and calculated the 
energy for one of the structures A or B, we would 
have obtained the value W’ = Q+1.5a. Hence the 


extra energy of the molecule resulting from 
resonance among the five independent structures 


‘L. Pauling, J. Chem. Phys. 1, 280 (1933). 


2(0—W)+3a/2 
2(Q—W)+3a/2 
i(0—W)+3a/2 


Ay 
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is 1.1055a. It is interesting to see how much of 
this extra energy is due to resonance between the 
two Kekulé structures and how much is con- 
tributed by the excited structures C, D and E. A 
simple calculation shows that 0.9a@ or approxi- 
mately 80 percent of the resonance energy comes 
from the Kekulé structures alone and only about 
20 percent from the three excited structures. 

In the following paper of this series® a value of 
about 1.7 v.e. has been found from thermo- 
chemical data for the resonance energy of 
benzene. Equating the negative of this quantity 
to 1.1055a, we calculate the value of a to be 
about —1.5 v.e. This value may not be very 
reliable, however, since it is based on the 
assumption that values of bond energies obtained 
from aliphatic compounds can be applied directly 
to aromatic compounds. 

The results of the calculation for benzene are 
summarized in Table I. They are identical with 
those obtained by Hiickel. 

The conclusions we draw regarding the struc- 
ture of the normal benzene molecule are the 
following. The principal contributions to the 
structure are made by the two Kekulé structures, 
resonance between them stabilizing the molecule 
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TABLE I. 


Total Resonance 
energy energy a:b 


Single Kekulé structure Q+1.5a 0 
Resonance bet ween two 
Q+2.4a0 0.9a 1:0 


Kekulé structures 
Q42.6055a 1.1055a@ 1:0.4341 


Resonance among all 
five structures 


The ratio a: b is the ratio of the coefficient of structures 
A and B to that of the singly-excited structures C, D, and E. 


to the extent of 0.9a or about 1.35 v.e. over a ring 
with three double bonds. In addition, however, 
the excited structures contribute appreciably 
both to the energy (0.2055a) and to the eigen- 
function.® In a sense it may be said that all 
structures based on a plane hexagonal arrange- 
ment of the atoms—Kekulé, Dewar, Claus, etc.— 
play a part, with the Kekulé structures most 
important. It is the resonance among these 
structures which imparts to the molecule its 
peculiar aromatic properties. 


Naphthalene 

With naphthalene the calculation is very much 
more involved. The number of canonical struc- 
tures is here forty-two (Fig. 2), so that the 


Fic. 2. The forty-two canonical structures contributing to the normal state of the 


*L. Pauling and J. Sherman, J. Chem. Phys., July, 1933. 
‘The Claus centric structure, an old-quantum-theory 
analogue of which was suggested several years ago by one 


naphthalene molecule. 


of us (L. Pauling, J. Am. Chem. Soc. 48, 1132 (1926)), is 
found to make a less important contribution to the normal 
state of benzene than do the Kekulé structures. 
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secular equation is of the forty-second degree. 
We have evaluated the matrix elements in the 
secular equation by the graphical method, but do 
not reproduce them because of their great 
number. For convenience the structures with the 
maximum number of bonds between adjacent 
atoms are called unexcited structures and those 
with one, two, etc., fewer such bonds first excited, 
second excited, etc., structures. From con- 
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siderations of symmetry, we can set a number of 
coefficients equal in the complete eigenfunction, 
and thereby reduce the secular equation to the 
fifteenth degree. To obtain an approximate 
solution of this we then assume that all of the 
first excited structures have the same coefficient, 
and similarly for all of the second and all of the 
third excited structures. The resultant secular 
equation, corresponding to the function 


(Q—W)+2a 
(1/2)(Q—W)+13a/4 
(35/8)(Q— W)+163a/8 
(27/8)(Q—W)+165a/8 
(5/8)(Q—W)+31a/8 


2(0—W)+7a 


(27/8)(Q—W)+165a/8 
(21/2)(Q— W) +159a/4 
(525/8)(Q—W) +2169a/8 
(633/8)(Q—W)+19110/8 
(57/4)(Q— W) +189a/4 


To solve this we first neglect entirely the second 
and third excited structures. We then solve the 
secular equation, which is now only a cubic, and 
evaluate the ratio of the coefficients, a:b: c. 
Finally, we assume that this same ratio holds in 


(1/2)(Q—W) +13a/4 
(17/8)(Q—W) +43a/8 
(19/2)(O—W) +169a/4 
(21/2)(Q—W)+159a/4 


(35/8)(Q— W) +163a/8 

(19/2)(Q— W) + 169a/4 
(557/8)(Q— W)+2035a/8 
(525/8)(Q—W)+2169a/8 

(95/8)(Q— W)+397a/8 

(5/8)(Q—W)+31a/8 

2(Q—W)+7a 

(95/8)(Q—W)+397a/8 
(57/4)(Q—W)+189a/4 
(19/4)(Q—W) —7a/4 


the complete eigenfunction, thus obtaining an- 
other cubic secular equation, the solution of 
which is an approximate solution of the original 
secular equation. The results are summarized in 
Table IT. 


TABLE II. 


Resonance 
energy 


Total 
energy a:b:c:d:e 


Single unexcited structure +2a 0 30:60:66 
Resonance among three unexcited structures +3.3703a  1.3703a : 0.8757 :0:0:0 
Resonance among unexcited and singly excited structures Q+3.9760a  1.9760a 513 : 0.3026 :0:0 


8 
Resonance among all structures Q+4.01533a 2.01533a 8513 : 0.3026 : 0.0919 : 0.0077 


It will be noticed that the singly excited 
structures contribute somewhat over 30 percent 
to the total resonance energy instead of only 
about 20 percent as in benzene. The doubly and 
triply excited structures, however, make only a 
very small contribution. 

Pauling and Sherman have found empirically® 
that the ratio of the resonance energy of naph- 


thalene to that of benzene is 1.90, while the 
present calculation gives the ratio of 1.82. The 
agreement is really quite satisfactory. It would 
have been somewhat, but probably not a great 
deal, better if the secular equation in the 
naphthalene calculation had been solved rigor- 
ously. The discrepancy must be attributed to the 
crudity of the method of calculation, and partly 
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also to the fact that, as mentioned above, the 
resonance energy calculated by Pauling and 
Sherman is not strictly comparable with that 
found by the present method. The value of a 
obtained from this calculation is —1.6 v.e., in 
good agreement with the value of —1.5 v.e. 
obtained from the benzene calculation. 

This treatment could be applied to anthracene 
and phenanthrene, with 429 linearly independent 
structures, and to still larger condensed systems, 
though not without considerable labor. It is 
probable that the empirical rule’ of approximate 
proportionality between the resonance energy 
and the number of benzene rings in the molecule 
would be substantiated. 


THE STRUCTURE AND STABILITY OF HypRO- 
CARBON FREE RADICALS 


The dissociation of a hexaarylethane into free 
radicals, 
CR3;92R3C — 


first observed by Gomberg about thirty years 
ago, can be considered formally to take place in 
two steps. First, the bond connecting the two 


halves of the original molecule breaks; and 
second, the free radicals formed thereupon 
stabilize themselves to a certain extent by means 
of a change in structure. Most theories which 
have been proposed to explain the existence and 
stability of free radicals have more or less tacitly 
assumed that the first of these steps is of primary 
importance and that the influence of the second 
upon the dissociation is negligible in comparison. 
Consequently these theories have concerned 
themselves very largely with attempts to explain 
why the ethane linkage in the undissociated 
molecule should be weakened by the presence of 
aryl groups. Of these various attempts, two stand 
out most prominently, one based upon the 
theory of “‘affinity-demand” (Valenzbeanspru- 
chung), and the other upon the consideration 
of steric influences. According to the theory of 
affinity-demand, the aryl groups use up so much 
of the available affinity of the methyl carbon 
atoms that there is very little left over for the 
ethane linkage. This is not a very satisfactory 
explanation, however, for several reasons. It does 
not show why aryl groups should have such a 
great affinity-demand—as compared, for ex- 


ample, with alkyl groups, nitro groups, etc.—nor 
does it enable one to predict which aryl groups 
should be relatively more effective and which 
relatively less effective in promoting dissociation. 
Furthermore, it leads to definitely incorrect 
predictions in some cases. Thus, tetraphenyl- 
methane would be expected to be rather unstable, 
while actually it is an extraordinarily stable 
compound, boiling at 431° without decompo- 
sition. 

The steric theory of free radical dissociation is 
based upon the assumption that the aryl groups 
are so large that they prevent the two halves of 
the undissociated molecule from coming close 
enough together to permit the formation of a 
strong bond. There are, however, some diffi- 
culties with this explanation as well. It is roughly 
true, to be sure, that the effectiveness of a group 
in promoting free radical dissociation runs 
parallel to its size, but why, for example, should 
the linear biphenyl group have a greater steric 
effect than the phenyl group, and why should 
there be an appreciable difference between the 
a and the 6-naphthyl groups? In addition the 
stability of tetraphenylmethane is about as 
irreconcilable with this theory as with that of 
affinity-demand. 

We may conclude from this discussion that it is 
apparently impossible to base a satisfactory 
theory of the free radicals upon the assumption 
that the C—C bond in the undissociated molecule 
is weakened by the influence of the aryl groups. 
We avoid the necessity of making this assump- 
tion, however, by considering that the energy 
necessary to carry through the first step of 
breaking the bond is the same in all cases, but 
that the subsequent stabilization of the radicals, 
resulting from the second step of changing the 
structure, may vary widely. Thus we attribute 
the dissociation not to the instability of the 
undissociated molecule but to the stability of the 
radicals formed in the dissociation. 

This idea has been developed by Burton and 
Ingold,’ but not upon a very satisfactory 
theoretical basis. From their work it is not clear 
why the postulated redistribution of charge 


should occur after, but not before, dissociation, 


7H. Burton and C. K. Ingold, Proc. Leeds Phil. Soc. 1, 
421 (1929); C. K. Ingold, Ann. Reports Chem. Soc. 25, 154 
(1928). 


a. 
4 
ite 
he 
it 
’ 
+ 
he 
ar 
ri 
ie 
= 
¥ 
0) 
al 
in 
77 
Me, 
1€ 
at 
y 
a 


368 


nor is the mechanism by which this redistribution 
stabilizes the radical definitely established. In the 
following pages we describe the extension to the 
free radical problem of the quantum-mechanical 
methods which we used with benzene and 
naphthalene. It is shown that a straightforward 
calculation of the resonance energy for the 
undissociated ethane and for the free radical 
leads to values of the dissociation energy of 
various substituted ethanes in excellent quali- 
tative and semiquantitative agreement with 
experiment. This theory of the stability of free 
radicals may be considered as a refinement and 
extension of Burton and Ingold’s general views, 
for which it provides a sound theoretical basis. 

Before going into our methods of calculation 
we shall first discuss briefly the nature of the 
experimental results which we wish to explain. In 
spite of the vast amount of work which has been 
done on free radicals in the last thirty years, the 
information available regarding the degree of 
dissociation of most hexaarylethanes is meager, 
and exceedingly inaccurate. Ziegler and Ewald* 
have carried out precise photometric determi- 
nations in the cases of hexaphenylethane and 
tetraphenyldi-8-naphthylethane. In other cases 
the only data at hand are derived from determi- 
nations of the apparent molecular weights by the 
cryoscopic method. That this method is unde- 
pendable is shown by the fact that it frequently 
indicates degrees of dissociation of 130 percent or 
even higher, and in the case of the slightly 
dissociated tetraphenyldi-8-naphthylethane it 
indicates a degree of dissociation two or three 
times greater than that found by Ziegler and 
Ewald. In some cases even this source of infor- 
mation is lacking, and then recourse must be had 
to such data as the intensity of the color of 
solutions, the temperature at which the color 
first appears, the number of times the color will 
reappear after having been discharged with air, 
the rate of absorption of oxygen, and so on. 
Obviously such methods can give only the 
crudest of results. As a further complication we 
have the fact that the degree of dissociation of a 
given ethane varies considerably with tempera- 
ture, with solvent, and with concentration, and 
frequently the determinations have not been 
made under comparable conditions. In spite of 


8 K. Ziegler and L. Ewald, Ann. 473, 179 (1932). 
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all these difficulties, however, we can say with 
some assurance that the tendency toward dis- 
sociation increases in the following order: any 
ethane with less than six aryl groups < dipheny]- 


difluoryl K hexaphenylethane < tetraphenyldi- 


B-naphthylethane <_ tetraphenyldi-a-naphthyl- 
ethane = diphenyltetrabiphenylethane < hexa- 
biphenylethane. We have intentionally omitted 
tetraphenyldibiphenylethane from this list since 
we do not consider that its degree of dissociation 
has been measured with sufficient accuracy to 
warrant its inclusion. It has been studied neither 
photometrically nor cryoscopically, and so the 
estimate of 15 percent dissociation usually given 
in the literature is not to be taken very seriously. 
In most lists like the above it occupies a position 
between hexaphenylethane and tetraphenyldi-6- 
naphthylethane. 

From such crude data as are to be found in the 
literature we can calculate approximate values of 
the equilibrium constants, and hence of the free 
energies of dissociation for the various hexa- 
arylethanes. From our quantum-mechanical 
treatment, on the other hand, we obtain only the 
heats of dissociation, for which, except in the 
single case of hexaphenylethane, we have no 
experimental data. Thus, in order that we may 
compare our results with those of experiment, we 
must make the plausible assumption that the 
entropies of dissociation vary only slightly from 
ethane to ethane. Then at a given temperature 
the heats of dissociation run parallel to the free 
energies and can be used instead of the latter in 
predicting the relative degrees of dissociation of 
the different molecules. 

In carrying out the calculations we use 
essentially the same procedure as in the case of 
benzene and naphthalene. As an additional 
simplification, however, we neglect entirely all 
the excited states of the molecule, since their 
contribution to the total energy is comparatively 
small, and since they would complicate the 
calculations tremendously if retained. Another 
slight modification of the procedure is neces- 
sitated by the fact that a free radical possesses an 
odd number of electrons, one of which must 
remain unpaired. This is taken care of formally 
by introducing a “phantom orbit” X with an 
accompanying “phantom electron” which is 
paired with the odd electron.‘ In the subsequent 
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calculations it is necessary merely to remember 
that all exchange integrals which involve the 
phantom orbit vanish. The rules for setting 
up the matrix elements and the secular equations 
are otherwise unchanged. 

We shall discuss in detail only one example, 
and for the sake of simplicity we shall take the 
(actually undissociated) sym. diphenylethane. 
Let us first consider what happens when the 
ethane dissociates. In the first step, the C—C 
bond breaks and there are formed two phenyl- 
methyl radicals, which however can resonate be- 
tween only the structures A and B of Fig. 3. 


Fic. 3. The five unexcited canonical structures contributing 
to the normal state of the phenylmethy] radical. 


AH for this reaction we take to be 3.65 v.e., the 
energy required to break a C—C bond in aliphatic 
compounds. The second step of the dissociation 
consists in a change in the structure of the 
radicals, involving a change in nuclear con- 
figuration from a tetrahedral arrangement of 


This can be simplified to a cubic equation since 
from considerations of symmetry it follows that 
the coefficients of ¥4 and of yz must be the same 
in the complete function, and similarly for wc 


(5/2)(Q—W) +19a/4 
(7/4)(Q—W) +25a/4 


On solving this we find for the energy W=Q 
+2.4092a, and for the wave function y=a(¥4 
t+¥z)+b(¥o+Wp+wWe) where a and b are in the 
ratio a: b=1:1.0279. The “free radical reso- 
nance energy”’ is thus W— W’=0.5092a, and AH 
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(11/2)(Q—W) 
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bonds about the methyl carbon to a_ plane 
arrangement of all the nuclei, and in electronic 
structure through the removal of the restriction 
of the resonance to the structures A and B, all 
five of the canonical structures being now 
allowed. AH for this reaction is just twice the 
difference between the energies of a single radical 
after and before the removal of the restriction. 
The energy change for the over-all reaction is 
then the algebraic sum of the energy changes of 
the two steps. 

To evaluate A// for the second step we must 
set up and solve the secular equations corre- 
sponding to the two different states of the 
radical. In the case where we consider only the 
structures A and B the secular equation is 
easily shown to be: 


(Q—W’)+ea 1(0—W’) +11a/8 
1(Q—W’)+11la/8 (Q—-W’)+a 
On solving this we find for the energy W’=Q 
+1.9a, and for the wave function =a(~4+yz) 
where a is a normalization factor. 


When we consider resonance between all five 
structures we obtain a quintic secular equation: 


(Q—W)+a 5(Q—-W)+7a/8 4(Q-W)+a 2(Q—W)+5a/4 
(Q—-W)+ea 2(0—W)+5a/4 4(Q-W)+a s(Q—W)+7a/8 
sQ-W)+7a/8 3(Q-W)+S5a/2  (Q—-W)+a |=0. 
1(Q—W)+a (Q-W)+a 2(0—W)+5a/4 
2Q-W)+5a/4  3(Q—W)+7a/8 2(Q—-W)+5a/4 (Q—-W)+a 


and wz. We shall however, make a further 
simplification by assuming that yp has the same 
coefficient as Yc and wz, and thereby reduce the 
secular equation to a quadratic: 


(7/4)(Q—W)+25a/4) 


for the dissociation is equal to 3.65 v.e.+2 
X 0.50920 = 3.65 v.e.+1.0184a. As we shall show 
later, the value of a is probably about —1.4 v.e., 
and so the heat of dissociation of sym. diphenyl- 
ethane is approximately 2.2 v.e., a value in good 
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accord with the vanishingly small degree of 
dissociation. In view of the uncertainty con- 
nected with the values of a and of the C—C bond 
energy we shall as a rule not calculate the heats of 
dissociation of the remaining ethanes, but instead 
shall use as a criterion of the tendency of a given 
compound to dissociate simply the coefficient of 
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a in the expression for the free radical resonance 
energy. It will be noted that on account of the 
negative sign of a the dissociation increases as 
this coefficient increases. 

The results obtained from this and from 
similar calculations for a number of other free 
radicals are summarized in Table III. In the first 


TABLE ITI. 


Radical 


W-W 


Phenylmethyl 
Biphenylmethyl 
B-Naphthylmethyl 
a-Naphthylmethyl 


Fluoryl 
Diphenyl methyl 


Phenylfluory! 
Triphenylmethyl 
Biphenyldiphenylmethy] 
Tribiphenylmethyl 
8-Naphthyldiphenylmethyl 
a-Naphthyldiphenylmethy! 


0.5092a 
0.5568a 
0.6309a 
0.7515 


0.8004a 
0.8409a 


1.0164a 
1.1078a 
1.1310a 
1.1545a 
1.1773a 
1.1715a 
1.2398a 


aoa 


F444 
OOD 


OO OOOO 
FF 
ODO OOOO 


column we have listed the radical under con- 
sideration; in the second, the energy W’ calcu- 
lated for the radical with resonance restricted to 
the structures in which the unpaired electron is 
on the methyl carbon atom; in the third, the 
energy W calculated for the radical with reso- 
nance allowed among all possible states; and in 
the fourth, the free radical resonance energy 
W-W’. 

Mention may be made of a special method 
developed for setting up the secular equation for 
free radicals. The direct calculation of matrix 
elements, even by the simple graphical method, 
would be extremely laborious for such a radical 
as tribiphenylmethyl, for example, for which 496 
structures are taken into consideration. Let us 
consider three superposition patterns, each con- 
sisting of one of the substituent groups in the free 
radical plus the methyl carbon and the phantom 
orbit, the corresponding Coulomb and exchange 
coefficients being qi, g2, gz and di, de, a3, re- 
spectively. It is easily proved that the Coulomb 
coefficient g and the exchange coefficient a for the 
entire superposition pattern have the values 


= 919293 
and 


@ = 019293 + 910293 + 919203. 


For a set of superposition patterns obtained by 
combining a set for each of the three groups in 
all possible ways, the sums of the coefficients (for 
equal weights of the corresponding structures) 
are 

Yq 
and 


As an example we may use these equations to 
calculate the coefficients summed over the 64 
matrix elements corresponding to the eight 
structures of triphenylmethyl with the unpaired 
electron on the methyl carbon. Reference to the 
secular equation for phenylmethyl shows that 
and 
= Ya3;=11/4, leading to the values Yg=125 8 
and Ya= 1425/16. 

The secular equations which must be solved 
for the values of W are frequently of very high 
degree, even after all of the simplifications 
arising from considerations of symmetry have 
been introduced. In such cases we have reduced 
the equations to cubics (sometimes to quadratics) 
by equating coefficients in a more or less arbitrary 


| 
.6218a 
.7310e 
6545a 
835 1a 
.9033a 
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manner. The calculated values for the free 
radical resonance energy are, therefore, uncertain 
to the extent to which they are affected by this 
not strictly legitimate simplification. We feel, 
however, that the errors introduced in this way 
are not very large, since we have found by actual 
calculation in two cases that the energies are not 
very sensitive to small changes in the coefficients. 
For example, the energy W of the triphenyl- 
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methyl radical is given in Table III as Q 
+6.8078a. This was obtained by solving rigor- 
ously the secular equation, which in this case 
could be reduced to a cubic. If, on the other 
hand, we equate all coefficients and make the 
equation linear, we obtain for W the value 
Q+6.8064a. In the case of the a-naphthyl- 
diphenylmethyl] radical the value of 0+7.9033a 
given for W in Table IV was derived from a cubic 


TABLE IV. 


Relative values 
of coefficients 


Number of 
structures 


Position of 


Radical odd electron 


Phenylmethyl methyl carbon 
phenyl group 
methyl carbon 


Biphenylmethy] 
biphenyl group 


oF Wh 


methyl carbon* 
methyl carbon? 
naphthyl group 


6-Naphthyl methyl 


anN— 


methyl carbon* 
methyl carbon?’ 
naphthyl group 


a-Naphthylmethy! 


methyl carbon 
biphenylene group* 
biphenylene group* 


Fluory! 


methyl carbon 


Diphenylmethyl 
phenyl group 


— 


methyl carbon 
group 
iphenylene group 


Phenylfluoryl 


00 


methyl carbon 
phenyl group (para) 
phenyl group (ortho) 


Triphenylmethyl 


methyl carbon 
group 
iphenyl group 


Biphenyldiphenyl methyl] 


methyl carbon 
group 
iphenyl group 


Phenyldibiphenylmethy1 


Tribiphenylmethyl methyl carbon 
biphenyl group 


8-Naphthyldiphenylmethyl methyl carbon 
phenyl group 
naphthyl group 
a-Naphthyldiphenylmet hy] methyl carbon 
phenyl group 
naphthyl group 


5 With the naphthyl group symmetrical as in structure I, Fig. 2. 
— the naphthyl group unsymmetrical as in structures II and III, Fig. 2. 
he distinction between these two types of structure is described below in the discussion of the fluoryl radical. 
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equation, obtained by equating the coefficients of 
all structures in which the odd electron is (1) on 
the methyl carbon atom, (2) on the naphthyl 
group, and (3) ona phenyl group. If the equation 
is reduced further to a quadratic by equating the 
coefficients of all structures of the last two types 
the energy W is changed only in the fifth decimal. 
The effect of these errors which, though small, are 
sometimes quite appreciable will always be to 
make the value of W too large (that is, the coeffi- 
cient of a too small). Since the secular equations 
for W’ can be solved rigorously in all cases, this 
will cause the radicals to appear less stable relative 
to the ethanes than they should. This effect should 
be particularly pronounced in the cases of the 
a- and the 8-naphthyldiphenylmethyl radicals 
because of the high degree of the secular equa- 
tions. Furthermore, in these radicals we have had 
to give the structures in which the naphthyl 
group has the symmetrical form (as in structure 
I, Fig. 2) the same coefficient as the structures in 
which the naphthyl group has the unsymmetrical 
form (as in structures II and III, Fig. 2), and this 
simplification seems to be particularly unjusti- 
fied. At any rate it is found to be pretty far off in 
the cases of the simple a and 8-naphthylmethyl 
radicals where the calculations have been carried 
through somewhat more rigorously. This error 
has been partially corrected by making the same 
approximation in the calculations of W’; that is, 
by putting all coefficients equal and solving 
linear instead of quadratic secular equations. 

In order to make clear the manner in which we 
have simplified the treatment of the various 
radicals, we give in Table IV a statement of the 
value and significance of the coefficients occurring 
in the function y in the individual cases. It will 
not be necessary to give a similar table for the 
functions y’ since, except in the cases of the a 
and the B-naphthyldiphenylmethy] radicals which 
have been discussed above, no simplifications 
were necessary for them, and the secular equa- 
tions were solved rigorously. 


Discussion of results 


The principal source of inaccuracy in the 
treatment given the free radicals is the neglect of 
all excited structures. In the case of benzene the 
three excited structures contribute 0.2055a to 
the resonance energy, which is about 20 percent 
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of the total resonance energy of the molecule. 
The same error is introduced in W’ for the 
phenylmethyl radical. The change in W produced 
by the inclusion of the nine excited structures as 
well as the five unexcited structures is not so 
large, so that W— W’ is decreased from 0.5092a 
to 0.4455a, a decrease of 0.06a or 12 percent of 
W— W’. Similar changes in the values of W— W’ 
would result from the consideration of excited 
structures for other radicals. In consequence, so 
long as excited structures are not included in the 
treatment the values of _W—W’ cannot be 
considered to be accurate to 10 percent. This 
does not prevent us from making use of the values 
given in Table III, however, in certain cases. 
Thus for each of the radicals included in the 
table the value of W—W’ can be considered to 
give the extra resonance energy stabilizing the 
radical with an accuracy of 10-20 percent. 
Moreover, in some cases differences in the 
W—W’ values for different radicals can be 
trusted to the same extent; namely, when one 
radical differs from another only by the addition 
of a group, so that all of the structures considered 
for the first radical are retained for the second, 
as well as other structures. Thus the difference 
of the W—W’ values for biphenylmethyl and 
phenylmethyl should give with considerable 
accuracy the extra stability resulting from the 
substitution of biphenyl for phenyl, and a similar 
comparison should be reliable for a-naphthyl and 
B-naphthyl, a or B-naphthyl and phenyl, fluory! 
and diphenyl, etc., but not for a-naphthyl and 
biphenyl, inasmuch as in these completely differ- 
ent groups the effect of neglect of excited 
structures may well mask the small differences 
under consideration. When this restriction in the 
interpretation of the W—W’ values is borne in 
mind, it is found that the agreement between the 
rather rough theoretical calculations and the 
existent empirical information is excellent, and 
provides sound evidence that resonance of the 
type considered is the primary factor in the 
stabilization of the hydrocarbon free radicals. A 
number of details in the comparison are treated 
in the following paragraphs. 

(a) The calculation for triphenylmethyl leads 
to a value of 3.65 v.e.+2.2156a@ for the heat of 
dissociation of hexaphenylethane. Equating this 
to the experimentally determined value of 0.5 
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v.e., we obtain as a result that a= —1.4 v.e. 
This is in fair agreement with the value of —1.5 
v.e. obtained from the calculation for benzene; 
the discrepancy may be due in part to the fact 
mentioned before that the evaluation of a from 
the benzene calculation is not strictly legitimate, 
and in part to the neglect of excited structures, 
which (from analogy with phenylmethyl) are 
expected to decrease the coefficient of @ in 
W—W’, and so increase the absolute value of a. 

(b) Except for comparisons of naphthyl and 
biphenyl, the amounts of dissociation indicated 
by the values of W—W’ for various radicals are 
in excellent qualitative agreement with experi- 
ment. It is undoubtedly significant that, as is 
demanded by experiment, the variations in the 
free radical resonance energies obtained by 
increasing the number of aryl groups is large 
compared with the variations obtained by 
changing the nature of the groups already present. 

(c) In order to make a quantitative com- 
parison of our results with experiment, let us 
calculate the ratio of the dissociation constants 
corresponding to a difference in the free radical 
resonance energies of 0.01a. If we let the sub- 
script 1 refer to the more highly dissociated 
ethane and 2 to the less highly dissociated, we 
have the following relation: 


RT In (K2/K;) 


Since a is equal to — 1.4 v.e. or — 32,000 calories 
per mole, this equation becomes 


RT = 640 cal./mole 
which at 20°C leads to 


In (K,/Ke)=1.1 or K,/K2=3.0. 


From Table III we see that the difference 
between the free radical resonance energies of 
tribiphenylmethy! and triphenylmethy] is 0.07a. 
Hence Ziegler and Ewald® 
found that at 20°C the value of the dissociation 
constant for hexaphenylethane in benzene solu- 
tion is 4.1 10-4 and consequently we calculate 
for hexabiphenylethane a value of K =2.2X 10 
X4.1X10-*=0.90. This value is probably too 
low as the compound is reported to be completely 
dissociated ; the error may not be large, however, 
since a dissociation constant of 0.90 would lead 
to 91 percent dissociation in 0.05M solution. 
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The solutions actually worked with have usually 
been more dilute than this, and so our calculated 
value may be within the experimental error. 

(d) The ability of the present method to 
distinguish between apparently small differences 
in structure is brought out forcibly by a 
comparison of the relative effectiveness of the 
a- and the 8-naphthyl groups on the one hand, 
and of the fluoryl and the diphenylmethy] groups 
on the other. The greater dissociating power of 
a-naphthyl is qualitatively correlated with the 
fact that for it there are seven unexcited struc- 
tures with the odd electron on a naphthyl carbon, 
whereas there are only six for B-naphthyl. The 
case of fluoryl is particularly striking. For the 
fluoryl radical it is possible to write six more 
unexcited structures than for the diphenylmethy] 
radical. These extra structures are all of the type 


x 


That is, they all have a double bond in the ortho 
linkage between the phenyl groups, and they all 
have the odd electron on one or the other of the 
phenyl groups. Hence it would seem that the 
radical should be stabilized by these extra 
possibilities of resonance, while the ethane should 
be left practically unaltered, and that conse- 
quently the dissociation should be increased. 
Actually, however, the opposite is the case. On 
putting in the phantom orbit and drawing all the 
bonds it is found that no matter where the 
phantom orbit is put with respect to the other 
orbits all the possible structures cannot be made 
canonical simultaneously. There will always be 
six non-canonical structures which must be 
represented as linear combinations of the 
canonical ones, and so the structures suggesting 
extra resonance are eliminated. Just which six 
structures are non-canonical will depend upon 
the position given the phantom orbit. If it be put 
in the symmetrical position indicated above, the 
non-canonical structures will be those in which 
there is a double bond between the two phenyl 
groups. Two of these structures turn out to be 
combinations of excited structures and so are 
neglected. The remaining four are combinations 
of one unexcited and one excited structure each. 
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These excited structures are also neglected, and 
the four unexcited structures are distinguished 
from the remaining eight by being given a 
different coefficient. This is the reason for the 
fact noted in Table IV that for the fluoryl 
radical there are three instead of two coefficients. 
Since, however, the coefficients for the two 
different types of structures proved to be nearly 
the same here, they were put equal to each other 
in the calculation for the phenylfluoryl radical. 
This explains why the fluoryl group should not be 
more effective in promoting dissociation than the 
diphenylmethyl group, but to see why it is 
actually less effective we must turn our attention 
to the rules by which the matrix elements are 
calculated. Since we have shown that we can 
eliminate all the extra resonance structures of the 
fluory] radical, the structures which we have left 
are exactly the same as those in the case of 
diphenylmethyl, except for the ortho bond 
between the two phenyl groups. Therefore the 
secular equation is the same except for the effect 
of the interchange integral between the atoms at 
either end of the bond, and this integral will 
always come into the matrix elements with a 
negative sign. Consequently the coefficient of a 
in any given matrix element for the fluoryl 
radical will be less than in the corresponding 
element for the diphenylmethyl radical. As a 
result of this, the coefficient of a will be decreased 
in both the energies W and W’, and since this 
decrease is greater in the former case than in the 
latter, the coefficient in the extra resonance 
energy W—W’ will also be decreased. 

(e) Although the theory described in this 
paper is able to explain qualitatively and semi- 
quantitatively the dissociation of a considerable 
number of free radicals, it proves to be definitely 
inadequate in some cases. For example, di-tert.- 
butyl-tetraphenylethane has been found to be 
appreciably dissociated at room temperature, 
although our calculation would give exactly the 
same value of the free radical resonance energy 
as in the case of the undissociated tetraphenyl- 
ethane, since the presence of the tert. butyl 
group does not affect the number of resonating 


structures nor the values of the matrix elements. 
The same considerations also apply to decapheny!I- 
butane, which has been found to be completely 
dissociated into pentaphenylethyl radicals. We 
consider that in these cases the steric effect may 
enter and be of dominating importance. There is 
no doubt that the tert. butyl and the tripheny]- 
methyl groups are relatively large, and their 
presence in these compounds may tend to keep 
the radicals apart. 

Furthermore, there are a number of interesting 
free radicals other than the hydrocarbons to 
which our treatment could be extended. In order 
to carry through similar calculations for com- 
pounds containing oxygen or nitrogen, however, 
we would have to introduce a number of ad- 
ditional parameters, corresponding to exchange 
integrals other than that to which we have given 
the value of a, which would make the calcu- 
lations more difficult and the results less valuable. 

Mention may be made of the interesting 
observation of Wieland that triphenylmethyldi- 
p-anisylamine does not dissociate appreciably 
except at temperatures above 140° 


N(C,H 
=(Ce6H3;)3C + N(C6H,OCHs)s 


although the free radicals formed can exist to an 
appreciable extent in the monomeric form even 
at room temperature—a fact which has excited 
considerable comment. The explanation does not 
involve a difference in resonance in this com- 
pound and the corresponding ethane and hydra- 
zine, but instead a difference in the energy of the 
carbon-nitrogen bond and the mean for the 
carbon-carbon and_ nitrogen-nitrogen bonds, 
resulting from additional ionic character of the 
unsymmetrical bond. This amounts to about 
0.5 v.e.,9 and so is large enough to prevent 
dissociation. 

On the other hand, the stability of positive and 
negative triarylmethyl ions is due to the same 
resonance effect as that of the neutral radicals 
discussed in this paper. 


® See the preceding and the following paper of this series. 
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Predissociation and the Crossing of Molecular Potential Energy Curves 


O. K. Rice, Chemical Laboratory, Harvard University 
(Received March 16, 1933) 


The calculation of the width and shape of a line which is 
broadened because of predissociation has been extended so 
as to include perturbations of intermediate size. By 
“perturbations of intermediate size” we mean perturbations 
which are not so great that lines (of given rotational 
quantum number) belonging to two adjacent vibrational 
levels are broadened so as to overlap appreciably; they may 
be any size up to that limit. These calculations have been 
applied with certain simplifying assumptions to the case 
where a potential curve giving molecule formation is 
intersected by a repulsive curve, the curves which cross 
defining the unperturbed energy levels and wave functions. 
It is found that the discrete lines are not only broadened, 
but they are shifted in position, slightly distorted in shape, 
and there is present a series of subsidiary maxima of the 


absorption coefficient. The exact amount of broadening of a 
line depends very greatly upon its energy relative to the 
energy at which the potential curves cross, and in general a 
line which is much broadened will also be much shifted in 
position. In molecules we may expect tc have isolated 
groups of rotational levels in which the central level is 
sharp, and the rotational levels on either side become more 
and more diffuse, finally fading out. The predissociation 
phenomena in iodine chloride are discussed on the basis of 
the above theory. In particular, two groups of sharp and 
diffuse levels, such as described, are considered and are 
respectively ascribed to the isotopic molecules, ICI® and 
ICI*", There is some difficulty in describing the phenomena 
in iodine chloride quantitatively, but qualitatively the 
experimental results appear to accord with the theory. 


HE question of what happens when two 

potential energy curves of a diatomic 
molecule cross has interested spectroscopists for 
some time. The Franck-Condon principle, so 
useful in explaining intensities in band spectra, 
says, when applied to the case of radiationless 
transitions such as occur for example in predis- 
sociation, that transitions are likely to occur from 
the electronic state corresponding to one poten- 
tial energy curve to the electronic state corre- 
sponding to the other if the total energy of the 
system is about equal to the energy at which the 
two curves cross. This principle has been used by 
Herzberg! and others to explain the variation in 
the widths of the rotation lines in various cases 
of predissociation. Nevertheless, little quanti- 
tative work has been done on the subject. It is 
the purpose of this article to investigate quan- 
titatively just what does happen in the neighbor- 
hood of the crossing point of two potential energy 
curves in a particularly simple case of predis- 


‘See Herzberg, Zeits. f. Physik 61, 604 (1930); Ergeb. d. 
exacten Naturwiss. 10, 237 (1931). Other references are 
given in the latter article. The question has been treated 
from another point of view by London, Zeits. f. Physik 74, 
158 (1932). In this connection see also, Neumann and 
Wigner, Phys. Zeits. 30, 467 (1929), and Zener, Proc. Roy. 
Soc. (London) A137, 696 (1932). 


sociation. We shall find that a given rotation 
level is not only broadened, but is somewhat 
distorted in shape, is shifted in position, and has 
a rather complicated secondary structure. An 
attempt will then be made to apply the results 
to the case of predissociation in iodine chloride, 
observed by Brown and Gibson.’ 


§1. GENERAL THEORY 


The phenomenon of predissociation, first dis- 
covered by Henri and Teves,’ was explained as 
a sort of molecular Auger process by Bonhoeffer 
and Farkas,‘ and considered from a quantum 
mechanical point of view by Kronig.' I discussed 
it further,® and gave a fairly rigorous derivation 
of the formuia for the amount of broadening of 
the rotation lines as a function of the perturba- 
tion matrix component of the disturbance causing 


the effect. In order to handle the present situ- 


2 Brown and Gibson, Phys. Rev. 40, 529 (1932). 

3 Henri and Teves, Nature 114, 894 (1924). 

4 Bonhoeffer and Farkas, Zeits. f. physik. Chemie 134, 
337 (1928). 

5 Kronig, Zeits. f. Physik 50, 360 (1928). See also 
Wentzel, Zeits. f. Physik 43, 524 (1927); Fues, Zeits. f. 
Physik 43, 726 (1927); Wentzel, Réunion Internationale de 
Chimie Physique, Paris, 1928, p. 121. 

6 Rice, Phys. Rev. 33, 748 (1929); 35, 1551 (1930). 
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ation, some of the approximations and special- 
izations made in this theory need to be removed, 
and to this task we now address ourselves. 

We consider a diatomic molecule, or a pair of 
atoms, and make use of an approximate wave 
function, Y, which may be written in the form 


y= OR. (1) 


Here © is the wave function for the electrons 
and the rotation, obtained by holding the dis- 
tance, r, between the two nuclei fixed. It depends 
therefore on the rotational and electronic coor- 
dinates, and on 7 only as a parameter. R is the 
radial part of the wave function and is a function 
only of 7, though it depends on the quantum 
numbers of 0. 

A © belonging to a given electronic and 
rotational state of the molecule, distinguished by 
a subscript, say k, will obey an equation of the 


form 
(Ho— Ux) O.=0, (2) 


where Hp is the Hamiltonian operator for the 
electronic and rotational coordinates, and U;, 
(a function of r) is the energy. 

A given R which goes with 0;, say Rim, may 
be written in the form Fi»/r, where Frm obeys 
the equation 


— U;.) Fim =9, (3) 


where 7?=87°M/h?, and M is the reduced mass 
of the system. 

We now write, in formal fashion, the wave 
equation which is obeyed by the approximate 
wave function Pim = as 


(IT — Exm)Wim =0. (4) 


But this cannot be an exact equation. At the 
very best, in. separating the electronic and 
rotational coordinates from the coordinate r we 
have made an approximation. We shall let an 
exact wave function y’ obey the equation 


(1+V—E’)y'=0, (5) 


where V is a perturbation term and E’ the per- 
turbed energy. 

We now come more specifically to the con- 
sideration of a particular case in which predis- 
sociation may occur. We will suppose that we 
have two potential energy curves U; and U2 
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with corresponding wave functions ©; and @,. 
We further suppose that U; is the repulsive type 
curve, or else that we are interested in an energy 
region above its energy of dissociation, while U, 
is the type giving molecule formation and we are 
below its dissociation energy. Any R, let us say 
Rea, or simply Ra, which goes with U2 will be 
an eigenfunction belonging to a discrete state, 
while those which go with Uj, designated in 
general by R, will form a continuous set. To 
normalize these we shall assume that 7 has a 
largest possible value which we shall call 7,,; this 
converts the continuous set of energy levels to a 
very close-spaced discrete set, the energy, «,, 
between two adjacent “‘continuous” levels being 
given by 

(6) 


where E, is the energy of the continuous state 
and U;(2) is the asymptotic value of U;. We 
normalize R, by setting 


Rerdr= Ffdr=1. (7) 
0 


Now the interaction, represented by the per- 
turbation V, causes, as is well known, an effective 
“‘broadening”’ of the discrete states. That is, in 
the exact wave function, ¥’, which can be set up 
as a linear function of discrete and continuous 
states, the coefficient of a given discrete state 
has an appreciable value over a certain finite 
energy range in the neighborhood of the original 
energy, E.. By solving the perturbation problem 
one is able to find the ‘‘width” of these broadened 
discrete states. We define certain matrix com- 
ponents, as follows, 


- for VO,dr, V22 V Ocdr, (8) 


Ve = f V = f @.* V O,d7, (9) * 


Veo = FwyFdr, vaa= Fim Fadr, (10) 
“0 0 


Usa = | adr, (11) 


* We assume 2, real. If it is not the necessary modifica- 
tion may be readily made. 
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the integrals (8) and (9) being taken over all 
possible volume elements, dr, in the coordinate 
space defining the ©’s, and the asterisk denoting 
the conjugate complex. If we assume 71; and v2 
absorbed into U; and Us, respectively, (we shall 
write V; to mean U,+2;; and similarly with V2) 
then it may be readily shown’ that the matrix 
components v,, and vga do not enter into the 
problem and the only ones we have to consider 
are those of the type d... 

We consider the case where the discrete states 
are so far apart in comparison to the broadening 
suffered by them that we need to take into 
account only one of them at a time.® This is, of 
course, a restriction on the size of the per- 
turbations, but the only one necessary to make; 
such perturbations may be called “‘perturbations 
of intermediate size.’’ We also make the assump- 
tion that there is only one set of continuous 
states involved. We can then write for a partic- 
ular say Wn, 

Vn SanWa + (12) 


Ya and y, being written for the complete wave 
functions of form (1) corresponding respectively 
to the discrete state and to the continuous states. 
If E, is the energy that goes with y,, then one 
may show, by application of the perturbation 
theory, that 


San(Ea— En) =0 (13) 


En En 


it being understood that the singularity in the 
integrand is handled in such a way that no 
contribution to the integral results from it. This 
gives us the values of the perturbed energy levels. 
For, with any given value of E, it allows us to 
find the value of 8,, and thus for any particular 


"Zener, Phys. Rev. 37, 560 (1931). Morse and Stueckel- 
berg, Ann. d. Physik 9, 581 (1931). This method of 
handling 2, avoids a troublesome difficulty previously 
encountered in the theory. 

* We consider only the case in which levels of different 
rotational quantum number do not interact with each 
other and so do not need to be considered at the same time. 
The energy between discrete states therefore refers to the 
energy between two adjacent vibrational levels. See 
reference 6 and §3. 


and, for all s, 


= E,) +0 =0. (14) 


These form an infinite set of simultaneous equa- 
tions, which are to be solved for the S’s, furnish- 
ing thereby a series of values of E,, consistent 
with non-zero solutions for the S’s. This problem 
I have solved on a previous occasion® on the 
assumption that v, did not depend on the value 
of the subscript s. As this assumption is too 
restricted for our present purposes, I will here 
present just a sufficient outline of the solution 
to show how it must be modified. From Eggs. 
(13) and (14) one may write 


(Ea E,) Vea? / (E, E,) =0. 


Now contributions to the summation in this 
equation will occur in two ways. There will be 
contributions from regions in which E, is very 
different from E, and a contribution from the 
region where E,—E,, is of the order of e,. If we 
let a, be the least positive value of E,—E,, and 
set B,=a,—e,/2, then it may be shown (see 
reference 6, especially p. 755 of the earlier 
article) that this latter contribution is —v,,47(2/€n) 
tan (78,./en), where v,,4 is the value of v,2 and e,, 
the value of «, when E,—E,, is very small. If we 
allow r,, to become very large and hence e, to 
become very small we may write the sum as an 
integral, getting 


(15) 


TBn ved? d(E,—Vi(%)) 
Ea—Entvna — tan =(), 
0 


16 
E,-E, és 


region we see how the perturbed levels are placed 
with respect to the unperturbed. It should be 
noted that £,/e, approaches a definite limit for 
any given value of E,, as e, goes to zero. 

We are really interested, however, in the 
absorption to the state with energy E£,, from the 
ground state and, indeed, for the special case in 
which, of the original eigenfunctions, only the 
one for the discrete state combines with the 
ground state. If this is the case, then we may 
take the absorption to the state E, as propor- 
tional S,,” and hence if we can find this quantity 
as a function of E, we shall have accomplished 
our purpose. Now from Eq. (12), since we 
must have /y,’dr=1 and remembering the 
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orthogonal properties of ~~ and the y,, we get 
Savr+ > »Ssn?=1. Substituting for S,, from (14) 


1 = San?+ San* Ved? / (E, Ey}. 7) 


Now starting with the state for which E,—E, 
has the smallest positive value E,—E,=€,/2+8,, 
the next value of E,—E, is 3e,/2+8,, etc., and 
similarly for E,<£, the values run ¢,/2—8,, 
3e,/2—B;, etc. It is seen that the series can be 
made to converge in as small a range of E,—E,, 
as we desire by making r,, large enough and 


En Und En ad Ved" )) 


This holds, without further restriction, for the 
effect of a single continuum on a discrete state, 
the perturbations being of ‘intermediate size.” 


§2. EVALUATION OF THE MATRIX COMPONENTS 


The purpose of this paper is to discuss the 
case of special interest which arises when the two 
curves V; and V2 cross each other as shown in 
Fig. 1. For the unperturbed state of the system 
we use the system which would be obtained if 
the two atoms forming the molecule exerted no 
force on each other,'® then we bring in the mutual 
potential energy V as the perturbation energy. 


® Dirac, Zeits. f. Physik 44, 594 (1927) arrived at the 
same formula for the probability of scattering of a particle 
by an atom which can absorb and reemit it. It appears, in 
fact, that our result can be derived from Dirac’s if we 
assume in his calculation that the particle which is reemit- 
ted is not the same one which is absorbed. Thus we could 
have absorption of a light particle, then predissociation 
(essentially equivalent to emission of another particle) the 
matrix component for the former process being different 
from the latter. If these two matrix components are 
different then it is no longer necessary for the matrix 
component for the predissociation (or reemission of the 
second particle) to be small (it can be of ‘intermediate 
size’’) which removes an apparent limitation on Dirac’s 
derivation. I have given the other deduction here, however, 
in the first place because it seems interesting to look at the 
problem from the point of view of stationary states, in the 
second place because it seems worth while to have a proof 
which avoids the use of the Dirac 6-function. See also 
Weisskopf and Wigner, Zeits. f. Physik 63, 54 (1930), 
especially the footnote, p. 64. 

10 This does not prevent us from bringing in the exchange 
phenomena in the proper way. We can also refer the 
electron coordinates to axes moving with the nuclei thus 
getting ordinary molecular states for our unperturbed 


O. K. 


RICE 


hence e, small enough. So we can, in the limit, 
take v,q outside the summation, as 7,7. The sum 
can be evaluated as in our former papers, and 
the first term becomes negligible in the limit as 
r,, goes to infinity (e,>0). We thus find® 


(1 +tan? — (18) 
En En 


Substituting now the value of tan (78,/e,) from 
(16) we get, finally,® 


(19) 


With this scheme V; and V2 consist entirely of 
the terms 71; and 222. 

We shall assume in making the calculation 
that for the important region of r near the point 


r 


Fic. 1. Potential energy curves illustrating the situation 
in §2. Schematic drawings of the eigenfunction, F,, for a 
discrete state of energy Ez, and the eigenfunction F, of a 
continuous state with energy E, are shown, drawn on the 
energy levels as zeros of ordinates. The value of v,« depends 
on the integral of the product of F, and Fa. It will be seen 
that as E, increases and r, thus decreases the phase rela- 
tionship of F, and Fy will change periodically and the 
value of v.~ will depend in a quasi-periodic way on ’,—"«- 
See Eq. (23) and Fig. 2a (noting Eqs. (26) and (27)). 
The relation between v,,~ and rn—7rq is, of course, the same 
as that between v,~ and r,—ra, the subscript ” here really 
being a special case of the subscript s. The situation de- 
picted here does not correspond very well with the sim- 
plifying assumptions made in §2, but the position of the 
wave functions will roughly correspond to that producing 
the maximum of A cos (A{rn—ra}) which occurs at \(?»—7¢) 
=—7 in Fig. 2a. 


states. In so doing, we bring in other perturbations due to 
the Coriolis forces, which may or may not be of importance. 
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of intersection, which we shall designate as 70, 
the value of v12 may be taken as constant" and 
the curves V; and V2 may be considered as inter- 
secting straight lines. Under these circumstances 
it is possible to evaluate v,~ and the integral 
that appears in (16) and (19). 

If the slope of the curve V;, at ro is V;’, then 
if the approximations indicated in the above 
paragraph hold good we may write for an eigen- 
function” of the Eq. (3) 


Fin= Nim cos | |}(r—rn)z}dz, (20) 
0 


where 
1m =TotEnm/ Vi! (21) 


and is the point at which V;=Eim, as we shall 


V.a= dzc 
Vi'| +| V2’|) 


The integral °dzcos (z*—xz) has been 
evaluated by Airy™ over a considerable range 
of the parameter involved and the range may be 
extended by the Wentzel-Kramers-Brillouin 
approximation method." A discussion with a 
table of values is given in Appendix III. 

We now turn to the evaluation of 


Os 


let the zero point of energy be the energy of the 
crossing of the potential energy curves. The 
plus or minus sign is to be used according as V’;’ 
is positive or negative. Ni» is the normalizing 
factor. 

If v32 is constant, and if V;’ is negative and 
V2’ is positive as in Fig. 1, we have 


F,Fadr, (22) 


the integration with respect to r being extended 
from © to —®, which is permissible as there is 
but little contribution to the integral, excepting 
from the region between the points r, and ra. 
As shown in Appendix I this may be reduced to 


3+ z (23) 
(1/| +1/| 


. 


ved? d(E,— )) 
£,—E, Es 
To a good approximation we can replace this by 
1 
f —— d(E,-E,). 
E,—E, 
As shown in Appendix II we thus get 


B. dz cos dt sin {AO (24) 
0 0 


E,—E,, és 
where By N Vi" | + | (25) 
and A= (3n?)8(1/| Vi’ | +1/| Vo"|)-, (26) 


r, being equal to ro+-E,,/ V;’. The evaluation of the second integral in (24) is discussed in Appendix 


III. If we also write 


cos (23-+2)ds = Acos(x) (27) 
0 


" This is similar to the assumption commonly made in 
calculating the absorption intensities in band spectra. 

See Watson, Bessel Functions, Cambridge, 1922, p. 
188. This solution of Eq. (3) has the general properties of an 
eigenfunction, inasmuch as it becomes small when r 
becomes small, and while it does not actually become zero 
for r=0, it isa sufficiently good approximation. F;»,, has the 
form (20), of course, only in that region where the potential 
curve is approximately a straight line. 

For the continuous curve the energy must be determined 
by the boundary condition at ho, and the normalizing 


factor is determined by the way this form joins on to the 
asymptotic form of Fy, (that is, Fy, in this case). The 
amplitude of the wave for large values of r is fixed by Eq. 
(7), and is the same for all energies. If the wave-length is 
short enough the value of N, which gives the amplitude 
near 7 will vary from one value of E, to another as (Ey 
—Vi(~))*. See, e.g., Nordheim, Zeits. f. Physik 46, 842 
(1927). 

13 Airy, Trans. Camb. Phil. Soc. 8, 598 (1848). 

14 See Kramers, Zeits. f. Physik 39, 828 (1926). 
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and sin (28) 
0 


then," from (19), (23) and (24), 


(Ea- 
) | +( 


€ 


This expression gives the absorption coefficient 
of the single broadened, shifted, and distorted 
predissociation line whose unperturbed energy 
is Ea(= V2'(ra—ro) by (21)) as a function of the 
total energy E,(=Vi'(r2—70) by (21)) of the 
energy level which is the final state in the absorp- 
tion process. 

The best way to appreciate the behavior of 
San? is to study a few concrete examples. S,,, it 


En) A sin(A{ — 


Acos(A{%n— Ta} ) 


is to be remembered, is a function of E,, since 
Eq is a perfectly definite energy for any discrete 
line; Sa,? plotted as a function of E,, will give 
the shape of the broadened discrete line. Let us 
now consider what happens when Ey has such a 
value that when E,=£,, we have Acos(A —1a}) 
at its first* zero point. To distinguish this special 
case we shall designate this value of Ez as E,,. 
In Fig. 2(b) we have plotted Sy,,” as a function 


2 \-2 


Asin [a 
10r 


Absorption Coefficient 
(Arb itrary Units) 


i 


to oo 


(b) 


of ME,—E.,)/Vi', the increments of which, it 
will be observed, are equal to the increments of 
We take as the value of B,. 
The ordinates are in arbitrary units. Just 
above the curve for Sza,,.2 we have plotted 
and against 


15 In the evaluation of this expression we note that B, is 
approximately independent of E,. See Eq. (6) and reference 
12. 


MEn— Ea,), | Vi'| 
Fic. 2. 


\(rn—Ya,) in such a position as to bring out the 
relationship between the various curves. Now 
suppose we consider another case in which Fu 
has a value somewhat displaced from its former 
value. Then the point at which E,=£, is 
removed by a certain amount from the value of 
E, for which As.({r,—ra)=0, and this is 


* The general features of the phenomena will be the same 
for any zero point. 
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different than the amount by which £, is dif- 
ferent from E,,. The latter quantity is, of 
course, equal to Ey—Ez,, while the former is 
equal to (Ea—Ea,)(| V1’ | + | V2’ | )/| | . For the 
value of E, where Aws(A{fn—ra}) is zero is 
E,,—(Ea—£a,) | V2"|, this being the point, 
as is readily verified, where 7,—ra equals the 
original zero-point value of 7,—r4,. Several in- 
tensity curves, as they appear in the immedi- 
ate neighborhood of the energy for which 
Acs(\{tn—Yra})=0 are exhibited in Fig. 3. 
These curves are based on the assumption that 
=| Ve’|, and again B,, is taken as —5Vj’/2. 
The maximum in each case corresponds to the 
infinitely intense and narrow line of Fig. 2, now 
displaced and broadened. Thus we see that the 
apparent position of the discrete line will, unless 
the line happens to be in the special position, 
E,,, be displaced toward the zero-point for 
by an amount AE,, where 


AEa=(Ea—Ea)(| Vi'| V2'|)/| V2'|, (30) 


the equality sign occurring for B, =. If Eq is 
not too greatly different from Ez., the rest of the 
intensity curve will not differ greatly from the 
curve in Fig. 1 except that it will be displaced 
by the same amount that the zero point of 
Aws(\\{tn—La}) is displaced, i.e., by an amount 
(Ea—Ea,) | Vy | / | | . As the value of Ea—Ea, 
grows greater the adjacent intensity maximum 
develops while the maximum near Eq, grows 
smaller. When E,—Ea, has such a value that 
for E,=Eaq is near its maxi- 


and the half absorption width AE, of the line 
will be given by 


AE, (32) 


as may be readily seen from (29) by subtracting 
the one value of E, where S,,2 has half its 
maximum value from the other. Now we have 
already noted that if B, increases (for example 
by increasing v,2) E,(max) tends to go to such 
values that Avos(A{7,(max)—ra}) approaches 
zero. It is readily seen why this happens from 
(31). As this happens E,—E,(max) tends to 
approach a constant value, in terms of which, 
B, and A,in(A{r,(max)—7ra}) it is possible 


(Ea— E,(max))/Bn—Acos(A {%n(max) — 7a} )A sin(Af7,(max) —ra}) =0 
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mum then both maxima in the intensity curve 
will be approximately of equal size. 

It is interesting to consider a little further the 
behavior of the expression (29) when the line is 
so narrow (as is actually practically the case in 
Fig. 3) that and A 


600} 


(a) (b) (c) 


Absorption Coefficient 
W 


-15 -.09 -.24 -.20 338 
NE Eq.) /\; | 
Fic. 3. The broadened predissociation line in various 
cases. (a) MEa— Ea) =0.15; (b) MEa—Ea)/|Vi'| 
=0.3; (c) A(Ea—Ea,)/| Vi'| =0.6. Units for the absorption 
coefficient arbitrary but same as in Fig. 2b. 


are nearly constant over the width of the line. 
Then the energy £,(max) at which the absorp- 
tion coefficient has a maximum is given very 
closely by 


(31) 


to evaluate from (31). 
Asin(A{?,(max)—ra}) also tends to approach 
a constant value as B, increases. From (31) and 
(32) we have 

2(Ea—E, (max) )? 


—ra}) | 


One thus arrives at the curious result that the 
larger B, and hence the perturbation matrix 
component v2. the sharper is the displaced 
discrete line, though at the same time its 
intensity decreases, as is readily seen from (29). 
It is a little difficult to see at just what point 
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this breaks down. Since the width of the lines 
becomes smaller as B, increases, it would at 
first sight seem that our criterion that two 
discrete lines should not be so broadened as to 
overlap appreciably would be automatically 
fulfilled. But the secondary structure of Fig. 2 
must be included in the breadth of any discrete 
line. As B, increases the secondary structure 
increases in prominence, the secondary maxima 
becoming stronger and narrower, and as soon as 
the secondary structure of an adjoining discrete 
line becomes as important in the neighborhood 
of that discrete line itself the calculation must 
of necessity break down. It is by no means 
certain, however, that this criterion for the 
validity of the treatment used is stringent 
enough. 

It should be stated that the spacing between 
the secondary maxima will quite often or perhaps 
usually be of the same order of magnitude as 
the spacing between adjacent vibrational levels, 
which means that the secondary maxima will 
probably not be of importance in practical spec- 
troscopy, though we conjecture they may have 
something to do with the phenomena in iodine 
chloride, as mentioned in §4. 

It is interesting to note that the spacing of 
the maxima of the secondary structure is closely 
related to what the spacing of energy levels of 
the upper adiabatic, non-crossing curve would 
be if it were the unperturbed curve, but the 
nature of the wave functions is different, and the 
phenomena involved could be properly under- 
stood only if the effect of all the discrete states 
could be studied. 


$3. VIBRATION-ROTATION LEVELS OF A 
MOLECULE 


The actual appearance of the term diagram of 
a diatomic molecule, to be expected on the basis 
of the theory presented, can be readily deduced 
from the above discussion. The theory is based 
on the assumption that a discrete energy level 
interacts with only one continuum; consequently 
it can be applied only to molecules in which 
there is a selection rule operating,”® such that a 


16 In general such a selection rule will operate, at least if 
Hund’s Case a holds. For the value of J and its projection 
on a given axis must be the same for the two states which 
interact, the only ambiguity remaining, with given 
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discrete state with given rotational quantum 
numbers interacts only with continuous states 
with the same rotational quantum numbers. 

If the rotational quantum number (for total 
angular momentum) is J, then in order to get 
the proper effective potential energy curve we 
must add to both V; and V2 the quantity 
J(J+1)/n’r’. Thus the point, 7o, of intersection 
of the potential energy curves is independent of 
J, but the energy of intersection is increased 
over that for J/=0 by an amount J(J+1)/n’7,°. 
For our present purposes we may ignore the 
change of slope of V; and V2 in the neighborhood 
of 7) and assume that the curves are corrected 
simply by adding the amount J(J+1)/7?r,?. 

Suppose a given rotational line is at such a 
position that =O when E,=E,. 
Then the line with the next higher value of J 
will have an energy such that when E,,=E, the 
quantity A.os(A{7%,—1a}) is not equal to zero; for, 
in general, 7,.~70, where 7, is the value of 7 at 
the minimum of the potential curve, and the 
energy of the line increases at a different rate 
than the energy for the crossing point 7». Thus 
in the molecule we will have one rotation level 
sharp, and it will be surrounded on either side 
by levels growing more and more diffuse. The 
diffuse levels will be displaced in position, and it 
may be readily verified that if r.<7ro the spacing 
between levels will be less than normal; in fact, 
one could even have reversal of the order, the 
levels of lower rotational quantum number 
having higher energy. 


§4. APPLICATION TO IODINE CHLORIDE 


We shall now consider, in the light of the 
theory just developed, the rather complex set of 
predissociation phenomena which has_ been 
found by Brown and Gibson? by examination of 
the absorption spectrum of iodine chloride. 

In describing these phenomena we shall follow. 
as far as possible, the standard notation used by 
Brown and Gibson. v” and v’ will be the re- 


potential energy curves, after J and its projection have been 
determined, being the orientation of the spin in one 
direction or the other referred to J, which gives rise to 
A-type doubling. But there will be a selection rule for the 
spins which prevent the two states due to A-type doubling 
from both interacting, at least appreciably, with any one 
state belonging to the other potential energy curve. 
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spective vibrational quantum numbers and J” 
and J’ the respective rotational quantum 
numbers of the initial and final states giving rise 
to an absorption line. B,, will be used for the 
quantity h?/8x°J,, where J,, is the moment of 
inertia for the v’ vibrational level. Brown and 
Gibson have given their results in the form of 
wave numbers of observed lines. As we shall be 
chiefly interested in energy levels, all values 
here given, unless otherwise noted, are values of 
energy levels, obtained from the observed 
spectral lines by taking into account the energy 
of the initial states; the values are in wave 
numbers (cm™') and are referred to the state, 
v’=1, J’ =0, as the zero of energy. » will 
denote the energy of the band head, quantum 
number v’, 7 will denote the distance between 
atoms, and r, the equilibrium distance for the 
excited (’) energy level. 

Brown and Gibson have found a number of 
bands corresponding to an absorption *IIp-—'>, 
which they have designated as System II. They 
have measured the positions of a number of 
band-heads, analyzed the rotational structure in 
so far as possible, determined v’ values by use 


of the isotope effect and given an estimate of the 
convergence limit. They have found that the 
higher members of the band show predissociation 
phenomena. Their results are summarized in 
Table I. The data for v’ = 1, 2 and 3 are obtained 


TABLE I. Constants for bands of system II. 


ve(cm) B,* Remarks 

17,093 Lines sharp 

17,282.5 0.0829 Lines sharp 

17,446.3 0.0800 Lines sharp J>50, 

diffuse J>50 
4 17,580 0.0748 Most lines diffuse 
Convergence 17,864 See Brown a Gibson, 
p. 5 


* B. is in such units that energies will be obtained in 
cm7}, 


directly, whereas that for v’ =4, I have obtained 
by indirect means to be described. 

No regular band-heads or band structure are 
observed for v'=4 or greater. Instead there are 
isolated groups of diffuse and sharp lines, which 
come in pairs, the pairing being due to the fact 
that there are P and R branches present. Each 
group consists of a central sharp line, surrounded 


by others of increasing diffuseness as one goes 
away from the sharp central lines. The lines of 
three such groups have been measured and by 
use of the combination differences their rota- 
tional quantum numbers have been obtained. In 
addition, there is a background consisting of a 
series of very faint diffuse absorption maxima 
which are referred to as System III. The super- 
ficial resemblance of these phenomena to those 
considered in §§2 and 3, is at once apparent. The 
groups of lines presumably correspond to lines 
such as shown in Fig. 3, while, though this is 
more doubtful, one may suppose that the faint 
maxima correspond to the maxima of Fig. 2b; 
there would possibly be enough states piled 
together at a band head to make these visible." 

Brown found by his analysis that the sharp 
energy levels of the two most prominent groups 
of sharp and diffuse lines, which were designated 
as Series I and Series II, had rotational quantum 
numbers J’ equal to 28 and 46, respectively. 
Both Series I and Series II presumably have 
v’’=1, v’=4. Now with any reasonable slopes 
for the two curves V; and V2 it is utterly impos- 
sible to account for two sharp energy levels 
appearing so close together, the periodicity given 
by the theory being of a different order of mag- 


1 Brown had a tentative explanation of these pre- 
dissociation phenomena. He noted that the faint maxima 
seemed to be arranged as though they were an independent 
series of band-heads, which he called System III. The 
potential energy curve for these band-heads he supposed 
to be the adiabatic curve formed by the crossing of a 
repulsive curve and the curve of System II. He supposed 
the sharp lines to belong at the same time to Systein III 
and System II, as their placing seemed to allow that 
assumption. He supposed the spacing of the lines to be 
that characteristic of the rotational levels of System III. 
From this he was able to get an apparent position of the 
minimum of the potential curve for System III, and the 
approximate point at which it was crossed by the potential 
curve. Van Vleck, however (see Brown and Gibson, p. 
542) preferred to consider the sharp lines as vestigial 
remains of System II, which is more in accord with the 
view we take here. One obvious objection to any literal 
acceptance of Brown's ideas and one which he apparently 
did not realize existed, is the fact that the crossing of the 
curves occurs at a point (see Brown and Gibson’s Fig. 4) 
which is nowhere near the Morse curve for System II. 
(See our Fig. 4 and reference 22.) It must be admitted that 
our theory runs into a somewhat similar, but I believe less 
pronounced difficulty. In any event, Brown's suggestion, 
while containing a certain element of truth, had no absolute 
theoretical significance, as he himself well realized. 
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nitude. It is also apparently impossible to 
account for the phenomena by changing the 
vibrational numbering assuming, e.g., that 
v’’=2 for one of the sharp lines and its sur- 
rounding group of diffuse lines, as this pushes 
the energy level beyond the dissociation limit 
of the group of bands, which appears from the 
analysis to be pretty well fixed. The only reason- 
able explanation I have been able to find lies in 
assuming that we have here a new and very 
curious kind of isotope effect, namely that 
Series II is due to the isotopic molecule® ICI*’. 
If this assumption is made, it is necessary to 
compare the combination differences (A:F’’) 
with the combination differences for the molecule 
ICl*7. These can be found from those given in 
Table I of Brown and Gibson’s paper by a short 
extrapolation, and a correction to take into 
account the difference in reduced mass. It is 
found that the agreement between the expected 
combination differences and those actually 
observed is good" if one assigns to the sharp line 
the value J’ =48 instead of 46. Since one knows 
that the value of B,. for v’’=1 for the ICI*® 
molecule is 0.1132 and the value for B,,, for the 
ICI*? molecule is, therefore, very approximately 
equal to 0.9576X0.1132 or 0.1084 (the value 
0.9576 being the ratio of the reduced mass of 
ICl** to that of ICI*7) one can, from the data 
given by Brown and Gibson, calculate the actual 
energies of the sharp and diffuse levels by adding 
to the observed frequency the proper value, 
By J" +1), and, in the case of Series II, 


18 This conclusion depends of course on all the assump- 
tions we have made. It is not claimed that we have here 
made an exhaustive analysis, but certain support is lent by 
other observations. Thus Brown and Gibson state, p. 536, 
that it was impossible to observe the isotope shift of Series 
I and II, from which we may infer that no ordinary isotope 
effect is to be found. And I am indebted to Professor 
Gibson for the information that the intensity of the lines of 
Series I is stronger than that of the lines of Series II, and 
that judging by eye they might easily be in the isotope 
ratio. Of course, the Boltzmann factor would tend to make 
Series II weaker than Series I, but this should be about 
compensated by the increasing multiplicity of the lines with 
larger J values. 

19 This is also indicated by the fact that when one 
calculates the energy of rotational lines from the P and R 
lines using the value of B,»- for the molecule ICl*’ one gets 
values which agree with each other in all cases to within 
0.04 
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subtracting 12.36, to allow for the isotope dis- 
placement of the lower band-head.”° This will 
give the energy referred to the J” =0 level of the 
v’’ =1 vibrational level. In most cases two values 
for a given energy have been obtained, and the 
average of the two has been used. The results 
of this calculation have been tabulated in the 
second column of Table II. If now we assume the 


TABLE II. Energy levels for series I and II and calculation of 
AEa/(Ea- Ea,). 


Series I, J’ = 25-32, series II, J’ =46-—50. 


Obs. 


Energy displ. Ea—Ea _, 

25uu 17,634.26 —8.18 3.94 1.325 1.97 
26un 636.98 —5.46 2.77 0.900 2.08 
27 639.62 —2.82 1.37 0.458 1.99 
28s 642.44 0. 

29 645.48 3.04 -—1.30 0.474 1.74 
30u 648.70 6.26 —2.57 0.965 1.66 
31uu 652.08 9.64 —3.82 1.47 1.60 
32 655.66 13.22 —5.03 2.00 1.52 
46uu 17,729.35 —8.99 5.22 1.55 2.37 
47u 733.82 -—4.52 2.66 0.785 2.39 
48s 738.34 0. 

49u 742.99 4.65 -—2.68 0.802 2.34 
50un 747.96 9.62 —5.19 1.62 2.20 


s=sharp, u=broad, uu=very broad. 


sharp line to be unshifted in position, we can 
calculate the unperturbed energy »,, for v’ =4, 
provided we can make an estimate of B,-. From 
the values of the latter quantity for v’=2 and 
v’ =3, I estimate that for v’ =4 it will be about 
0.0765 for ICI** and this will give 0.0733 for ICI”, 
neglecting the effect of the isotope shift of the 
band-head. We thus find for v’=4, J’=0, ICI*® 
the value 17,580.3 cm~ and for v’ =4, J’ =0, ICI" 
the value 17,566.2 cm™ giving an isotope shift 
of —14.1. We can calculate the isotope shift 
from the formula w,:(p—1)(v’+4) where is 
the frequency of the v’ vibrational level and 
p—1=(M3;— M37) /(M35+ M37) the M’s being 
the reduced masses. w,’, judging from Table I, 
is about 120. This gives for the isotope shift 
—11.7. The discrepancy would be increased a 
little by taking into account the effect of the 
displacement of the band-head on for 
this isotope, as this would tend to increase 3,’ 


20 | cannot check the value of 12.19 given by Brown and 
Gibson, but the difference is not of much significance. 
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a little; nevertheless, the agreement seems to be 
reasonably good, considering the uncertainties 
involved. Thus if we used the, to be sure, rather 
low value of B,-=0.073 for ICI*, to which 
Brown was actually led by the assumption that 
both groups of lines belonged to ICI**, we get 
for the isotope shift —9.3, already on the other 
side of the calculated value. The value for B,, 
which I have tabulated in Table I lies between 
these extremes. And it must not be forgotten 
that the assumption, which we have made (and 
which Brown also made), that the sharp lines 
are unshifted in position, is based on the approx- 
imate calculations of §2; but it is, in any event, 
probably safe to assume that the two sharp 
levels are displaced approximately the same 
amount, which would leave the isotope displace- 
ment the same. 

If, however, we continue to assume that the 
lines are unshifted, we are now in a position to 
calculate the vaiue of 7 which we have called 7 at 
which the intersection of the two potential energy 
curves occurs. For it seems probable that the 
sharp line in each case corresponds roughly to 
the first zero* in where = Ea. 
If this is so, then the energy of the sharp level of 
Series I, minus the energy of the intersection of 
the potential curves should be equal to the same 
quantity for the sharp level of Series II. The 
energy of intersection will be the same for the 
two isotopes for J’=0, but needs to be corrected 
for given values of J’ by an amount J’(J’+1)B; 
where B; = h?/82?Mr,? and is different for the two 
isotopes. If we let B; be the value for the molecule 
IC|*, then the difference in the correction for the 
two cases will be 48 49x0.958 B;—28 X29 B,, 
and this must be equal to the difference in 
energy of the two sharp levels, which is 95.9. 
We thus find that B;=0.0666 and 3.03. 

We can derive some further information on 
the assumption that the sharp levels of Series I 
and Series II correspond to the first zeros of 
Acos(X{ for the respective isotopic mole- 
cules. The energy for the level J’ = 28 is 17,642.4 
cm™ and when corrected by the subtraction of 
the amount J’(J’+1)B; or 28X29X0.0666 it 
becomes 17,588.4. Now Brown states that in the 


* This amounts to saying that in the case of each isotope 
the unperturbed energy of the sharp line corresponds to 
the energy we have called Eq, in §2. 


385 


v’=3 band the rotation levels can be followed 
up to about J’=52 and in this neighborhood the 
levels gradually become diffuse. It is difficult to 
get any comparison between the breadth of 
these levels and those of Series I and II from 
Brown and Gibson’s paper. Likewise, it is dif- 
ficult to know the values of Acos(A {7,(max) — ) 
for the individual levels of Series I and II. If 
these two pieces of information were at hand we 
would be able, from (32), to get Acos(A{7,(max) 
—ra}) and hence \}7,(max)—r.a} for the v’=3, 
J’=52 level. I make the assumption that the 
broadening of the level at v’ =3, J’=52 is about 
the same as that of the broadest of the lines of 
Series I and II yet visible, and it seems quite 
likely that A.os(A{7n(max)—ra}) is still quite 
small for all the lines of Series I and II. The 
value of A.os(A{7n—1a}) has become pretty small 
for \(r,—1a) =4 and we take this as the value of 
A{r,(max)—ra} for =3, J’=52. This would 
make the range AX(r,(max)—7,) of the quantity 
between the level v’=3, J’=52 
and the first sharp line of Series I or II equal to 
about 7.5. Now the energy of v’ =3, J’=52 can 
be calculated by adding the quantity J’(J’+1)B.. 
to the energy of the band head; it must then be 
corrected by subtracting J’(J’+1)B;; the final 
result is 17,483.2cm~'. Therefore the range 
— ra) =4 to A(7,(max) — ra) = —3.5 cor- 
responds to a range in energy of about 17,483 to 
17,588 Now since and 
ra=tot+Ea/V2' we can calculate the quantity 
1/| Vi'|+1/|V2’| from the relation Ad(r,(max) 
—ra)=7.5, using Eq. (26), and putting in the 
numerical values.27 We thus find V4’! 
+1/| V2’| =0.0086. The units of V;’ and V2’ 
are cm™! per A. Since the slope of the Morse 
curve” for V2 is around 1200 cm™ per A in the 
neighborhood of the crossing point of the poten- 
tial energy curves, it is seen that this gives for 
V1’ a value of about — 125 cm™ per A. Since the 


atomic levels of iodine and chlorine are such that 


21 EF, can, in each case, be taken as equal to E, (max) for 
these purposes. Remember, also, V;’ and V2’ have opposite 
signs. 

22 Morse, Phys. Rev. 34, 57 (1929). The Morse curve for 
V2 constructed from the data of Table I is given by the 
expression V2=16,758+1106(1—e7*-7%"-"-))? and is shown 
in Fig. 4. The value of r., calculated by extrapolation of the 
B, values with the aid of the work of Dunham, Phys. Rev. 
41, 721 (1932), is found to be 2.62A. 
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the curve V; must be approaching an asymptotic 
value some 500 cm less than the energy of 
intersection,” it appears that this value is rather 
unreasonably small. 

A rather more serious difficulty appears when 
we attempt to use Eq. (30) to obtain the value 
of | Vi'|/| Ve’|. For the value found is incon- 
sistent with the above conclusion. This is shown 
in Table II. In compiling this table we have 
taken the energy of the sharp level as Eq, in 
both Series I and Series II. In Column 3 is 
given the observed difference of energy between 
the given level and the sharp level (Obs. Displ.) 
and in Column 4 is given AE, which is obtained 
from Column 3 by subtracting from it [J’(J’+1) 
—Jo(Jo+1)]B,, where Jo is the J’ value for the 
sharp level, this being the expected difference in 
energy if the lines were not displaced. To get 
Ea—Ea,, corrected for the change in energy of 
the point of intersection, we may use the obvious 
relation 


Ea— Ea +1) — Jo(Jo+1) Bi); 


the value of E,—E., is given in Column 5. 
‘Column 6 gives the value of AE,/(Ea—£a,)—1 
which from (30) is a minimum value for 
|Vi'|/| V2'|. This turns out not to be quite 
constant, but this lack of constancy would seem 
to be a comparatively minor difficulty. The 
difference in the values of this ratio for Series I 
and Series II may be due to the actual difference 
in this ratio due to the change in slope of the 
curves V; and V2 when the rotational term is 
added. In making the calculations of Table II 
I have used the values B,=0.0748 and 
B;=0.0666, ICl*, with the corresponding values 
for ICI”. The value of B,, is chosen so as to 
make the isotope effect for the vibration level 
v’ =4 come out about right. 

In seeking to explain the discrepancies between 
the values of | Vi’|/| V2’| calculated in different 
ways, a number of possibilities come to mind. 
In the first place, it may be that the curve V2 
begins to deviate somewhat from a Morse curve 
around v’ =4, and has a somewhat smaller slope. 
This is suggested by the low value of B,, for 
v’=4, and the rather low energy of the band- 
head which we have found, in spite of the fact 


3 This is shown in Fig. 4. See Brown and Gibson, p. 537. 
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that the values of B,- for v’=2 and 3 have just 
about the relation to each other that one might 
expect for a Morse curve.”* It should be borne in 
mind that rather sudden changes can occur in 
the law of force between two atoms in a given 
electronic state, and this, in fact, is known to 
occur in another potential energy curve” of ICI. 
It hardly appears, however, that we can account 
in this way for all of the discrepancy, though it 
might be part of the explanation. Another pos- 


sible factor is a breaking down of the assumptions 


that v12 is independent of r, and that the curves 
can be taken as straight lines near the point of 
intersection. Thus one would expect the value 
of V,’ as determined in Table II to be rather 
the value of V/)’ at an energy near Ey, while the 
other method of determining V,’ might be 
expected to give a sort of average value between 
the point of intersection and E,. So if the curve 
V, were quite flat in the neighborhood of the 
point of intersection and then bent up rather 
suddenly at somewhat smaller values of 7, as in 
Fig. 4, this might produce the effect observed. 
Or possibly the flat curve V; intersects the other 
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Fic. 4. Morse curve and hypothetical repulsive curve for 
iodine chloride. 


24 The value of 7% is also about what one might expect 
from the Morse curve, but it may well be too small, due to 
the approximations we have made in calculating it. 

*5 Curtis and Darbyshire, Trans. Faraday Soc. 27, 85 
(1931). 
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side of the curve V2. This would certainly 
produce complications of an unpredictable sort. 
It is also possible that the assumption that only 
one discrete energy level need be considered at 
a time is not quite correct. Finally, it is not at all 
certain but that we would do better in this 
instance if we used the adiabatic non-crossing 
curves to define our unperturbed states rather 
than the crossing curves. An attempt will be 
made to find approximate solutions, using the 
adiabatic curves as the unperturbed curves, in 
order to see whether this assumption fits in 
better with the iodine chloride data. 

Before concluding, mention should be made of 
certain other sharp and diffuse lines which do 
not go with the regular bands and which do not 
belong to Series I or II. Particularly noticeable 
is another group of lines called Series III by 
Brown and Gibson for which apparently v’=5, 
and which is remarkable in that it contains many 
more lines than Series I or Series II. It does not 


387 


seem to be worth while to attempt any detailed 
discussion of Series III. 

We have already noted the existence of faint 
absorption maxima in the same region of the 
spectrum as Series I or Series II, and have stated 
that it is tempting to suppose that they have 
something to do with the secondary maxima of 
Fig. 2b. It must be said however that their 
spacing is much closer than would be expected. 
One of these diffuse maxima, which occurs at 
17,577 cm~, Brown and Gibson identify with the 
band-head for the v’ =4, v’’=1 band. This how- 
ever, does not seem reasonable, as from Table II 
we may judge that the head of this band would 
be displaced to around 17,600. 

In conclusion, I think it may be said that while 
it does not seem possible to account quanti- 
tatively for the predissociation phenomena in 
iodine chloride, progress has been made in 
understanding and interpreting them, and afford- 
ing a basis for their discussion. 


APPENDIX I—EVALUATION OF INTEGRAL, EQ. (22) 


If we denote the quantity |3y?V;’|} by Cx, then we may write from (20) (omitting the normalizing factors) 


F.Fidr= {dr ds {at cos Ci(r—r.)z ] cos [+ C2(r—ra)t]. 


We may extend the integrations with respect to z and ¢ from — « to ~, providing we divide through twice by 2. Doing 
this, and making a double application of Formula 596 of Peirce’s Short Table of Integrals, the expression becomes 


ds di{cos [(C2t— Ciz)r] cos [22 +8+Cir.s— Crrat] 


+4 cos [28—#— Cort + Corat]— cos [22 Cire— Corat — Cort +Cirz]}. 


It may be readily seen, by substituting —¢ for ¢ and —r for r in one of them, that the last two terms cancel in the inte- 
gration. We now perform the integration with respect to r, but instead of integrating from — * to ~, which gives an 
indeterminate result, we integrate between very large but arbitrary limits, let us say*® —r, and r;. This gives us 


We now change over to the variables C.t— C\z and t+z. The Jacobian of this transformation is 1/(C;+ C2). If we integrate 
with respect to C2t—C,z, holding t+z constant, it will be noted that if r; is great enough all the contribution to the 


integral will come from the point where C2t—C,z is zero. We can thus use the known formula ii po sin mx dx =r, setting 


the variables in the cosine function equal to the values they have at Cot— C,z=0. If we let ¢+z=y, then, when Cot—Ciz=0, 
we have 23+£+ Cyr,2—Corat = + C2) 3+ CG). If we make the further substitution 
#=9(C+C,*)(C:+C2)-%, we get for the final value of the integral (remembering the integrand must be multiplied by the 
Jacobian of the transformation) 


[™ ds cos 
which substituted into (22) gives an expression which may readily be seen to be equivalent to (23). 


**It is convenient, but by no means necessary, to have the positive and negative limits the same. 
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APPENDIX II—EvALUATION OF INTEGRAL, EQ. (24) 


We make the substitution r,,—7,=x, and use the definition of \ (Eq. 26). It is seen that (E,— E,,)~d(E,— E,) =dx/x and 
substituting from (23) into the left-hand side of (24) we have for the integral (omitting the constant factor from the 


portion of (23) in front of the integral sign) 


J ds J, dt dx x7 cos cos rae]. 


A trigonometric transformation gives 


at dx x7) cos [23 +t) —Ax(z +4) ]+ cos —Ax(s—2)] 
dx x cos [23 17a) (2 +2) —Ax(s+2)]. 


Another trigonometric transformation gives 


a [as at ax x cos [23 +2) ] cos [Ax(z+é)]+ sin sin [Ax(s+2)], 


which may be readily integrated, yielding (2/2) J. ds sin +2) ]. When the integrand is expanded 


as the sine of the sum of two angles, 2°+2(r,—7a)z and #+(r,—ra)t, it yields one term which vanishes in the integration 
because the integration from 0 to ~ cancels that from — ~ to 0, while the other part gives 


Airy™ has tabulated valuesof the integral J, cos (23 — xz)dz 


(in a slightly different form) over a range of values of x 
from —7.57 to 7.57. This is ample for our purposes for 
negative values of x. On the positive side we need to go 
further, but this can be done by means of Kramers’ 
asymptotic formula, 


J. (2!—x2)dz = cos 


The value of J, ” sin (z8—xz)dz can be obtained for small 


values of x in the same way that Airy obtained the values 
for the other integral, while for large values of x an 
asymptotic value of it may be obtained by Kramers’ 
method." Since Kramers gave but a bare outline of the way 
his results were obtained, we shall give the calculation for 
the sine integral in some detail. 

We make use of contours in the complex plane as 
indicated in Fig. 5. The lines OP and OR make angles of 
4/3 and —7z/3, respectively, with the real axis, and the 
contours W;, W2, Ws; approach asymptotically the lines 
OP, OR, and OQ, the real axis, as shown. We shall indicate 


integrals by subscripts. Thus I, p Means the integral taken 


from 0 to ~ along OP, f. _ means the integral taken be- 
1 


tween infinite limits along W, in the direction shown by the 
arrow, etc. 
Now it is readily seen that 


2.” sin exp 
exp (23+ xz)dz. 


By the use of Jordan’s lemma we may transform the 
integrals from 0 to ix and from 0 to —i~= to integrals 


de dt sin cos [2+A(rn— rae]. 


JJJ—NUMERICAL EVALUATION OF THE INTEGRALS 


along the lines OP and OR, respectively. Making the 
substitution y=z exp (—iz/3) or y=z exp (ix/3) in the 
result, we have 


exp (2+ 2)dz =exp (+i7/3) 
xf exp [—y'+xy exp (+i7/3) jdy. 


For small values of x the integral on the right may be 
evaluated by expanding the factor exp (xy exp (+i73)) asa 
series in x and integrating term by term. One gets in this 
way a series similar to the one used by Airy to evaluate the 
cosine integral. It is extremely tedious to use, however, for 
values of greater than 3, and for such values the asymptotic 
form may be used. 
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The asymptotic form will now be obtained. We note that, 
by Cauchy’s theorem 


exp xz)dz 


The left-hand side of this equation is, as we have noted, the 
quantity we wish to evaluate. It may be obtained from the 
integrals on the right-hand side, as follows. If x is large 
f . and f w, ™ay be evaluated by the method of 


steepest descents. Setting the derivative of 2°+-xz equal to 

zero we see that a saddle point occurs at z=i(x/3)! for 

J and at z= —i(x/3)} for Expanding by Taylor’s 
2 

theorem in the neighborhood of the first of these points we 

see that close to the saddle point we may write 


x2 = 


and it is apparent that the curve of steepest descent for 
exp (z+ xz) will have a slope of 7/4. Integrating along the 
curve we get 


lim exp (2+ xz)dz 
tim Sor p (2-+xz) 


= (3x)-* exp 
Similarly 


lim exp (23+ xz)dz 
sim Sir p (2°+xz) 


exp 
Therefore 


lim J, exp (23+-xz)dz 


lim sin (s*—x2)ds 


= sin [2(x/3)!— +f,” exp (@°+x2)de. 


Now toevaluate J er exp (z*+-xz)dz for small values of x 


we could expand exp xz as a series and integrate term by 
term, but this series will also be difficult to use for x >3. If 
x is very large it is obvious that the integral approaches the 


value — J, > exp (—xz)dz = —1/x. This suggests that we ex- 


pand exp 2° as a series and integrate term by term, but if we 
do this we get the asymptotic non-convergent series, —1/x 
+--+. When x=3, the value of 
1/x is 0.333 and the value obtained from the first two 
terms is 0.259, while the value obtained from the other 
expansion for the integral is 0.296. When x=4, we have 
1/x=0.250 and 1/x—3!/x4=0.227, while when x=6, we 
have 1/x=0.167 and 1/x—3!/x'=0.162. At x=6, then, 
~1/x will be a reasonably good value for the integral, 
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TABLE III. Values of J, ” cos (23 — xz)dz and J, ” sin (z3— xz)dz. 


| | | | 
& & & & 
wn n 
x = x 
—4.0 0.020 5.7 —0.240 0.67 
+3.5 0.038 6.0 0.121 0.69 
—3.0 0.067 0.354 6.3 0.463 0.55 
—2.5 0.112 0.408 6.6 0.720 0.28 
—2.0 0.182 0.461 6.9 0.827 —0.07 
—1.5 0.280 0.507 7.2 0.767 —0.43 
—1.0 0.419 0.534 7.5 0.548 —0.735 
—0.5 0.581 0.521 7.8 0.202 —0.908 
0.0 0.773 0.446 8.1 —0.184 —0.900 
0.5 0.963 0.289 84  —0.528 —0.708 
1.0 1.110 0.039 8.7 —0.742 -—0.370 
1.169 —0.296 90 —0.765 0.030 
2.0 1.078 —0.674 9.3 —0.587 0.393 
25 0.811 —1.021 9.6 —0.250 0.619 
3.0 0.384 —1.239 9.9 0.158 0.642 
3.3 0.078  —1.27 10.2 0.518 0.450 
3.6 —0.237 —1.21 10.5 0.724 0.094 
3.9 —0.530 —1.04 10.8 0.707 —0.320 
4.2  -—0.756 -—0.78 0.469 —0.660 
45 -—0.889 —0.47 11.4 0.080 —0.817 
48 —0.900 —0.12 11.7 —0.337 —0.730 
5.1 —0.788 0.24 12.0 —0.637 —0.426 
54 0.560 0.50 


especially as it is, for large values of x, essentially a 
correction term. For values of x between 3 and 6 we have 
simply made an interpolation. 


The asymptotic formulas for J ” sin (z3— xz)dz and 


J. ” cos (z’— xz)dz hold well for values of x greater than 3. 


Over the range (for x>3) in which Airy has made his 
calculations we have calculated sin [2(x/ 
from his values by making use of the trigonometric relation 
between sine and cosine. This gives bad values when the 
sine is nearly zero, but one can interpolate through those 
points, and, in general, it seems to give better results than 
direct use of the asymptotic value. Outside the range of 
Airy’s calculations the asymptotic formulas were used. I 
have tabulated the values of the two integrals in Table 
III. These values are not highly accurate, but are prob- 
ably good enough for most purposes. The error may 
possibly be as great as 0.03 or 0.04 in the case of the 
sine integral for the worst part of the range where the 
asymptotic approximations were used; where the series or 
interpolation from Airy’s values were used the errors in the 
worst cases should be under 0.01, and this is probably the 
upper limit of error for any part of the cosine curve. 
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The Deflection of Molecular Rays in an Electric Field: The Electric Moment of 
Hydrogen Chloride 


.IMMANUEL ESTERMANN AND RONALD G. J. FRASER* 
(Received March 14, 1933) 


A comparison of the method for determining electric 
moments of molecules by deflection of molecular beams by 
inhomogeneous electric fields with the dielectric constant 
method shows that the molecular beam method has certain 
advantages. This method determines the electric moment of 
gases directly and is able to detect the effect of higher 
rotational states of non-gyroscopic molecules. Moreover, a 
deviation of the axis of the dipole moment from that of the 
axis of rotation by other than 90° should be detected. The 
additional energy due to an electric field is proportional to 


the square of the field intensity and is a function of the 
quantum numbers m and j, if the electric moment is 
perpendicular to the axis of rotation. Otherwise there will 
be a linear dependence on the intensity of the electric field. 
An experimental method of the Rabi type for studying the 
HCI molecule is described. The patterns produced by the 
undeflected and deflected molecules are given and from 
these patterns the electric moment of the HCI molecule is 
found to be approximately 1.95 X10-* e.s.u. Editor. 


I. INTRODUCTION. SCOPE OF METHOD 


URING the last decade, the method of 
molecular rays has been notably developed 
as a tool in research, especially in the Institute 
of Physical Chemistry at Hamburg. The method 
aims at the production of directed, collision-free 
beams of neutral molecules, moving with thermal 
velocities in vacuo; and the study of, for example, 
their scattering in gases, reflection and diffraction 
at solid surfaces, or deflection in magnetic and 
electric fields. The essential feature of the 
method is its extreme directness; for, on the one 
hand, the unidirectional nature of the beams 
usually makes the effect of given external con- 
ditions proximately interpretable, while on the 
other hand their collision-free character elimin- 
ates the disturbing factor of molecular inter- 
actions. It is these features in particular which 
often make the method of peculiar value in the 
study of problems of molecular structure. 
The application of the method which concerns 
us here is to the behavior of molecules in an 
electric field. An electric dipole in an inhomo- 


* The experiments described in this paper were begun by 
one of us (R. F.) while an 1851 Exhibitioner, in the 
Institute of Physical Chemistry at Hamburg (1928-29); 
and continued by the other (I. E.) in the Department of 
Chemistry of the University of California at Berkeley 
(1931-32), while holding a Fellowship of the Rockefeller 
Foundation. 


geneous electric field is acted on by a pondero- 
motive force; hence if a ray of dipole molecules 
SID, selected by the image aperture J from the 
bundle of rays emerging from the source aperture 
S, is directed through an inhomogeneous electric 
field F, it suffers a deflection, the extent and 
character of which, as determined by the detector 
D, allow a determination of the molecular dipole 
moment, and in certain cases the deduction of 
details of molecular structure. (Fig. 1.) 


Fie. 1. 


The standard method of determining molecular 
electric dipole moments, due to Debye,' depends 
on measurements of the dielectric constant of the 
substance, either as a gas, or dissolved in a non- 
polar solvent. The accuracy of the dielectric 
constant method is in most cases far in excess of 
that of the method of molecular rays at the 
present stage of its technical development. On 
the other hand, the deduction of the value of the 
dipole moment of the individual molecules from 
measurements of the dielectric constant of the 


1cf. P. Debye, Polar Molecules, New York, 1929. 
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gas or solution in bulk involves the use of often 
complicated statistical reasoning; which must 
take account, not only of the macroscopic inter- 
action of the electric medium with the field; but 
also of molecular interactions, about which 
assumptions must be made which are often of 
doubtful validity, especially at high densities of 
the fluid. It is just in the characteristic absence 
of statistics from the interpretation of the results 
of an electric deflection experiment that the 
molecular rays method again shows its special 
heuristic advantages. 

We enumerate the following more practical 
advantages. 

(1) The dielectric constant method is subject 
to the limitations imposed, on the one hand by 
the considerable pressures required for reason- 
ably accurate dielectric constant measurements 
in the gas phase, on the other by the solubility of 
the substance it is desired to examine in a suitable 
nonpolar solvent. It may therefore happen, for a 
given substance, not only that the temperature 
required to obtain the necessary gaseous pressure 
is inconveniently high per se or else leads to 
decomposition of the substance, but also that the 
substance is insoluble in nonpolar solvents. On 
the other hand, the pressures required in the 
source of a molecular beam are at the most of the 
order of a millimeter; and again, the applicability 
of the method is clearly independent of the 
solubility of the substance in a particular solvent. 
Thus the molecular rays method may prove 
applicable to the measurement of dipole moments 
in cases where the dielectric constant method 
fails completely; a notable example is that of 
pentaerythritol and its derivatives.” 

(2) Nevertheless the discovery that sym- 
metrically built molecules may often possess 
large permanent dipole moments has led one to 
recognize that the relative motions, of vibration 
or rotation, of the component parts of a com- 
plicated molecule may exert a decisive influence 
on its dipole moment. In certain cases, moreover, 
a dependence of the dipole moment on the tem- 
perature may be expected.? Now Debye’s theory 


*I. Estermann, Zeits. f. physik. Chemie B2, 287 (1929); 
Estermann and Wohlwill, ibid. 20, 195 (1933). 

* cf. Sack, Naturwiss. 8, 337 (1929); also Wolf and Fuchs’ 
article Sterischer Bau und elektrische Eigenschaften, in 
Freudenberg’s Stereochemie, Leipzig, 1932. 
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of the variation of the dielectric constant with 
temperature, which is the basis of the measure- 
ment of dipole moments by the standard method, 
presupposes a temperature-independent dipole 
moment. Thus a possible temperature variation 
of the dipole moment could reveal itself only as 
a departure of the temperature variation of the 
dielectric constant from that to be expected on 
Debye’s theory. Apart from the fact that in 
general only a small range of temperature is 
available, other complications (e.g., association) 
impose a severe restriction on the range of the 
dielectric constant method, and reliable measure- 
ments of the temperature variation of the dipole 
moment by this method have been reported so 
far for isolated cases only.* On the other hand, 
since the molecular rays method measures the 
dipole moment directly, it should be immediately 
applicable to the study of this question, and 
moreover in particular to just the type of mole- 
cule which would in any case prove difficult of 
attack by the dielectric constant method on the 
grounds of (1) above. So far, no experiments of 
this kind have been reported. 

(3) It will appear in the next section (p. 393) 
that if the angle between the axis of the dipole 
moment of a molecule and the axis of rotation is 
only a few degrees other than a right angle, the 
deflection of a beam of such molecules in an 
inhomogeneous electric field is very markedly 
greater than is the case when the axes are strictly 
perpendicular to one another. Thus a merely 
qualitative inspection of the scale of the deflec- 
tion pattern might suffice to give a decision in 
cases where an intramolecular configuration is 
in doubt. 

(4) It has been shown by the application of 
wave mechanics that in the case of non-gyro- 
scopic diatomic molecules the contribution of the 
permanent dipole moment to the dielectric con- 
stant arises entirely from molecules in the lowest 
energy state 7=0.° In such case, higher rotation 
states escape detection by the dielectric constant 
method. On the other hand, the method of 
molecular rays is capable of detecting, and has 


already detected in a qualitative way, the exist- 


‘ Zahn, Phys. Rev. 40, 291 (1932); Greene and Williams, 
ibid. 42, 119 (1932). 

5 Van Vleck, Phys. Rev. 29, 727; 30, 31 (1927); also 
Debye, reference 1, p. 153. 
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ence of these states.® Further, we shall show (p. 
391) that at low temperatures direct evidence of 
the splitting of the rotational energy levels in an 
electric field (quadratic Stark effect) may be 
expected to follow the application of the molec- 
ular rays method under the appropriate con- 
ditions. The experiments on hydrogen chloride 
described in the third section of this paper are a 
first step in this direction. 

Thus the method of molecular rays offers 
several advantages in the study of the behavior 
of molecules in an electric field. Technical dif- 
ficulties have hitherto made progress in this 
direction relatively slow; but the work which 
has been done, and the results which have been 
obtained, are sufficient to justify further develop- 
ment which shall make it possible to follow up 
some of the promising lines of advance we have 
indicated. We propose therefore in the remainder 
of this paper to outline the principles which must 
govern such a development and to illustrate 
their application to a number of problems which 
have already successfully been attacked, with 
special reference to the behavior of the HCl 
molecule in an electric field. 


II. THEORETICAL 


We shall restrict our discussion in the main to 
a simple model of a polar diatomic molecule, 
considered as a rigid dumb-bell; the conclusions 
arrived at need not therefore be immediately 
applicable to polyatomic molecules. 


1. The time-averaged moment 


If a polar molecule, moment uy, is shot into an 
electric field E, the instantaneous change of 
energy is of amount 


AW=une:£, (1) 


where yz is the component of yu in the direction 
of E; but the molecule executes a uniform tem- 
perature rotation whose period is short compared 
with the time the molecule takes to traverse the 
field in a deflection experiment; hence wz changes 
sign every half period, the time-averaged energy 
change AW vanishes to the first order, and the 


6 Wrede, Zeits. f. Physik 44, 261 (1927); Estermann, 
The Dipole Moment and Chemical Structure, ed. P. Debye, 
London 1931, p. 15 ff. (Leipziger Vortriige, 1929). 


392 I. ESTERMANN AND R. G. J. FRASER 


molecule is unaffected by its passage through the 
field. 

The presence of the field, however, distorts the 
motion of uniform rotation of the dipole, and 
there results a time-averaged moment 7=8-E£, 
proportional to the field strength, in the direction 
of the field. Thus to the second order 


AW =}6E?. (2) 


If the center of gravity of the dipole be taken as 
the origin of an arbitrary system of coordinates 
xyz, then the components of the force acting on 
the dipole are clearly 


F,=7-0E/dx; F,=-0E/dy; (3) 


since 

According to quantum theory, AW in (2) can 
take only certain discrete values determined by 
the rotational quantum number j7=0, 1, 2, ---, 
and the equatorial quantum number m, an 
integer defined by the relation —j<m<j. That 
is B is a function B(j, m) of the quantum numbers 
j and m. The values of 7=8(j, m)-£ for the 
rigid dumb-bell molecule were obtained inde- 
pendently by Mensing and Pauli’ and by 
Kronig.® They are 


1 ( 3m? -t) 

if jx£0, 

= 1802] y2E/h? if j=0,] 


where J is the moment of inertia of the molecule 
and h is Planck’s constant. 

For large values of 7 (high temperature), (4) 
takes the classical form 


wE 
i= (3 cos? ¢—1) = (3 cos? ¢—1), (5) 
Tw? 


where w is the (unperturbed) angular velocity of 
rotation, ¢ is the angle between the axis of 
rotation and the field, and e¢, is the rotational 
energy. This expression was first obtained by 
Kallmann and Reiche®; it can be derived very 
simply by considerations due to Stern.” 


7 Mensing and Pauli, Phys. Zeits. 27, 509 (1926). 

8 Kronig, Proc. Nat. Acad. Sci. 12, 488 (1926). 

* Kallmann and Reiche, Zeits. f. Physik 6, 352 (1921). 
1 See Fraser, Molecular Rays, Cambridge, 1931, p. 156 ff. 
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The result (5) still holds for more complicated 
molecular types, provided that their structural 
symmetry is sufficient to allow their temperature 
motion to be referred to a single axis of rotation; 
and provided also that the direction of the dipole 
moment is normal to this axis. In a complicated 
molecule, however, the latter proviso may not 
hold; in which case it is easy to show that if the 
angle between the axes of dipole moment and 
rotation deviates by only a few degrees from 
normality, then in spite of temperature rotation 
there will remain to the first order a dipole 
moment, independent of the field strength. 
Hence the deflection of a beam of such molecules 
in an inhomogeneous electric field will be .a 
comparatively large first order effect, not a 
second order effect as in the case of the dumb- 
bell model (cf. 1(3), p. 391 above). 


2. Theory of a deflection experiment 

Two successful forms of the deflecting field 
have been devised: (a) the beam is sent parallel 
to a charged wire*"'; the directions of E and 
dE/ds are the same (say 2); F,= F,=0 in (3); 
(b) the beam is sent at an angle between the 
plates of a parallel plate condenser; E and dE/ds 
are at right angles.!* (See Fig. 2.) 


Fic. 2. Experimental arrangements for securing inhomo- 
geneous electric fields. 


In the first case, the deflection s of a single 
molecule at the exit of the field is given by 


(6) 


where / is the length of path in the field, and e, 
is the energy of translation corresponding to the 
velocity v. 


In the second case, 


tan 6/e, (7) 


“ Estermann, Zeits. f. physik. Chemie B1, 161 (1928). 
® Rabi, Zeits. f. Physik 54, 190 (1929). 
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to a first approximation; where @ is the angle of 
incidence of the beam, and EF? is averaged in the 
y direction over the distance y’=/ cos 6. (See 
Fig. 2b.) 

Eqs. (6) and (7) refer to the deflection of a 
single molecule, having a given velocity and 
value of @ (or 8). In order therefore to examine 
more closely the character of the deflection suf- 
fered by a beam of molecules, we must consider 
how the value of v and @ permitted under the 
temperature conditions of the experiment by 
Eqs. (4) or (5) are distributed amongst the 
molecules of the beam. We begin with the case of 
high temperature. 

High temperatures. Substituting the value of 


i (or B) from (5) in (6) and (7), we get 
3 cos? @¢—1 3 cos? 


s= 


where 
C’=(u2/16)-E?-1 tan @ 


are constants of the experiment. Thus the 
deflection s suffered by a given molecule is 
decided by particular values of ¢,, €,, the proba- 
bilities of whose occurrence at high temperature 
are governed by a Maxwell distribution law; 
and of ¢, all values of which are equally probable, 
since the beam is collision-free. 

If «,=e,=kT and ¢=0, the corresponding 
deflection is 


So=2C/(kT)?=2C'/(kT)?, (8) 


for the field dispositions (a) and (b) respectively; 
and the deflection s of any molecule referred to 
So as norm is thus s=a-So. 

It is convenient to calculate the intensity 
distribution in the deflected beam in terms of the 
quantity o=s/so rather than s. Fig. 3 shows the 


| 


——» 


Fic. 3. Calculated intensity distribution in a deflected 
beam. 
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result of such a calculation, assuming a parent 
beam of negligible width; in the figure the 
probability P(c)do for the occurrence of deflec- 
tions lying between o and o+doe, obtained by 
numerical and graphical integration over the 
individual probabilities P(e,)de,,  P(e,)de,, 
P(¢)d¢@, is plotted against 

It should be noted that the molecules are still 
assumed to be rigid in calculating the curve of 
Fig. 3; actually the electron configuration is dis- 
torted by the field, resulting in a polarization 
a-E of the molecule in the field direction, the 
deflection due to which is superimposed on the 
deflection arising from the time-averaged mo- 
ment. 

The intensity distribution in the component 
of the deflection pattern arising from the 
electron polarizability is easily shown to depend 
on the equations 


P,(a)do -1/03-do, | 


(9) 


where P,,(c)do is the probability for the occur- 
rence of deflections between o and o+daz, and So 
is obtained from (6) or (7) with a substituted for 
B and kT for e,. Thus an actual deflection pattern 
is the summation of the curve of Fig. 3 and that 
expressed by (9). It frequently happens, how- 
ever, particularly if the field conditions are 
properly chosen, that the contribution of (9) to 
the summation curve is negligible. 

The feature of the intensity distribution 
revealed in Fig. 3 which is of chief importance 
from the experimental standpoint is the fact 
that the maximum of the curve lies very close 
to the position of the undeflected beam; thus the 
effect observed in a deflection experiment is a 
diminution of intensity at the position of the 
undeflected trace, coupled with a broadening of 
the beam, rather than a lateral deviation of the 
beam as a whole. The way in which this diminu- 
tion of intensity can be made the basis of a 
quantitative evaluation of the dipole moment is 
exemplified for HCI in Part III. 

Low temperature. The deflection suffered by 
a molecule in a given state j, m is obtained by 
substituting the appropriate value of (or 


13 Feyerabend, Hamburger Staatsexamensarbeit; see 
Fraser, reference 10, p. 162 ff. 


from (4) in (6) or (7). In this case there is only 
one continuously variable quantity, namely e,; 
and we define so(j, m) for each state j,m as the 
deflection of a molecule in that state having an 
energy of translation «,=kT. Then for each 
state in turn we can calculate the probability 
P;,(a)do for the occurrence of deflections 
oj, m=S(j, m)/So(j, m) lying between o;, » and 
a;, m+do. Now the probability that a molecule is 
in the state j, m is 


where £; is the energy of rotation in the state /; 
the summation extends over all states, the energy 
differences of the states having the same j but 
different m being neglected. Thus the intensity 
distribution in a deflected beam at low tem- 
perature, containing molecules in a number of 
energy states, is obtained by plotting the 
P;, »(a)do vs. o curves for the individual states, 
properly weighted in accordance with (10), and 
constructing the summation curve. An example 
of such a curve, for the case of HCI at 120°K, 
is given in Fig. 4; for the sake of clearness, the 
individual P;, m(«)do curves are shown for the 
three lowest rotation states only. As in the corre- 
sponding high temperature curve (Fig. 3) a 
parent beam of negligible width has _ been 
assumed, and the ordinates chosen to make the 
area under the summation curve unity. It will 
be seen that the lowest energy state 7=0, which 
alone contributes to the observed dielectric con- 
stant, gives rise to attracted molecules (+e) 
only in a deflection experiment. Thus the repelled 
molecules which have in fact been observed (at 
high temperature) are direct qualitative evidence 
of the presence of the higher rotation states (cf. 
1(4) above). In order to identify the individual 
j, m states, however, it will be necessary to carry 
out the deflection experiment with a velocity 
selector; otherwise as seen in Fig. 4, only a single 
maximum can be expected. 


III. EXPERIMENTAL 


1. Intensity measurements 

It will have been evident, from consideration 
of the character of the deflection patterns of 
Figs. 3 and 4, that in order to evaluate the dipole 
moment, or in fact to make any quantitative 
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Fic. 4. Calculated intensity distribution in a beam at low temperatures taking account of the distribution of molecules 
in different rotational states. 


deductions therefrom, a quantitative method of 
measuring the intensity distribution in the 
deflected beam was absolutely essential. The 
earlier work in this field had all been done using 
the condensation target as detector; hence the 
results could only be regarded at best as semi- 
quantitative. It was necessary to find suitable 
quantitative methods of detection for two chief 
classes of substances: on the one hand, easily 
condensible organic substances, where examples 
of the problems enumerated under I(1) and I(2) 
more particularly arise; on the other, gases and 
volatile substances which it is possible to examine 
at low temperature. For the first class of sub- 
stances, a special form of detector making use of 
the heat of condensation was worked out and 
applied to p-nitranaline; for the second class 
exemplified here by HCl, it was only necessary 
suitably to adapt an already existing method for 


“ Wohlwill, Zeits. f. Physik 80, 67 (1933). 


the detection of gases, namely to allow the beam 
to enter by a slit an otherwise closed vessel con- 
taining a hot wire gauge, and to measure the 
resulting equilibrium pressure.!® 


2. The deflecting field 


It is much easier to obtain comparatively large 
deflections with the charged wire field arrange- 
ment than with the parallel plate condenser. 
Nevertheless, the Rabi field is to be preferred for 
quantitative measurements on the following 
grounds: 

(1) It is impossible except by recourse to 
beams of very small cross section to obtain a 
high inhomogeneity near a charged wire which 
shall at the same time be constant over the 
cross section of the beam. 

(2) It is difficult to estimate the values of E 
and dE /ds (Eq. (6)) near the wire on account of 


16 Knauer and Stern, Zeits. f. Physik 53, 766 (1929). 
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the departure of practical arrangements from the 
form of the ideal cylindrical condenser; on the 
other hand, EF? (Eq. (7)) is readily obtained, if 
the plates are properly designed, from the theory 
of the parallel plate condenser. 

(3) The field near a charged wire varies as the 
inverse cube of the distance from its center; 
hence the contribution of the electron polariza- 
tion a-E to the deflection pattern can easily 
become unduly prominent. 

(4) The convenience of the method of esti- 
mating the dipole moment from the diminution 
of intensity at the position of the undeflected 
beam can actually make a small deviation of the 
maximum of the deflection pattern from that 
position an advantage. 


To HCL Reservoir 


Qqy Ss 


very short male cone, fitting the corresponding 
female seating Gj in the central housing of the 
apparatus. In this way, the apparatus is divided 
into three chambers: the source chamber, the 
collimator chamber, and the observation cham- 
ber. Each chamber is separately evacuated by 
means of a butyl phthalate diffusion pump 


Fic. 5. Experimental arrangement of the apparatus. 
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For these reasons, the parallel plate condenser 
field was chosen for use in measuring the dipole 
moment of hydrogen chloride. 


3. Apparatus 


The apparatus in its final form is shown in 
Fig. 5. The beam system, consisting of the source 
slit Ss; the fore slit Fs; the image slit Js; the 
parallel plate condenser C; and the detecting 
system, a pair of Pirani gauges Pg, is mounted on 
a rail R. The upper surface of the rail is ground 
true, and the slits are adjusted parallel to this 
surface, by means of a rider carrying a knife 
edge, before inserting the beam system into the 
apparatus. 

The image and fore slits are carried each on a 


To hig 


specially designed for high pumping speed."* The 
vacuum envelopes for source chamber and 
observation chamber are of Pyrex, and are con- 
nected to the central housing by the ground 
joints J; and Js, respectively. The central hous- 


16 Estermann and Byck, Rev. Sci. Inst. 3, 482 (1932). 
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ing is of Monel metal, as are all metal parts, this 
material having been found to resist the corrosive 
action of HCI rather well. 

The beam is defined -by the slits Ss and Is, 
which are 0.016 mm and 0.04 mm wide, respec- 
tively. The function of the fore slit Fs is not to 
define the beam, but to step down the vacuum 
from about 5.10-* mm in the source chamber to 
about 10-* mm in the observation chamber, the 
pressure in the central chamber being about 10~° 
mm; its width is therefore comparatively large, 
about a tenth millimeter. The source slit is given 
a wedge-shaped form in order to reduce the 
Knudsen resistance between it and the fore slit, 
the distance between them being only 7.3 mm. 

The detecting system diverges markedly in 
design from that used in previous work on gas 
beams.'® There the detector slit was connected 
via a narrow tube to a glass-walled gauge outside 
the apparatus, the gauge being immersed in 
liquid air to give increased sensitivity. Now HCl 
is strongly adsorbed on glass walls, and moreover 
it is too readily condensible to make cooling 
profitable; therefore the gauge was made of 
Monel metal, and was placed inside the appa- 
ratus. In this way it was possible to make the 
volume of the gauge very small, thus reducing 
the incidence of adsorption. 

The detector slit Ds, in the form of a canal 3 
mm high, 0.04 mm wide, and 1.5 mm deep, 
communicates directly with the gauge; immedi- 
ately above it, but outside the range of the 
beam, an exactly similar canal communicates 
with a second gauge, the function of which is the 
compensation of pressure and temperature vari- 
ations.!® The two gauges are made from the 
solid: two parallel holes are drilled in a block of 
Monel metal; these holes are closed at each end, 
at the front by plugs carrying the canals, at the 
rear by glass insulators Gi carrying each a 
platinum hook. Near the front end of each 
gauge a second glass insulator enters from the 
side, carrying the platinum leads for the gauge 
filament F. The latter, of platinum strip 0.05 mm 
wide, and 0.003 mm thick, is looped double over 
the platinum hook as seen in Fig. 5. 

The block housing the gauges is mounted on 
a dovetailed slide Dt, and can be traversed 
across the beam by means of the screw S 
operated by a screwdriver carried on the ground 
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joint J;. The position of the detector slit Ds can 
be read off on the graduated drum Dr, each 
scale division corresponding to a movement of 
the detector slit of 0.01 mm. 

The gauges form two adjacent arms of a 
Wheatstone bridge; thus stray pressure changes 
occurring in both gauges simultaneously are not 
registered by the galvanometer. A pressure 
change of 10-* mm in the measuring gauge gave 
a galvanometer deflection of 3 cm at a scale 
distance of 130 cm; the galvanometer sensitivity 
was 10-7 volt/mm at one meter. To eliminate 
zero creep, timed readings were taken at inter- 
vals of 15 sec. with the beam alternately allowed 
to enter the detector slit and cut off therefrom 
by the electromagnetically operated shutter Sh, 
working in the central chamber. 

The deflecting field is produced at the entrance 
of the parallel plate condenser C. One plate is 
mounted, with provision for adjustment, directly 
on the rail R, which is earthed via the central 
housing; the other plate is carried on the first 
by a glass U, and is connected by the lead L toa 
high-voltage source consisting of a transformer 
and one wave rectifying set (Kenotron and 
Leyden jar). The plates are 40 mm by 10 mm, 
with a gap of 2.3 mm. The angle between the 
edges of the plates and the beam was 7°. The 
edges are carefully rounded to follow in outline 
as far as possible the equipotentials at the 
entrance of an ideal parallel plate condenser; in 
this way it was possible to calculate the average 
field E? with confidence from the theory of the 
parallel plate condenser. 

The gas was fed to the source slit, which can 
be cooled if necessary by the Dewar Dw through 
a heavy flexible lead of copper Cc, from a 
reservoir in which the pressure of HCI was held 
constant at about 2 mm, this reservoir being 
fed in turn from a high-pressure reservoir 
through an adjustable leak. Before starting a 
run, however, it was found necessary to fill the 
apparatus with HCI at some centimeters pressure 
for several hours, in order first to saturate the 
walls of the gauges with the gas. 

The intensity distribution in the beam was 
measured first without field, then with field, and 
finally again without field. Since the intensity in 
the parent beam is not absolutely constant, the 
intensity at a given point in the deflected beam 
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was compared with the average value of the 
intensity at that point observed in the bracketing 
runs without field; but only those complete runs 
in which the parent intensity remained constant 
to within 10 percent were accepted. . 


4. Results 


An example of the results obtained is given 
in Fig. 6, which shows the intensity distribution 
in an undeflected beam, and that in the same 
beam after passage through the electric field. 


1.0 
o without field 4 
x with field on 
30.5 
0.02 0.04cm 


“0.04 -002 
& 


Fic. 6. Experimental distribution of molecules in a de- 
flected and an undeflected beam. 


The dipole moment of the HCl molecule is 
evaluated from the intensity diminution at the 
position of the maximum of the undeflected 
trace. The argument is as follows: The width of 


P (a) 
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the undeflected beam is not negligible in com- 
parison with that of the deflected beam, as was 
assumed in calculating the ideal curve of Figs. 3 
and 4, but has the form shown in Fig. 7. Let 
f(oo) be the equation of the curve representing it; 
then the equation of the deflected beam is clearly 


F(o0) = [ (11) 


where P(c) is the shape of the deflected beam 
arising from a parent beam of negligible width, 
JS P(c)do being normalized to unity. In par- 
ticular the ordinate at the center of the deflected 
beam, assuming it coincident with the center 
oo=0 of the undeflected beam, is 


F(0)= (12) 


and is obtained by multiplying each ordinate of 
the undeflected beam by the appropriate P(c) 
and then integrating. 

In order to evaluate yu, therefore, the abscissae 
of the experimental curve for the undeflected 
beam are reduced to the same abscissae a as the 
curve of Fig. 3 for a certain assumed value of 5» 
(see Fig. 7): the ordinates being reduced so that 
Imax is unity. Multiplication of each ordinate of 
this reduced curve by the appropriate P(c) 
gives the dotted curve of Fig. 7, the area under 
which (cf. Eq. (12)) represents I’max, the inten- 
sity at the maximum of the deflected curve for 


ity| oF the undetlected beam. 


ja 
Ri 
II 
— 
ace 
det 
= 
cul 
\ ims 
10 bee 
| 
Fic. 7. 


the particular value of s» assumed. In this 
way I'max could be plotted as a function of 
so= (u2/8(RT)?)-E*-1tan @ for a given set of 
experimental conditions. 


Oo 
— Umax / Tmax 


/ 


0 1 2 3 4a 5 
— So (x107% em) 


Fic. 8. 


Such a plot under the conditions of the present 
experiments is given in Fig. 8. The value of so 
corresponding to a given value of I’max/Imax iS 
read off from the curve, and yu evaluated. Table I 
gives a summary of the results. 


TABLE I. 
tan 


1=4.72 cm; 6=83°; T=290°; E? =0.6E? e.s.u.; Condenser 
gap d=2.3 mm; V=0.23E e.s.u. 


Run V I'max 

No. kilovolt 104-5, 
| 35 0.95 14 1.69 
II 32 0.97 10 1.56 
Ill 35 0.92 24 2.21 
IV 35 0.93 21 2.07 
V 40 0.90 31 2.20 


Mean value: p= e.s.u.* 


*The limit of error given above does not of course 
account for systemmatic errors which may arise, e.g., in 
determining the form of the curve for the undeflected beam. 


Two assumptions have been made in the cal- 
culations, which must be justified. (1) The max- 
ima of the deflected and undeflected beams have 
been assumed coincident. Now the mean value 
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So from Table I is 2.10-* cm; the value of o at 
the maximum of the probability curve of Fig. 3 
is 0.04; hence the displacement o- 59 of the max- 
imum of the deflected curve from that of the 
undeflected curve is s~ 1y, and can be neglected. 
(2) No account has been taken of the electron 
polarizability. The polarizability a of HCl as 
measured by Zahn!’ is 1.7010-* e.s.u.; with 
the field strengths used a-E~2-10-*; insertion 
of this value in Eqs. (9) shows that the effect 
of the polarizability on the height of the deflected 
maximum is less than 3 percent, and it can 
therefore be neglected. 

Great accuracy is not claimed for the present 
experiments. They should rather be regarded, 
together with the parallel investigation on easily 
condensible substances referred to on p. 395, asa 
first attempt to obtain quantitative data on 
dipole moments from molecular ray experiments. 
It should be readily possible, for example by 
using narrower beams and stronger fields, greatly 
to improve the accuracy of the method." 

It is intended to repeat the experiments under 
improved conditions, and moreover with a cooled 
source slit, in order to realize the case of low 
temperature dealt with in Section II of this 
paper, where the influence of the quantum dis- 
tribution of the rotational states becomes 
noticeable. 

We are glad of this opportunity to express our 
several thanks to Professor O. Stern and Pro- 
fessor G. N. Lewis, who successively gave the 
above work the ready hospitality of their 
respective laboratories. 


17 Zahn, Phys. Rev. 24, 400 (1924). 
18 Owing to the expiration of the Fellowship held by one 
of us (I. E.), the experiments had again to be interrupted; 
but we considered the results so far obtained to be suffi- 
ciently interesting methodically to justify publication. 
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Activation Energies for Reactions of Atoms in Different States 


VOLUME 1 


G. K. ROLLEFSON AND J. C. Potts, Chemical Laboratory, University of California 


A slight modification of Eyring’s method enables the 
activation energy for the reaction of iodine monochloride 
with normal or with excited chlorine atoms to be calculated. 
The calculations indicate that normal chlorine atoms are 
much more reactive with respect to iodine chloride than are 
chlorine atoms excited to the 2P} state. Experimental 
evidence supporting these calculations is obtained from the 


(Received April 4, 1933) 


results of Rollefson and Lindquist with mixtures of 
hydrogen, chlorine and iodine monochloride. 

In the modified method, the atomic interactions—which 
give rise to the activation energy—are obtained from the 
potential energies of those molecular states which may be 
adiabatically formed from atoms in the particular elec- 
tronic states to be considered. 


SHORT time ago Rollefson and Lindquist! 
concluded from their experiments on mix- 
tures of hydrogen, chlorine and iodine mono- 
chloride that the reaction chains yielding hy- 
drogen chloride were started by chlorine atoms 
in the *P; state. Recent work,” especially that 
by Rodebush and Klingelhoefer, has shown that 
the reaction between chlorine atoms and hy- 
drogen has a heat of activation sufficiently large 
that there should be time enough for the two 
lowest states of the chlorine atom to reach 
statistical equilibrium. Therefore, no matter 
what source of chlorine atoms was used, hy- 
drogen chloride formation should be observed. 
Rodebush and Klingelhoefer suggest that the 
results (of the experiments of Rollefson and 
Lindquist) with iodine chloride are due to 
chlorine atoms reacting more readily with some 
other substance than hydrogen. The truth of 
this is apparent since iodine chloride is a powerful 
inhibitor for hydrogen chloride formation. It 
has been suggested! that this may be caused by 
a reaction such as ICI+Cl=Cl.+I. To account 
for the observed phenomena fulfillment of a 
second condition is necessary: normal chlorine 
atoms must be removed more readily than 
excited ones. This second requirement is readily 
fulfilled by postulating that ICI+Cl=Cl.+I is 
very much faster than ICI+CI*=Cl.+I. This 


1 Rollefson and Lindquist, J. Am. Chem. Soc.:52, 2793 
(1930); 53, 1184 (1931). 

2 Rodebush and Klingelhoefer, J. Am. Chem. Soc. 55, 
130 (1933); H. v. Hartel, Zeits. f. physik. Chemie B11, 96 
(1930). 
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may seem to be very improbable since it requires 
an excited atom to be less reactive than the 
corresponding unexcited atom. We shall show 
that considerations based on thermodynamics 
and the quantum theory of activation lead to 
just this conclusion. 

The use of quantum mechanics in the calcu- 
lation of activation energies has been presented 
in great detail in papers by Eyring and 
Polanyi*: ‘; the second reference deals especially 
with the reactions of the halogen atoms. How- 
ever, in Eyring’s calculations, no distinction is 
made between the reactions of halogen atoms in 
the ?P; state and those in the ?P; state. The cal- 
culations are all made for the halogen atoms 
in the ?P; state. 

The approximate values of the interaction 
terms in the equation 


potential energy 


are usually obtained from the potential energy 
functions for the corresponding diatomic mole- 
cules. Now if we consider the reaction ICI+Cl* 
=I+CICI*, the interaction terms are to be 
derived from the energies of the molecules ICI, 
CICI*, and ICI*. If Cl* is in the 2P; state and the 
other two atoms are in the ?P} state, then we 
must obtain our interaction terms from ICI ('2), 


3 Eyring and Polanyi, Zeits. f. physik. Chemie B12, 


279 (1931). 
4 Eyring, J. Am. Chem. Soc. 53, 2537 (1931). 
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Cle (IIo), and ICI The reasons for this 
choice of states are: (1) ICI ('Z), case under 
consideration. (2) Cle (*IIo), a 2P3Cl atom and a 
2P;Cl atom cannot form a normal molecule. 
(3) ICI (Io), atom and a?P;Cl atom cannot 
form a normal moleculé. This treatment permits 
us to differentiate between the reaction of an 
IC] molecule with a chlorine atom in the *P; 
state and one in the ?P; state, for in the latter 
case all the interactions would be obtained from 
the terms. 


20r 


E Ckilocalories) 
T 


1 2 3 4 bs 6 
Nuclear Separation in Angstroms 

Fic. 1. Potential energy diagram for ICI. The states 
marked '¥, *I1,, and Ot are formed by an iodine atom in the 
*Py state and a chlorine atom in the 2P} state; the *IIp state 
is formed by a normal iodine atom (?P3) and an excited 
*(P;) chlorine atom. The diagram for the chlorine molecule 
is similar® except that transitions involving the II, and O+ 
states have never been observed experimentally. 


In order to explain the experimental results of 
Rollefson and Lindquist it is not necessary to 
make the activation energy calculations for both 
of the reactions, 


See Fig. 1. 

* Providing no electronic change takes place as the atoms 
approach. 

"Weitzel, Handbuch der Experimental Physik, Er- 
ganzungswerk I, 383-385 (1931). Akademische Verlags- 
gesellschaft, Leipzig. 


*Mullikan, Phys. Rev. 36, 669, 1440 (1930); Rev. Mod. 
Phys. 4, 17, 70 (1930). 


I2p;Cl +Cl = I:p3+Cl (1) 
+Clep, = (2) 


Reaction (1) is exothermic by an amount equal 
to the difference between the heats of formation 
of IC] and Cle, (from normal atoms), approxi- 
mately 7 kilocalories. We calculate the corre- 
sponding heat of activation to be about 4 kilo- 
calories. To calculate the heat of reaction (2) we 
use the heat of formation of iodine chloride (from 
normal atoms) and the heat of formation of the 
3]Iy state of Cle. This latter’ is only 9.6 kilo- 
calories which makes reaction (2) endothermic 
to the extent of 40 kilocalories. According to the 
first law of thermodynamics it would be necessary 
to supply at least this much energy if the reaction 
is to proceed. Hence this is a minimum value for 
the heat of activation and without further cal- 
culation we may say that reaction (2) is very 
slow compared with reaction (1). Expressed in 
other words, by activating the chlorine atom to 
the ?P; state we have rendered it much less 
reactive toward ICI. This does not mean that 
if we introduce a large number of *P; chlorine 
atoms into some iodine chloride vapor that no 
reaction would occur; we must remember that 
collisions would eventually transfer the chlorine 
atoms to the ?P; state and then the reaction 
would proceed according to (1). From the 
experimental results of Rollefson and Lindquist 
we conclude that such energy transfers have a 
probability sufficiently low so that many of the 
excited atoms survive long enough to react with 
hydrogen. 

At the present time the possibility of extending 
these considerations to other reactions is limited 
by insufficient knowledge of the higher electronic 
states of most molecules. It may be that in those 
cases where there is a marked discrepancy 
between the experimental activation energies 
and those calculated by Eyring that the dis- 
crepancy would be greatly reduced if the calcu- 
lation, in which the higher states are considered, 
could be made. 
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The Vibration Spectra and Structure of the Cyanogen Halides 


W. WEsT AND MARIE FARNSworTH, Department of Chemistry, Washington Square College, New York University 
(Received March 13, 1933) 


The Raman spectra of the cyanogen halides in the liquid 
state or in alcoholic solution contain 3 frequencies , v2, and 
v3. In the order chloride, bromide, iodide, »;= 2201, 2187, 
2158 ve = 729, 580, 470; v3 =397, 368, 321; the order of 
intensities of Raman transitions is »; >v2>v3. The relative 


values of the frequencies indicate a linear arrangement of 
the nuclei in the lowest electronic state of these molecules, 
and the data are in accord with the structure X-—C=N 
rather than X—N=C. 


HE following report contains the results of 

experiments designed to establish the 
fundamental vibration frequencies of the unsym- 
metrical triatomic molecules, the cyanogen 
halides, by investigating the Raman spectra of 
these compounds, and the evidence on the 
debated matter of their structure afforded by a 
knowledge of these frequencies. 


EXPERIMENTAL 


Cyanogen chloride was easily obtained in a 
pure condition by passing chlorine into the 
double cyanide of zinc and potassium according 
to Mauguin and Simon’s modification of Held’s 
method.! The bromide was a commercial product, 
while the iodide was prepared by the action of 
mercuric cyanide on iodine in ether solution.’ 
The chloride and bromide were examined in the 
liquid state. The chloride, whose boiling point 
is 15°C, was condensed in a well cooled scattering 
tube, which was sealed off when filled. The 
spectrum was examined at room temperature 
(18°C). The bromide, which melts at 52° and 
boils at 61°, was examined in a tightly stoppered 
tube 2 inches long, placed in a bath of acetone 
vapor (B.P. 56°). 

The iodide melts only under pressure. There is, 
however, sufficient decomposition of the sub- 
stance at the melting point, about 146°, to color 
the liquid dark red, and it was impossible to 
obtain the spectrum of the pure liquid. The 
substance was investigated dissolved in methyl 
alcohol, in which it is very soluble in the cold 
without undergoing chemical change. 


1 Mauguin and Simon, Ann, Chim. 15, 18 (1921). 
? Linnemann, Annalen 120, 36 (1862). 


The source of light was a Hanovia 220 volt 
quartz mercury arc, focussed by a large glass 
condensing lens on the scattering tube, and the 
spectrograph was a Hilger constant deviation 
instrument with a dispersion of about 100 cm™ 
per mm in the blue and violet. Spectra were 
photographed both with the lines of wave-length 
4047A and the group about 4358A as exciting 
lines and with the former line removed by means 
of a filter composed of two plates of Noviol 0 
glass. A solution of p-nitro toluene in alcohol is 
also very suitable for removing 4047 without 
appreciably dimming 4358, and was used in the 
later experiments. The substances gave relatively 
strong Raman spectra with comparatively little 
background. 


TABLE I, 

CICN BrCN ICN 
Source Source Av cm~!|y Source Av cm” 
22541 e 397 |22572 e 368 |22617 e 321 
22501 a 2202 {22516 a 2187 |22549 a 2156 
22321 b 2195 |22360 e 580 |22468 e 470 
22263 d 732 |20754 e 2186 |20780 e 2158 
22219 e 729 
22129 e 809 
20848 c 2191 
20786 d 2201 
20736 e 2200 


The results are summarized in Tables I and II, 
and microphotometer curves are reproduced in 
Fig. 1. The column headed source refers to the 
line in the Hg spectrum causing the Raman line. 
The wave-lengths of Hg lines concerned are 4, 
4047; b, 4078; c, 4339; d, 4348; e, 4358a. The 
figures in parenthesis after the vibration fre- 
quencies in Table II give a rough estimate of the 


relative intensities of the Raman transitions. The | 
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feeble Raman line in the spectrum of CICN of 
Av 809 cm~ seems most probably to originate in 
a transition involving the oscillation of funda- 
mental frequency 397 in which the vibrational 
quantum number increases by 2. 


TABLE II. Fundamental vibration frequencies of cyanogen 


halides. 
vy; vecm! v3 
CICN 2201 (10) 729 (5) 397 (3) 
BCN 2187 (10) 580 (5) 368 (2) 
ICN 2158 (10) 470 (2) 321 (1) 
DISCUSSION 


Although the cyanogen halides have been 
known for over a century, there is still uncer- 
tainty as to their constitution, the choice being, 
in the symbolism of organic chemistry between 
the structures 


X-C=N and X—N=C= or 
(1) (2) (3) 
Nitrile Carbylamire 


Gutman’ has also proposed’ a cylic formula 


which, however, with trivalent halogen, seems 
inherently improbable and is excluded by the 
results we have obtained. 

The weight of opinion on the side of organic 
chemists seems to incline to the carbylamine 
formula (2), modified by arguments that the 


* Gutman, Ber. 42, 3623 (1909). 
* Chattaway and Wadmore, J. Chem. Soc. 81, 192 (1902). 
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Qu: 


CNBr 


Fic. 1. Raman spectra of cyanogen halides. 


CNI 


chloride is a carbylamine, the iodide a nitrile, 
and the bromide a tautomeric mixture of the 
two.’ The cyanogen halides in fact show the 
chemical behavior neither of typical compounds 
containing a well-defined nitrile group nor yet 
of those containing a nitrogen-halogen linkage— 
they are much more stable than any compound 
which we have been able to find reference to as 
definitely containing a nitrogen halogen linkage. 
Nor, it must be admitted at the outset, do the 
results we have obtained on the fundamental he 
vibrations of the compounds, permit a com- oF 
pletely definitive answer to be made as to their 7 
linkages, though there is no suggestion that the | 
linkages are not of the same nature in all three 4 
compounds. 

In the first place it is to be noted that the 
comprehensive work of Dadieu and Kahlrausch 
on organic C—N compounds shows that the — 
great differences in chemical behavior which exist 
between authentic nitriles and authentic car- 
bylamines is not well expressed by the distinction ‘ Ps. 
between a triple and a double bond between C mide 
and N as indicated in formulae (1) and (2).° oer 
The Raman spectra of the alkyl nitriles show a 
remarkably constant frequency of 2245 cm™ 
which is undoubtedly due to the vibration 
between the C and N atom, and which is about 
the magnitude to be expected for the vibration 
of the triple bonded group C=N. The alkyl 
carbylamines have a different, but not markedly 
different frequency at 2180 cm~; and it is im- 
possible to ascribe such a small change in vibra- 
tion frequency to a change in linkage represented 


5C. V. Zappi, Bull. Soc. Chim. 47, 453, 537 (1930); 49, ; 
397 (1931). Be! 

6 Kahlrausch, Der Smekal-Raman-Effekt, p. 319 et seq. 
(Julius Springer, Berlin, 1931). 
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by the difference between a triple and a double 
bond; the bond C=N would in fact have a 
frequency near 1700 cm~. It is possible to write 
electronic formulae for the two classes of com- 
pounds which are completely analagous. 


H:N:C: and H:C:N: 
carbylamine nitrile 

and these are probably a better representation 
of the state of binding in the unexcited state of 
the 2 molecules than the conventional bond 
formulae. 

Returning to the evidence provided by a 
knowledge of the vibration frequencies of the 
cyanogen halides on the alternative structures, 
one finds that, if not absolutely conclusive, it is 
strongly in favor of the nitrile constitution for all 
three compounds. The carbon-nitrogen frequency 
in the chloride, 2201 cm~, is decidedly higher 
than the value for the aliphatic carbylamines, 
2180 cm-'. True, it is lower than the value 2245 
cm~', characteristic of the aliphatic nitriles, but 
it is likely enough, and in accord with experience 
in similar cases, that substitution of a heavy Cl 
atom for the one or two C atoms in the alkyl 
group which alone influence the C—N vibration 
should lower the value of this frequency, as is 
shown by its progressive diminution in the 
bromide and iodide. The observed value of the 
frequency for CN I is only 12 cm™ lower than 
that for methyl carbylamine whereas substi- 
tution of a heavy iodine atom for the CH; group 
would have been expected to diminish it much 
more considerably, probably by about 80 or 90 
cm~!. In the conclusion that all three halides 
have the same structure we are in agreement 
with that of Badger and Sho-Chow Woo’ who, 
as a result of their study of the continuous ultra- 
violet absorption spectra of these compounds, 
advance cogent reasons based on a discussion of 
binding energies for preferring the nitrile 
structure. 


SHAPE OF THE MOLECULE 


It may first be relevent to refer to the con- 
clusive proof given by Barker and Kya Nam 


7 Badger and Woo, J. Am. Chem. Soc. 53, 2542 (1931). 


M. FARNSWORTH 


Choi’ in their rotational analysis of the infrared 
absorption bands of HCN that this molecule is 
linear. 

In the following discussion we shall assume 
that the linkage is X—C—N where X is the 
halogen. A triatomic molecule will have 3 x 3—6 
=3 degrees of vibrational freedom if bent and 
3xX3-—5=4 degrees if linear. In the latter case, 
two of the frequencies, those corresponding to 
the bending of the molecule in two perpendicular 
planes, are equal if the binding forces are axially 
symmetrical. Three fundamental frequencies will 
therefore exist, all of which appear in the Raman 
spectra in the case of the cyanogen halides. 
Anticipating the result that the frequencies 


observed for these compounds are in accord with | 


a linear molecule, we may represent the vibra- 
tions as follows: 
<_ > v2 
The solution of the mechanica! problem of the 
motion of 3 point masses under the action of 
“valence” forces has been given by Lechner® who 
finds the condition that for a linear arrangement 
of the masses, m1, m2 and m3. 


V3 


Vi ve 
—+—>2k, where k=1— 
me? 
m me Me2M3 
M23 = 
m+mMe Ms 


The values for v;/v2+v2/v; and 2k, respectively, 
are: chloride, 3.36, 2.60; bromide 4.00, 2.74; 
iodide 4.84, 2.80. The observed frequencies are 
compatible with a linear arrangement of the 
atoms. 

It is of some interest to calculate the force 
constants, the restoring force for 1 cm displace- 
ment of the atoms in the various vibrations 
when the latter remain simple harmonic. Accord- 
ing to Lechner, these are given for the unsym- 
metrical linear three-point-mass model by 


ny +n? =f 12/ Miz +f. 23/ 23, 


8 Barker and Choi, Phys. Rev. 42, 777 (1932). 
® Lechner, Sitz. Wien, Akad. Math.-natur. Classe 141, 
2a, 291 (1932). 
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fie fos 
Miz M23 m2? 


ys. 


If v is the frequency of vibration in cm™ and 
m is given in atomic weight units (Cl=35.46, 
etc.), m? is defined as equal to 47°c?y?/N= 
5.863 X 10-*v?, where N is Avogadro’s number. f, 
the restoring force for stretching displacements 
and d, the restoring force for the bending dis- 
placement are then given in dynes per cm.” 

d and yp, are defined by the relation 


Ma Me 23 S23" 


where Sj2, So3 are the distances between m, and 
and mz and m3, respectively, and 
= do3593*=ds*. diz is the restoring force for a 
bending displacement of m, m3; being fixed and 
similarly for 

The values of the force constants are given in 
Table III. The calculation of the bending forces 


TABLE III. 
fir fos d (approx.) 
CICN 5.14105 16.65 X 105 0.17 X 10° 
dynes/cm dynes/cm dynes/cm 
BrCN 4.16 16.84 0.15 
ICN 2.60 16.80 0.11 


requires a knowledge of the interatomic distances 
which are unknown in these compounds. If as a 
first crude approximation the distances S12. and 
Ses are assumed equal, siz=Se3=s, and 1/p4 
=1/u12+2/m2+1/p23, whence the above values 
of d are obtained. 

Comparison of the frequencies v; and v2 in the 
cyanogen halides with those exhibited by other 
compounds containing the C=N and C~ halogen 
linkages shows that these molecules are further 
illustrative of the preservation of “bond fre- 


" Kohlrausch, Der Smekal-Raman-Effekt, pp. 154, 170. 


quencies” in polyatomic molecules. The fre- 
quency »; has approximately the same value, 
about 2200 cm in all three of the halides, and 
all compounds containing the linkage C=N, 
including the diatomic CN molecule known from 
spectroscopy, have a frequency not very different 
from this value. In the same way, the frequency 
ve=729 for cyanogen chloride is close to the 
strong frequency 710 in CH;C1; v2for the bromide 
is nearly the same as the frequency 594 in CH;Br, 
and the value 469 for the iodide is comparable 
with the frequency 522 in CHI. It would there- 
fore be possible to make the statement that the 
vibration spectra of the cyanogen halides are 
partially composed of C=N and C~ halogen 
bond frequencies, and to identify these fre- 
quencies with »; and v2, respectively. Although 
the calculation of the force constants is in terms 
of a dynamical model in which the frequencies 
v, and v2 involve all of the masses, in consequence 
of which it becomes impossible to attach meaning 
to C=N and C~ halogen vibrations, the experi- 
mental fact of the persistence of ‘‘bond vibra- 
tions” in polyatomic molecules nevertheless 
suggests the reasonableness of identifying the 
constant fe3 as the ‘‘force constant for the C=N 
vibration” and the constant fi2 with the ‘‘force 
constant for the C-halogen vibration.’’ In fact, 
the constant feos is very close to the constant 
which is obtained from the experimental value 
of the frequency in the normal state of the 
diatomic CN molecule, 2069 cm, on the 
assumption of a simple harmonic oscillation, 
namely 16.3X10-> dynes/cm. The very small 
variation of fe; in the three halides leads to the 
conclusion that the carbon-nitrogen linkage is 
the same in all three compounds, while the steady 
diminution of fiz with increasing atomic weight 
of the halogen indicates a weakening in this 
linkage which is in accord with the chemical 
behavior of the compounds. F 

Our thanks are due to Mr. Edward Shriver of 
this department who very kindly obtained the 
microphotometer curves. 
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An improved apparatus for the study of Raman spectra 
has been developed. The use of filters including nitro- 
benzene for removing violet light, sodium nitrite for 
removing ultraviolet and copper nitrate for the removal of 
red light, has been discussed. The Raman spectra of 


n-butane, chlorobenzene, bromobenzene and iodobenzene 
have been reinvestigated and several new lines found in 
each. The Raman spectrum of alpha-chloropyridine has 
been studied for the first time and found very similar to 
that of chlorobenzene. 


HE problem of the structure of the benzene 
ring has long been of interest to chemists 
and has attracted the attention of many in- 
vestigators in the field of Raman spectra. In the 
five years since the opening of this field, at least 
thirty-five papers dealing with the Raman 
spectrum of benzene have appeared, and although 
the most careful and painstaking work has been 
done on this compound, authentic new lines are 
still being reported. In view of this fact and the 
importance of the problem of the structure of the 
benzene ring, it has seemed worth while to apply 
the best current technique to a repetition of the 
earlier work on the simpler aromatic derivatives. 
There seems to be a reasonable hope of obtaining, 
with the present methods, a series of spectra 
complete enough for the solving of the major 
points in the analysis of the modes of internal 
vibration of these molecules. 

With these facts in mind, we have photo- 
graphed the Raman spectra of a group of 
aromatic compounds and present, in this and a 
subsequent article, the results obtained. In 
addition to the data themselves, there are certain 
improvements in technique which may also be of 
interest; especially those designed to facilitate 
the study of compounds of high melting point 
and unusual fluorescence. 


I. APPARATUS 


In connection with this work, an illuminating 
apparatus of great power has been constructed. 


* From the dissertation submitted by John W. Murray 
to the Board of University Studies of The Johns Hopkins 
University in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy. 


(Figs. 1 and 2.) The illuminator is built around a 
Pyrex tube A of 48 mm inside diameter. The 
annular space between this tube and the Raman 
tube B, which is mounted at its axis, may be 
used for filters and for heating or cooling when 
necessary. The light is furnished by four quartz 
capillary mercury arcs G, of the type described 
by Langer and Meggers,' operating in running 
water. The water jacket can C, which also serves 
as a reflector, is of chromium plated brass and 
has the form of a combination of four ellipses as 
in the apparatus described by Anand.” This can 
is supported and held in shape by a set of hori- 
zontal plates D whose centers are cut out to fit 


4 


Fic. 1. Raman apparatus. 


1 Langer and Meggers, Bur. Standards J. Research 4, 711 
(1930). 
2B. Anand, J. Sci. Inst. 8, 258 (1931); 9, 324 (1932). 
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RAMAN SPECTRA OF RING COMPOUNDS. I 


Fic. 2. Diagram of reflectors. 


it and these in turn are bolted to two vertical 
plates E and the whole is mounted on a ring by 
means of the rods F. The arcs are aligned at the 
outer foci of the ellipses by adjustment of the 
racks HT and J which hold the bases of the arcs 
and the transite handles K in position. A large 
fraction of the light is thus reflected into the 
Raman tube which lies at the common focus of 
the set of ellipses. Scattered light passes down- 
ward through the window and the stop L and is 
reflected by the mirror M to the spectrograph 
which is placed at a distance of about 8 feet in 
order to eliminate reflected light from the walls 
of the tube. The spectrograph used in this work 
is a Bausch and Lomb glass prism instrument 
with F 10 lenses and has a dispersion of 17A per 
mm at 4000A and 50 at 5000A. 

For liquids, a tube of the Wood type may be 
used as shown in Fig. 1. Two types of tubes have 
been designed for the study of high-melting 
compounds which must be distilled continuously 
during the run on account of decomposition. 
(Fig. 3, I and II.) The first type is simpler but 
bulkier. The whole tube is of Pyrex glass and is 
heated with a winding of nichrome wire. The 
bulb N is heated by a separate circuit at a higher 
temperature so that distillation takes place from 
N to O. The bend below O is sooted to furnish a 
black background. In the second type, the two 
tubes are concentric and the depressed part 
serves as a condenser while the horn Q is sooted 
to give the black background. The distillate falls 
into the funnel at the top of the inner tube. When 
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Fic. 3. Raman tubes for continuous distillation (I and I1) 
and cell for use with three liquid filters (III). 


distillation is not required, the tube used is like 
type I with the distilling arm omitted. In either 
case, the tube is mounted inside of a straight 
jacket tube which replaces the Raman tube B in 
Fig. 1. The jacket may be lined with mica for 
thermal insulation. 

The filters are arranged in annular layers in the 
space between A and B, as in the apparatus 
described by Dadieu and Kohlrausch.* The cell 
shown in Fig. 3, III, permits the use of three 
layers of liquid filters, two in the cell and one 
outside. The cell may be simplified in an obvious 
manner when only two layers are needed. Glass 
filters may be used with this apparatus either as 
tubes or as sets of four plates. 

The use of filters is practically necessary for 
the location of the weaker lines of the spectrum. 
Three distinct purposes are served : simplification 
of the spectrum, removal of continuous back- 
ground from the mercury arc, and protection 
against decomposition and fluorescence from 
short wave radiation. For excitation with the 
blue mercury lines, a one mm layer of nitroben- 
zene will completely suppress the violet lines 
while transmitting most of the blue. It darkens 
with use but compares favorably, in this respect, 


3See Kohlrausch, Der Smekal-Raman-Effekt, Springer, 
Berlin (1931). 
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with quinine sulfate which is in general use for 
this purpose. A 2 mm layer of a saturated sodium 
nitrite solution will protect the nitrobenzene 
against ultraviolet light so that its transmission 
of the blue light is but slightly reduced by an 
eight hour exposure. Pfund‘ has recently recom- 
mended the use of sodium nitrite solution for the 
purpose for which nitrobenzene is used here. It 
has the advantage of stability but its relative 
transmission of blue to violet light is far inferior 
to that of nitrobenzene. It is, however, very 
useful for the purpose described above and also, 
in lesser concentrations, for removing ultraviolet 
when exciting with the violet lines. An excellent 
combination for excitation with the blue lines is: 
2 mm of saturated sodium nitrite solution, 1 mm 
of nitrobenzene, and 10 mm of saturated 
praseodymium nitrate solution to remove the 
continuous background. The nitrobenzene may 
be diluted with alcohol and crystal violet may be 
added to remove the blue-green region. For 
excitation with the violet lines, a 2 mm layer of 
Corning Red Purple Ultra glass is useful for 
removing the blue. This compares favorably with 
iodine in carbon tetrachloride which is often used 
for this purpose. A layer of sodium nitrite is used 
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with one of these filters and is more efficient than 
Noviol glass. When excitation with the green line 
is necessary, as with some colored or fluorescent 
substances, a combination of potassium chro- 
mate, copper nitrate, and neodymium (or 
didymium) nitrate gives almost complete isola- 
tion of the 5461 mercury line. Excitation by this 
line is sometimes of value in deciding questions as 
to the origin of lines of low displacement. 
. The apparatus described above will record the 
strong Raman lines of benzene in thirty seconds. 
Exposures of several hours with strong filters are 
generally used in order to obtain the weaker 
lines. 
II. RESULTS 

The Raman spectrum of n-butane was in- 
vestigated for the purpose of obtaining data for 
thermal calculations. Since this work was done, 
results of an earlier investigation by Bhaga- 
vantam® have appeared which are in good 
agreement with those reported here. The material 
used was 99 percent commercial n-butane and 
was distilled twice before using. Chlorobenzene, 
bromobenzene and iodobenzene have been re- 
examined and a number of new lines have been 
found. These compounds have been studied by 
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Fic. 4. Raman displacements of compounds studied. 


4A. H. Pfund, Phys. Rev. 42, 581 (1932), 


5S. Bhagavantam, Ind. J. Phys. 6, 595 (1931). 
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several earlier workers. The material used was 
specially purified for specific heat studies by Mr. 
D. E. Teets of this laboratory. Alpha-chloro- 
pyridine has been studied for the first time. This 
substance was prepared by Dr. Lyman C. Craig 
of this laboratory. The results obtained are 
recorded in the following tables and are presented 
graphically in Fig. 4. 

The procedure advocated by R. W. Wood,’ of 


taking two plates of each substance using 


RING COMPOUNDS. I 409 
different exciting lines, has been employed. The 
violet and blue mercury lines have been used in 
all cases and the green line in some. In the cases 
in which it has not been necessary to remove the 
ultraviolet to prevent decomposition, Raman 
lines excited by the 3650 group appear with the 
lines from the violet. This is very evident in the 
photograph of n-butane shown in Fig. 5. The use 
of several exciting lines makes it possible to 
decide questions arising from excitation by the 


i T 


Fic. 5. Raman spectra. 


. n-Butane excited by \=4358 
. n-Butane excited by \=4047 
. a-Chloropyridine excited by \=4358 
. a-Chloropyridine excited by ’=4047 
Chlorobenzene excited by \=4358 
. Chlorobenzene excited by \=4047 


* Pringsheim and Rosen, Zeits. f. Physik 50, 741 (1928); 
Petrikaln and Hochberg, Zeits. f. physik. Chemie B3, 217, 
405 (1929); Dadieu and Kohlrausch, Monats. f. Chem. 52, 
379 (1929); 57, 488 (1931); Dadieu, Kohlrausch and 
Pongratz; Monats. f. Chem. 61, 426 (1932); Langseth, 
Nature 124, 92 (1929); Fujioka, Sci. Pap. Inst. Phys. 


7. Bromobenzene excited by \=4358 
8. Bromobenzene excited by \ =4047 
9. Iodobenzene excited by \=4358 
10. Iodobenzene excited by \ =4047 
11. Iodobenzene excited by \=5461 


Chem. Res. Tokyo 11, 205 (1929); Cleeton and Dufford, 
Phys. Rev. 37, 362 (1931); Bhagavantam and Venkate- 
swaran, Proc. Roy. Soc. (London) A127, 360 (1930); Wood, 
Phys. Rev. 38, 2168 (1931). 

7 Wood, Phys. Rev. 38, 2168 (1931). 
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weaker companions of the 4358 and 4046 
mercury lines. The fact that the violet mercury 
lines themselves are not visible in the plates from 
which the data for the blue excitation are taken 
removes the possibility of error from this source 
in assigning origins to the Raman lines. This is of 
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great importance in the region of 12-1400 cm” 
displacement. Fig. 5 shows the photographs 
obtained. 

The values recorded are averages from several 
plates in most cases. In Tables I-V, the wave 
numbers of observed lines are given in the first 


TABLE I. n-butane. 


v 


9 


bey 


Ap 


22557 
22510 
22150 
22105 
22045 
22002 
21987 
21967 
21958 
21880 
21868 
21843 
21826 
21794 
21789 . 
21766 


24479 
24457 
24422 


24356 


24282 
23916 
23877 
23748 


23729 


23689 
23649 
23557 
23252 
22600 


or OFN 


438 
428 
788 
833 
2660 
2703 
951 
2738 
980 
1058 
(1070) 
2862 
2879 
2911 
1149 
2939 


> 


OK 


(338) (0), 428 (3), 787 (2), 832 (6), 951 (1), 980 (2), 1058 (3), oy (0), 1149 (1), 1293 (16), 1450 (36), 2670 (0), 
2702 (0), 2739 (1), 2864 (5), 2879 (10), 2912 (3), 2941 (8), 2963 (3). 


TABLE II. Chlorobenzene. 


v 


my 
9 


24307 
24287 


24091 


24001 
23964 
23962 
23938 
23813 
23702 
23680 
23618 
23581 
23549 
23530 
23512 
23491 
23430 


23352 


23336 
23260 


aae 


AV‘ eae 


21780 
21763 
21700 
198 21678 
275 21637 
299 21614 
331 21572 
415 21564 
418 21539 
465 21536 
615 21500 
614 
701 21459 
702 21447 
742 21410 
830 21377 
1004 21354 
1002 20609 
1022 19869 
1003 18110 
1023 18005 
1084 17907 
1121 17697 
17317 


198 (7), 275 (0), 299 (2), 331 (0), 419 (6), 465 (1), 614 (3), 703 (5), 743 (1), 832 (2), 1004 (10), 1023 (7), 1084 (5), 
1122 (2), 1159 (3), 1174 (2), 1324 (0), 1366 (0), (1402) (0), 1438 (1), 1475 (1), 1561 (1), 1584 (6), (2329) (0), 


3005 (0), 3027 (0), 3069 (10), (3141) (0), 3166 (0). 


4 AG 
1293 
2876 
2907 
2936 
1450 
2670 
2702 
2739 
2864 
2879 
2912 
2941 
2963 
| 
1175 
{ 3005 
3027 
3068 
1324 
1366 
(3141) 
“a 3166 
(1402) 
i 1438 
1479\ 
1580) 
3069 
1585 
1561 
1584 
(2329) 
303 
611 
991 
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TABLE III. Bromobenzene. 


I 


24388 
24286 


24229 


24199 
24089 
24032 
23869 
23840 
23701 
23681 
23631 
23542 
23527 
23512 
23491 
23443 
23253 
23248 
23124 
23113 
22753 
22719 


22681 


oo 


22622 
22585 


22323 


22265 
22200 
22143 
22107 
22035 
22018 
21994 
21972 
21936 
21915 
21866 
21834 
21778 
21761 
21701 
21664 
21640 
21637 
21616 
21553 
21502 


SK SSS 


oo 


185 (8b), 260 (16), 317 (10), 353 (1b), 616 (3), 673 (5), 738 (1), 795 (0), 831 (1), 1002 (10), 1023 (5), 1072 (4), 
1160 (3), 1177 (2), (1274) (0), 1301 (0), 1322 (0), 1442 (1), 1477 (1), 1580 (55), 1614 (0), 1667 (0), 2293 (0), 
2333 (0), 3004 (0), 3063 (10), 3153 (1), 3203 (0). 


TABLE IV. Jodobenzene. 


v 


v 


~ 


24437 
24247 
24090 
24048 
23966 
23904 
23859 
23704 
23685 
23640 
23586 
23544 
23514 
23497 
23454 
23353 
23329 


23263 


23231 
23200 
23130 
23103 
22770 
22720 


FS Araaas’s AAAS 


22671 
22623 
22550 
22383 


22324 


22282 
22210 
22100 
22037 
22028 
21995 
21976 
21938 
21919 
21878 
21845 
21780 
21759 
21702 
21681 
21647 
21619 
21579 
21564 


> 


MOR 


21500 
21469 


21458 
21425 
21365 
21323 
21278 
20612 
19943 
19878 
18565 
18146 
18080 
18039 
17693 
17652 
17316 
17295 
17154 
16735 
15248 
15175 


> 


168 (5), 228 (2), 268 (10), (319) (0), (388) (0), 614 (2), 656 (5), 728 (1), 838 (1), 1001 (10), 1019 (6), 1062 (2), 
1158 (3), 1287 roy i319 46), 339 (0), 449 (0), 1$73 (5), (0), (506 (0), (2326) (0), 


(3000) (0), 3059 (10), 3141 (1). 


a 
Ay Exc. I > Exc. I Ap 
| 317 10 316 || 21496 1442 
al 3067 |_| 1b 1477\ 
ve 615 1 
3008} 3 on} 21453 3063 
317 673 21415 1580 
616 738 || 21362 3154 
673 795 21358 1580 
(836) 831 21324 1614 , 
676 1004 21271 1667 a © 
- 1004 1021 20645 2293 4 
1024 1001 20605 2333 - 
1074 1023 19876 3062 
1163 1002 || 18974 — 666 Le i 
1178 1023 18559 (—251) — 
1004 1072 18123 185 
1025 1161 18044 264 
1073 1160 17989 319 ‘Sa 
1452 1177 17315 993 ee 
—310 3004 17292 1016 
1581 (1274) || 17250 1058 ee 
185 1301 
2 1322 tin 
3152 
314 3203 i. 
i Ex. Ap Exc. I Av > Exc. Ap 
268 10 267 1438 
269 0 315 
615 0b (388) 1570 i 
657 0 612 3058 
9 614 15 ae 
a2 | 2v 1573 
657 656 1615 
1001 728 1660 
1020 838 (2326) 
1065 1002 (2995) 
—648 (1011) 1 3060 
1161 1000 —257 
1002 1019 162 
1019 1000 228 = 
1062 1019 269 De 
1163 1060 615 a 
1187 1150 656 bo ‘ 
1474 3003 1154 — 
—262 1257 1573 ae 
1575 3058 3060 Be 
—165 1319 3163 
5b 168 1359 
3bv 218 3141 
: 
Mi 
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TABLE V. Alpha-chloropyridine. 


Exc I Av ? Exc I Ap v Exc. q Ap 
24388 d 0 (317) 23246 d 0 (1459) 21890 h 8 1048 
24272 d 1 433 23130 d 2b 1575 21852 h 3 1086 
24084 d 1 621 23113 h 1 —175 21817 h 3 1121 
23976 d 2 729 22744 h 8 194 21784 h 2 1154 
23796 e 0 720 22621 h 1 517 21698 h 1 1240 
23708 d 4 997 22564 g 1 431 21650 h 2 1288 
23655 d 3 1050 22509 h 8 429 21648 d 10 3057 
23615 d 1b 1090 22318 h 4 620 21625 d 6 3080 
23582 d 1b 1123 22270 g 0b 725 21577 h 0 (1361) 
23550 d 1 1155 22210 h 3 728 21573 d 1b 3132 
23519 e 1 997 22160 h 1 778 21544 d 0 (3161) 
23466 € 1 1050 22052 d 1 2653 21512 h 0 (1426) 
23414 d 0b 1291 22046 Pd 1b 993 21478 h 0 (1460) 
23268 h 0 (—330) 22002 g 1 993 21463 e 1 3053 
21981 d 1 2724 21439 e 1 3077 
21967 h 1 971 21366 h 5 1572 
21945 h 10 993 19851 h 2b 3087 


194 (8), 317 (1), 429 (8), 620 (4), 728 (3), 778 (1), 971 (1), 993 (10), 1048 (8), 1086 (3), 1121 (3), 1154 (2), 


1240 (1), 1288 (2), (1361) (0), (1426) (0), (1460) (0), 1573 (5), 2653 (1), 2724 (1), 3057 (10), 3080 (6), 3132 (1), 


(3161) (0). 


column, symbols for the exciting line in the 
second, the intensity and character of the line in 
the third, and the displacement assigned to it in 
the fourth. The symbols for the exciting lines are 
taken from Kohlrausch.* 


a27388cm™ 424705 f 23039 k 18308 
b 27353 e 24516 = g 22995 
c 27293 h 22938 


The intensities given are merely visual estimates 
of the relative blackness of the lines and are not 
strictly comparable from one substance to an- 
other. The letter 5 signifies broad, v shaded to the 
violet, and r shaded to the red. The values 
appearing in parentheses in the body of the tables 
are uncertain due to extreme weakness or 
interference of background. Those in parentheses 
in the summaries are uncertain from the above 
causes or on account of doubt as to their origin. 

The Raman frequencies are summarized at the 
end of each table and the intensities are given in 
parentheses after the displacements. The values 
are probably accurate to within 3 cm~ or better 
for the stronger lines but the error for the very 
weak ones may be as large as 10 cm. The plates 
were measured on a Gaertner comparator and the 
wave numbers determined from a large dis- 
persion curve constructed from an iron com- 
parison spectrum. The accuracy of the lines 
excited by the green line is less on account of the 
low dispersion in that region. 


III. Discussion 


The results reported here for n-butane confirm 
those of Bhagavantam except that the line he 
reported at 960 has been resolved into two lines 
at 951 and 980. The lines at 320 and 1067 were 
not well substantiated in this work though 
traces of them were found. 

Dadieu and Kohlrausch® reported a line for 
chlorobenzene at 484 in their first paper on this 
substance but did not report it in a later paper in 
which the work was repeated.® A line has been 
found at 465 which might be the one in question 
considering the error in measurement of such 
weak lines. The line at 3133 reported in the latter 
paper was not substantiated in this work though 
a trace of it was found. All other previously 
reported lines have been confirmed. New lines 
have been found at: 275, 331, 1174, 1324, (1402), 
1438, 1475, 1561, (2329), 3005, 3027. The 
justification for the assignment of a line to 1475 
is not apparent from the table and is based on 
the fact that the observed line at 21459 cm" is 
stronger than the line at 21410 which would not 
be the case if the former were due to the 1584 line 
excited by the 4339 line alone. This same 
reasoning applies to the other compounds in 
which this line is found. 


8 Dadieu and Kohlrausch, Monats. f. Chem. 52, 379 
(1929). 

* Dadieu, Kohlrausch and Pongratz, Monats. f. Chem. 
61, 426 (1932). 
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In the case of bromobenzene, Fujioka® reports 
a line at 222 which is derived from the observed 
line at 22716 and is assumed to have been 
excited by the 4358(h) mercury line. A similar 
line was found in the present work but it was 
deemed more probable that it was excited by the 
4348(g) mercury line and is so reported in this 
paper. The line at 1370 reported by Dadieu and 
Kohlrausch was not found though several other 
weak lines were found in this region. The reso- 
lution of the strong line at 3063 into two lines 
reported by Fujioka was not satisfactorily con- 
firmed though a trace of it was found in one 
plate. In other respects, the previously reported 
lines have been confirmed. New lines have been 
found at: 353, 738, 795, 831, 1274, 1301, 1322, 
1442, 1477, 1614, 1667, 2293, 2333, 3004, and 
3203. 

In the case of iodobenzene, a line at 462 was 
reported by Dadieu and Kohlrausch® in the 
earlier paper but was not reported by them in a 
later paper.’° This line was not found in the 
present work. Ali other previously reported lines 
have been found. New lines have been found at: 
228, 319, 388, 728, 838, 1179, 1257, 1319, 1359, 
1440, 1472, 1615, 1660, (2326), (3000) and 3141. 
The line at 228 is not distinguishable from the 
268 line excited by 4348(g) when excited by the 
4358(h) line and is masked by the mercury 
radiation when excited by the violet line. For 
this reason, its identity has not been established 
before though the line has been observed. This 
question has been settled by excitation with the 
green line. The incompleteness of the previous 
work on this substance is largely due to the 
decomposition caused by ultraviolet light. This 
difficulty has been overcome by protecting it 
with sodium nitrite. 

The marked similarity between the mono- 
halogen derivatives of benzene is very evident 
from Figs. 4 and 5. Nearly all of the lines which 
do not correspond throughout the series are so 


' Dadieu and Kohlrausch, Monats. f. Chem. 57, 488 
1931). 
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weak that they may well have been missed in the 
cases where they are absent. The line at 1122 in 
chlorobenzene appears to be missing in the 
bromo and iodo compounds; but this is not 
certain for it may coincide with the position of 
the 1159 line excited by 4348(g). In this case, 
however, the intensity must have fallen off 
sharply. The weak companion of the 1584 line in 
chlorobenzene is not found in the other com- 
pounds but two very weak lines appear on the 
other side of it at 1615 and 1660. 

The lines at 198, 299, 419, 465 and 703 in 
chlorobenzene move toward the exciting line by 
very appreciable amounts in passing to the 
bromo and iodo compounds and thus appear to 
correspond to motions in which the carbon to 
halogen bond plays an important part. The lines 
at 614, 832, 1004, 1023, 1159, 1174, 1438 and 1475 
remain nearly constant and hence pertain to 
motions in which this bond is not very much 
involved. The lines at 743, 1084, 1584, 3069 and 
3166 are very slightly lowered on increasing the 
mass of the substituent and thus appear to 
involve this bond to a slight extent. The re- 
mainder of the lines are so weak that too much 
confidence should not be placed in apparent 
shifts though the lines in the region of 1200 to 
1400 appear to have a genuine shift to lower 
values. It is of interest in connection with the 
lines at about 2300 that a line at 2291 has 
recently been found for benzene by Wood and 
Collins." 

Many of the new lines found for the benzene 
derivatives are of particular interest in that they 
were predicted by Mr. D. E. Teets of this 
laboratory from studies of the vibrating me- 
chanical models which will be published in the 
near future. 

The obvious similarity of the spectra of 
chlorobenzene and alpha-chloropyridine would 
appear to indicate that the ring structures of the 
two compounds are of the same type. 


1! Wood and Collins, Phys. Rev. 42, 382 (1932). 
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The infrared spectrum of ozone is discussed. Gerhard’s 


three fundamental frequencies, 528 cm~, 1033 cm™, and 
1355 cm~! are shown to be in somewhat imperfect agree- 
ment with experiment. They are nevertheless sufficiently 
reliable to justify a spectroscopic calculation of the free 
energy of ozone, in view of the unreliability of the experi- 
mental data. Ozone is found to be much less stable than 
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experiment had shown; Sxs3..=57.0 cal./deg./mole, as 
compared with S2:;;=86.5 based on e.m.f. measurement. 
The maximum partial pressure of O; in an equilibrium 
mixture of O, Oz, and O; at a total pressure of 1 atm. is 
1.0 10~* atm. at 3400°K. The 1 percent ozone obtained by 
heating filaments to 2300°K in liquid oxygen must have 
been formed in the cold regions from atomic oxygen. 


BOUT thirty years ago, strenuous efforts 
were made to determine the free energy of 
ozone. The various e.m.f. measurements which 
were attempted met with the general difficulties 
that the electrode reactions were not known and 
not reversible.* Many attempts to detect ozone 
in very rapidly cooled oxygen failed. When it was 
realized that the high rate of decomposition of 
ozone even at 200° doomed such experiments in 
advance, Fischer and Braehmer‘ burnt hydrogen, 
carbon monoxide, and other gases under liquid air 
or oxygen, and also heated Nernst filaments in 
liquid oxygen. They obtained in this way solutions 
of ozone in concentration up to 1 percent; this 
represents an upper limit to the equilibrium con- 
centration at the filament temperature (2300°K) 
since there may have been a considerable frac- 
tionation effect as the oxygen evaporated. It will 
be seen from what follows that the ozone found 
in this work must have been formed almost 
entirely from atomic oxygen, and that no exper- 
iment of this sort can possibly give the equilib- 
rium concentration of ozone. Recently, Wulf® has 
decided that an empirically estimated entropy 
value is more reliable than the very high “om 
given by the foregoing experiments. 

The existing spectroscopic data for ozone are 


! Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

* Associate physical chemist, U. S. Bureau of Mines, 
Pittsburgh Experiment Station, Pittsburgh, Pennsylvania. 

3 Jahn, Zeits. f. anorg. Chem. 60, 332 (1908). 

4 Fischer and Braehmer, Ber. 39, 940 (1906). 

5 Wulf, J. Am. Chem. Soc. 54, 156 (1932). 
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rather scanty, but they seem to be sufficient to 
yield entropy and free energy values far superior 
to those available from any other source. The 
early work on its infrared absorption’ showed 
bands at 1.1, 3.7, 4.8, 6.6, 7.6, 9.9, and 11.35y. 
Recently, Gerhard’ has examined the 4.7, 7.39, 
9.6, and 11.384 bands with a high-dispersion 
grating spectrometer; the resolution was appar- 
ently just insufficient to yield individual lines. 
The 4.74 (2108 cm~') and 7.394 (1355 
bands show a strong zero branch. The band at 
9.64 is apparently complex, with a zero branch 
band at 9.474 (1055 cm~') superimposed on a. 
doublet band at 9.67 (1033 cm~'). The weak 
band at 11.38 shows a single peak. Since it does 
not resemble the others, and does not fit the 
frequency scheme, it is tentatively assigned to 
nitrogen pentoxide by Gerhard; the method of 
preparation used is such as to make this impurity 
relatively unlikely. Gerhard adopts the three 
fundamental frequencies 528 cm, 1033 cm™ 
and 1355 cm~, as ve, v3 and 1, respectively. The 
experimental conditions did not permit study in 
the region of v2. The appearance of the bands 
suggests rather definitely that there are three 
different moments of inertia. It is therefore 
reasonable to consider the molecule an isosceles 
triangle. For this case the three normal modes of 
vibration have different frequencies. Analysis 


6 Ladenburg and Lehman, Ann. d. Physik 21, 305 (1906); 
Warburg and Leithauser, ibid. 23, 209 (1907); 28, 313 
(1909). 

7 Gerhard, Phys. Rev. 42, 622 (1932). 
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shows’ that in two of them (v; and v2, according 
to the convention adopted) the electric moment 
vibrates along the Y axis (the symmetry axis) 
and for v3 perpendicular to this axis. Now Den- 
nison® has shown that when the electric moment 
moves along the least moment of inertia, the 
bands show a zero branch; when it moves along 
the middle moment of inertia, a doublet struc- 
ture. Thus, if the three fundamental frequencies 
have been correctly selected, the assignment 
made by Gerhard supposes that the 528 cm 
band will show a zero branch, and that the 
triangle has an apex angle less than 60°. If 528 
cm~! has a doublet structure we would have to 
take v;=1033 ve=528 and v3=1355 
cm~'. The middle moment of inertia would lie 
along the figure axis, and the apex angle would 
be greater than 60°. The 1055 cm~ band, which 
is assigned to 2v2, should then show a doublet 
structure, instead of a zero branch. We cannot 
consider 1055 cm! a fundamental, and 1033 
cm-! the overtone of 528 cm~'; for if we did, 
528 cm by elimination would become a doublet, 
and »;. But, in view of the alternation shown by 
Dennison to exist for v3, 2vs would be a zero 
branch, in contrast to the observed structure. 
Thus, if the three fundamental frequencies have 
been correctly selected, there seems to be no 
alternative to Gerhard’s assignment of them, and 
his conclusion that the apex angle is less than 
60°. There is certainly no obvious way in which 
different fundamentals can be selected. 

The comparison of these frequencies with the 
theoretical ones for a simple model, however, is 
not very satisfactory. We take for the potential 
energy 


(1/2) (Kyi? 


where x;; is the relative displacement of atoms 
i and j; the apex atom is numbered 3. It can 
then be shown that of the three possible assign- 
ments of the observed fundamentals to the 
different modes of vibration, only the one made 
by Gerhard gives a real solution. It gives two 
solutions: for one the apex angle is 59° 16’, 
K/K’=0.2198, and the ratio of least to middle 
moment of inertia is 0.9706; for the other, the 
apex angle is 117° 16’, K/K’=1.662, and the 


* Dennison, Phil. Mag. 1, 195 (1926). 
* Dennison, Rev. Mod. Phys. 3, 280 (1931). 
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ratio of least to middle moment of inertia is 
0.1245. Reference to the diagrams given by 
Dennison shows that for the former case the 
separation into zero and doublet bands would 
not have developed to a noticeable extent, while 
for the latter case the doublet bands would have 
a very coarse structure; neither case corresponds 
well with the observed structures. Furthermore, 
the ratio K/K’ in both cases has an unexpected 
value.!° 

Jakowlewa and Kondratjew'! have measured 
some of the ultraviolet bands of ozone. They 
propose frequencies different from those of 
Gerhard, and favor a linear form for the molecule. 
Consideration, however, shows that their data 
are not well suited to analysis. The bands are not 
located with a certainty of better than 5 cm™, 
and there are so many bands that the chance 
coincidences outnumber the real ones several fold. 

The following calculations have been based on 
Gerhard’s frequencies. The symmetry number is 
taken as two, in agreement with his conclusion. 
The value for the moment of inertia has been 
computed, assuming O-O distances of 1.2, 1.4 
and 1.4X10-* cm. It is unlikely that the error 
in F because of a possible incorrect assignment 
of the frequencies exceeds 500 cal. at 1000°K, 
2500 cal. at 2000°K or 5000 cal. at 5000°K. 
We assume harmonic oscillators; the resulting 
error in F is perhaps 1000 cal. at 5000°K. The 
only other value for the symmetry number which 
seems at all possible is three. The product of the 
three moments of inertia will scarcely be wrong 
by more than 25 percent, unless the apex angle 
is considerably greater than 120°, and the ratio 
of least to middle moment of inertia less than 
0.1; as has been pointed out, the appearance of 
the bands is against this. We assume that the 


1 Since this paper was written, Badger and Bonner, 
Phys. Rev. 43, 305 (1933), have criticized Gerhard’s 
frequency assignment on grounds very similar to those 
mentioned above. On this basis, and with reference to 
Wulf’s analysis of certain electronic bands, Proc. Nat. 
Acad. Sci. 16, 507 (1930), they favor »; = 1050, ».» =440 and 
v3=1355 cm~!. These frequencies do give one reasonable 
solution, an apex angle of 81° 0’ with K/K’=0.619. The 
agreement with the infrared bands is definitely less good 
than that of Gerhard’s assignment, but the existing data 
scarcely justify a definite decision. 

11 Jakowlewa and Kondratjew, Phys. Zeits. d. Sow. 1, 
471 (1932). 
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TABLE I. Free energies, equilibrium constants and equilibrium concentrations for ozone. 


Logio K 


Equilibrium partial pressure at 


/T (Oz)3 (Oz)(O) 1 atm. total pressure 
(Os) (Os) Oz O O; 
250 47.121 — — — 
298.1* 48.558 28.490 — 11.899 1.000 4.0 3.21079 
300 48.828 28.285 —11.829 1.000 7.8107 5.2 x 10-29 
350 49.972 24.793 — 9.150 1.000 1.110-* 1.61075 
400 51.142 22.149 — 7.159 1.000 5.0 107° 7.11073 
450 52.200 20.092 — 5.606 1.000 2.0 10-26 8.11072! 
500 53.167 18.447 — 4.341 1.000 1.61073 3.6 
550 54.061 17.101 — 3.338 1.000 3.6107! 7.9X 1078 
600* 54.890 15.980 — 2.484 1.000 3.41079 1.0107 
700 56.408 14.214 — 1.144 1.000 4.410716 6.1 X10-" 
800 57.759 12.891 — 0.128 1.000 9.6107" 1.3107" 
900 58.992 11.858 + 0.659 1.000 6.4107 1.4107" 
1000 60.116 11.031 + 1.291 1.000 1.81071 9.3 X10-" 
1250 62.573 9.541 + 2.431 1.000 7.8X 1078 2.9107" 
1500 64.689 8.537 + 3.186 1.000 4.41076 2.9x10~ 
1750 66.533 7.816 + 3.724 1.000 7.8X10- 1.51078 
2000 68.073 7.291 + 4.148 0.999 7.21074 31x19 
2500 70.812 6.536 + 4.722 0.985 0.015 2.91077 
3000 73.096 6.032 + 5.108 0.888 0.112 7.8X107 
4000 76.777 5.400 + 5.595 0.237 0.763 4.6107 
5000 79.687 5.022 + 5.894 0.017 0.983 2.2 X 10-8 


* The equilibrium constants for O:=2 O and the percentage dissociations given by Johnston and Walker at 298.1° and 
600° are incorrect; their (F°—Eo°)/T values, however, are apparently correct at these temperatures. The method of 
graphical interpolation they have used to obtain (F°—E,°)/T at intermediate temperatures is apparently less accurate 
than the procedure of Giauque and Overstreet, J. Am. Chem. Soc. 54, 1731 (1932), but is good enough for the purposes of 


the present paper. 


normal state is '> and that there are no other 
electron levels low enough to require consider- 
ation; this may be incorrect. Some of these 
errors would increase the equilibrium concen- 
tration of ozone, others decrease it. The calcu- 
lated values are probably correct within a factor 
of two, and certainly within one of five. They 
may be much better than this, but at the worst 
are a great deal better than any experimental 
results. 

The heat of formation of one mole of ozone 
at constant volume is given as —34,220+180 
cal., by the recent work of Gunther, Wassmuth 
and Schryver.” From this we calculate that 
AE,® = 34,513 cal. 

The methods used in the calculations have 
been given in a previous article.'* The entropy of 
ozone at 25°C is found to be 57.0 cal./mole/deg. 
This is somewhat higher than the estimate of 
52.6 made by Wulf, but much lower than the old 
value of S273 = 86.5 based on e.m.f. measurements. 


12 Gunther, Wassmuth and Schryver, Zeits. f. physik. 
Chemie A158, 297 (1932). 
'3 Kassel, J. Am. Chem. Soc. 55, 1351 (1933). 


Table I gives the (F°—E,°)/T values for 
ozone, the equilibrium constants for O3 = (3/2)O2 
and for O;=O2+O, and the partial pressures of 
O, O2 and O; at equilibrium at a total pressure 
of 1 atm. These latter values are also shown in 
Fig. 1. The (F°—E,°)/T values for O and O: 
are taken from Johnston and Walker.'* The 
values for Oz do not take account of the low- 
lying !A state predicted by Mulliken’® and 
Hickel,'® and confirmed by specific heat meas- 
urements.!”7 The inclusion of this level would 
make O2 more stable with respect to all reac- 
tions. An approximate correction has been made 
for this level, and the new equilibrium values 
have been calculated. The curves for O and O2 
are not shifted enough to show in Fig. 1; the 
ozone curve is lowered somewhat, but the effect 
is not large. In view of what has been said 


14 Johnston and Walker, J. Am. Chem. Soc. 55, 172, 187 
(1933). 
% Mulliken, Phys. Rev. 32, 186 (1928). 
16 Hiickel, Zeits. f. Physik 60, 442 (1930). 
17 Lewis and von Elbe, Phys. Rev. 41, 678 (1932). 
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Fic. 1. Equilibrium partial pressures of O, O2 and O; at 1 atm. total pressure. 
The upper curve for ozone is based on the (F°—E,°)/T values for oxygen 
given by Johnston and Walker; the lower curve includes an approximate correc- 
tion for the !A state of oxygen; both ozone curves should be raised by an amount 
ranging from 7.6 percent at 2000° to 3.1 percent at 5000°, to correspond with the 
values of Table 1; such a shift is considerably less than the probable error of the 


values. 


previously with regard to possible errors, it is 
clear that there is no temperature at which the 
ozone partial pressure is as much as 5X10-* 
atm., at a total pressure of 1 atm. The ozone 
found by Fischer and Braehmer must have been 
formed from atomic oxygen, for which the 
equilibrium concentration at 2300°K is 0.0052, 
comparing favorably with their estimate of 1 
percent ozone. The calculations also show that 
there is no ozone contribution to the apparent 
specific heat of oxygen at high temperatures, and 
thus make it more certain that Lewis and von 
Elbe were justified in associating the observed 


increase of specific heat with excitation to the 
1A level of Oc. 

It has not seemed worth while to tabulate the 
other thermodynamic functions; quite good 
values for H and S may be obtained by graphical 
or analytical differentiation; a very satisfactory 
procedure is to fit F/T at three successive tem- 
peratures to F/T=A+BlnT+C/T. This same 
equation will give specific heats good to about 
0.5 cal./deg. The only reason for giving F/T 
and log K to three decimal places is to permit 
such calculations as these, which could not be 
made if fewer were given. 
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The Mechanism of Substitution Reactions 


A. R. Otson, Department of Chemistry, University of California 
(Received April 15, 1933) 


For many substitution reactions of the type AB+C 
—AC+B, the strength of the bond between A and B is 
too great to account for the rate of the reaction. For such 
reactions addition or simultaneous addition and dissoci- 
ations must occur. From the quantum mechanical con- 
ception of bonds, it is shown that a spatial inversion of A 


occurs for every such substitution. This theory is applied in 
detail to some experiments on the Walden inversion, and a 
definite correlation between rotational and configurational 
changes is obtained. Its use is studying molecular rearrange- 
ments and cis-trans isomerization is sketched. 


HE assumption that all substitution re- 
actions, AB+C-AC+B, can be divided 
into two mutually exclusive classes, i.e., (1) those 
in which dissociation of AB is the primary step, 
and (2) those in which the formation of the 
complex ABC is the primary step, is predicated 
on the postulate that each type of atom or 
molecule can be characterized by a definite size. 
However, dissociation and association are not 
instantaneous events, but are processes which 
require finite time, and therefore we must regard 
the above classes as limiting cases of a more 
inclusive group in which dissociation and 
association may be simultaneous phenomena. 
The transition from the ordinary addition 
product to the simultaneous association and 
addition can be illustrated by the gaseous 
reactions 


Cl, (1) 
(2) 
H.+ (3) 


For all of these reactions, the heats of dissociation 
of the diatomic molecules is so large that in each 
case we can neglect primary dissociation. 
Rollefson and Eyring have applied the theory 
developed by London? to reaction (1). They 
found that 4.5 kilocalories activation energy was 
required to form the addition product, and 8.6 
kilocalories was required to dissociate it. Further- 


1G. K. Rollefson and H. Eyring, J. Am. Chem. Soc. 54, 
170 (1932). 
2 F. London, Zeits. f. Elektrochemie 35, 552 (1929). 


more the separation of the chlorine atoms in the 
original molecule increased only very slightly as 
the third atom approached. We would therefore 
be justified in assuming that in this case addition 
preceded dissociation. 

For the hydrogen reaction Farkas* found that 
atomic hydrogen was a catalyst for the con- 
version of para hydrogen into ortho hydrogen. 
London’s theory was applied to this reaction by 
Eyring and Polanyi.‘ They found that it re- 
quired an activation energy of 11 kilocalories. 
They showed also that as the hydrogen atom 
approaches the molecule along the molecular 
axis, the atoms composing the original molecule 
begin to separate. The molecule H; needs about 
1.5 kilocalories to decompose it, and so here 
again an addition compound could exist, but toa 
much more limited extent than in the chlorine 
case. 

For the third reaction, the formation of an 
addition compound in the ordinary sense is 
impossible for the potential energy diagram of 
the reaction shows no region of stability for 
H.Br. We are, therefore, dealing with simul- 
taneous dissociation and addition. This evanes- 
cent addition product has this, however, in 
common with its more stable prototypes: The 
system requires much less activation energy if 
the bromine atom approaches the hydrogen 
molecule along the molecular axis than if it 
approaches it from any other direction. 


3L. Farkas, Zeits. f. physik. Chemie B10, 419 (1930). 
‘H. Eyring and M. Polanyi, Zeits. f. physik. Chemie 
B12, 279 (1931). 


418 


q 
2 
3 
js 
% 


d in 
ida 
ynal 
ige- 


mie 


MECHANISM OF SUBSTITUTION REACTIONS 419 


In the simple systems which we have con- 
sidered, experimental proof of this statement 
cannot be obtained, but if the theory is extended 
to reactions such as 


AB+C-AC+B- 
AB+C>AC+B 
AB+C+—AC+Bt 


where A in particular may be a complex radical, 
we will be forced to certain conclusions, some of 
which can be subjected to experimental test. 
The strict application of London’s theory to 
such complex systems is out of the question. 


Not only are important data lacking, but the 
calculations involved in a quantitative solution 
of the problems are, at the present time, hope- 
lessly burdensome. We can, however, get a 
qualitative picture from the type of potential 
energy diagram that would be expected for 
certain changes in the system. Thus let AB be 
the molecule R,R2R;CCl in which the bonds 
between the R’s and C are strong enough so that 
Cl can be identified with B. By plotting the 
energy of the system, e, against the distance 
separating R:R2R;C from Cl (or between the 
corresponding ions) we obtain a curve of the type 
shown in Fig. 1, where 7 is the equilibrium 
separation 


e— 


Fic. 1. Type of potential energy curve for the dissociation of R:R2RsCX into R;R2RsC and X, 
or into the corresponding ions. 


According to Pauling’ the average value of 
Dy for the carbon chlorine bond in substituted 
hydrocarbons is about 80 kilocalories per mol. 
This refers, however, to dissociation into a 


*L. Pauling, J. Am. Chem. Soc. 54, 3570 (1932). 


neutral atom and a radical. For dissociation 
into ions, D would be almost twice as much due 
to the difference between the ionizing potential 
of the radical and the electron affinity of 
chlorine. If the substance is in solution, cor- 
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rections for hydration and the dielectric constant 
of the solvent would have to be applied. This 
latter correction may be much smaller than is 
ordinarily supposed, for most of the work of 
separation is done at very small distance and 
here the effective dielectric constant seems to 
approach unity, as is shown by the fact that the 
frequencies of Raman lines are almost inde- 
pendent of the nature of the solvent. Fortunately 
for those substances which are to be discussed, 
we were able to obtain experimental® proof of the 
assumption of a large value of the heat of 
dissociation into ions, not only for some chlorine 
compounds, but also for some bromine com- 
pounds. According to the Arrhenius equation, 
the speed of a reaction in which dissociation is a 
primary step is proportional to e~?/””, In our 
experiments aqueous solutions of d-chlor succinic 
acid, and of /-phenyl brom acetic acid were 
prepared. To the first solution sulfuric acid was 
added and to the second perchloric acid was 
added. This was done to repress not only hy- 
droxide ion, but also the ions of the halogen acids 
so as to reduce either inter- or intra-ionic addition. 
Sulfate ion and perchlorate ion are much less 
prone to form addition complexes than other 
negative ions as was shown by Brénsted, 
Kilpatrick and Kilpatrick.? This procedure does 
not eliminate all addition compounds for, as will 
be shown later, substances which have unbonded 
electron pairs can add. The chlor succinic acid 
solution was kept at 85°C for four hours and the 
phenyl brom acetic acid solution was kept at 
room temperature for five hours. In both cases 
the change in rotation was negligible, thus 
showing that primary dissociation does not occur 
to an appreciable extent. To this it might be ob- 
jected that the free ions such as CSH;CtHCOOH 
retain their configurations, but were this the 
case, chloride and bromide ions could have no 
effect except to repress ionization. However, the 
racemizing action of these ions has been com- 
pletely demonstrated not only in this laboratory, 
but also by Holmberg® and by Kekule.’ The 

6 The original experimental work mentioned in the 
various parts of this article will be published separately in 
collaboration with Miss M. J. Young and Mr. F. A. Long. 

7 J. N. Broénsted, Mary Kilpatrick and Martin Kilpatrick, 
J. Am. Chem. Soc. 51, 428 (1929). 


8 B. Holmberg, J. prakt. Chem. 196, 553 (1913). 
® Kekule, Ann. Chem. 130, 21 (1864). 
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reaction of such compounds, therefore, cannot 
proceed by primary dissociation, and so we must 
assume an addition compound or simultaneous 
addition and dissociation. 

The assumption of primary addition com- 
pounds in substitution reactions dates back at 
least to Kekule.* Both Fischer’ and Werner! 
independently assumed addition reactions to 
explain the Walden inversion. Werner, further- 
more, postulated simultaneous dissociation and 
addition. Both authors assumed that the entering 
group might or might not take the place of the 
dissociating group. 

A more definite statement was made by Lewis.” 
“There seems to be but one possible way of 
accounting for this peculiar behavior (Walden 
inversion). Let us consider a carbon atom 
attached to the four radicals R;, Re, Rs and Rg, 
and let us assume that a fifth group, Rs, becomes 
temporarily attached to the carbon atom near to 
the face of the tetrahedron which is opposite to 
R,. Aslight shift of the kernel might make it now 
the center of a new tetrahedron with corners at 
Re, Rs, Rg and R;, while R; would become 
detached from the molecule. Then if the radical 
R;, in the new molecule were to be replaced by 
the radical Rj, the resulting molecule would be 
the mirror image of the one with which we 
started. In this explanation it is not necessary to 
assume that the five radicals are attached to the 
carbon for any appreciable period of time, 
indeed it might be assumed that the R; leaves at 
the same instant that the R; becomes attached to 
the carbon atom.” 

Chemists and theoretical physicists agree that 
a molecule such as R;,R2R;CCI should be repre- 
sented as a tetrahedron. There is, however, this 
important difference; to most chemists the bond, 
i.e., the bond between carbon and chlorine, is 
thought of as extending between the carbon atom 
and the chlorine atom. But Pauling” and Slater" 
have shown by means of quantum mechanics 


that the carbon eigenfunction which is involved ° 


in the carbon chlorine bond extends beyond the 


10 EF. Fischer, Ann. 381, 123 (1911). 
1 A, Werner, Ber. 44, 873 (1911). 
2G. N. Lewis, Valence and the Structure of Atoms and 
Molecules, p. 113, Chemical Catalog Company, N. Y., 1923. 
13, Pauling, J. Am. Chem. Soc. 53, 1367 (1931). 
4 J.C. Slater, Phys. Rev. 37, 481 (1931). 
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MECHANISM OF SUBSTITUTION REACTIONS 


carbon atom on the side away from the chlorine. 
The carbon atom, therefore, has a tendency to 
form a bond at each face of the tetrahedron. The 
stronger the bond is to the apex, the weaker the 
bond-forming power is to the opposite face since 
the strength of a bond is related to the amount 
that the eigenfunctions overlap. It must also be 
remembered that the chlorine in such a com- 
pound has, to some extent, the negative charac- 
teristics of chloride ion, as is shown, not only by 
its eventual ejection from the molecule as 
chloride ion, but also by dipole measurements, 
x-ray diffraction pictures of dichlor ethane, etc. A 
particle with unbonded electron pairs such as a 
chloride or a hydroxide ion will therefore require a 
smaller energy of activation to react with 
R:R:R;CCI if it approaches the carbon face 
directly opposite to the chlorine. It is immaterial 
whether there is addition or simultaneous addi- 
tion and dissociation. The important point is 
that the group which is to be displaced de- 
termines a unique path for the entering group, 
such that the system requires less energy than it 
would for any other path. This, then, leads to 
complete inversion of configuration for a reaction 
which takes place in one step. 


Experimentally the existence of this face 
centered bond is indicated by the work of 


_ Fic. 2. Type of potential energy curve for the dissocia- 
tion of RiR:R:CX into R:R2R:C and X, or into the 
Corresponding ions in the presence of X or X ion. 
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Menke." In order to interpret the x-ray diffrac- 
tion pattern of liquid carbon tetrachloride, 
Menke assumed a structure in the liquid in which 
a chlorine from one molecule of carbon tetra- 
chloride was face centered on the adjoining 
carbon tetrachloride tetrahedron. 

The potential energy diagram of the system 
with the chloride ion added would be of the type 
shown in Fig. 2 where 7 and ¢ would depend on 
the distance AB separating the two chlorines. 
Some of these heats of activation for halogen 
acids and halide ions are now being measured in 
this laboratory. They are about 25 kilocalories 
per mol. 

If there is more than one substance which can 
add we must consider sets of independent 
reactions. Thus in an aqueous solution of 
l-phenyl brom acetic acid, the initial reaction 
probably is 


i-C;H;CHBrCOOH 
(4) 


As the reaction proceeds the reaction 


1-CsH;CHBrCOOH 
+ Br--d-C;H; CHBrCOOH+Br- (5) 


becomes increasingly important. If (5) is very 
fast compared to (4) a racemic mixture only can 
result. If (4) is very fast compared with (5), 
l-mandelic acid never could predominate. Never- 
theless, such cases are known. Thus McKenzie 
and Clough'® report that when /-phenyl chlor 
acetic acid is treated with water or aqueous 
sodium hydroxide, /-mandelic acid results, but if 
it is treated with silver carbonate d-mandelic 
acid is secured. An explanation which suggested 
itself to us and which was in accord with the face 
centered addition theory, was that a third 
independent reaction was occurring which in- 
volved two inversions. Thus if the face of the 
asymmetric carbon atom opposite to the chlorine 
were absorbed on some surface, the carbon- 
chlorine bond would be weakened, causing the 
dissociation of chloride ion and _ inversion. 
Hydroxide ion could now add, again causing 
inversion and releasing the molecule from the 
% H. Menke, Phys. Zeits. 33, 593 (1932). 


16 Alex McKenzie and G. W. Clough, J. Chem. Soc. 
(London) 95, 777 (1909). 
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surface. The net result would be the production 
of /-mandelic acid from /-phenyl chlor acetic 
acid. Thus 


I-CsH;CHCICOOH 
+wall—d-wall complex+ Cl- 


d-wall complex 

No optically active phenyl chlor acetic acid 
was available to us but the corresponding brom 
compound was. This had been previously in- 
vestigated by McKenzie and Walker,’ who 
found the action of silver carbonate on the brom 
compound to be similar to its action on the chlor 
compound. Water, however, produced a mandelic 
acid which was slightly dextro rotatory. Ac- 
cordingly, we dissolved some phenyl brom acetic 
acid in water, divided the sample into two parts, 
and to one part added powdered silica. Both 
samples were then shaken in a thermostat kept at 
25° for several days. The clear sample checked 
the observation of McKenzie and Walker, but 
that to which the silica had been added was levo 
rotatory. 

So far, there are then two ways in which a 
racemic mixture results; (a) reaction (5) is very 
fast compared to reactions (4) and (6); (b) 
reactions (4) and (6) have the same speed; in 
which case the velocity of reaction (5) is not 
critical. For phenyl brom acetic acid it is obvious 
that reaction (5) is slow enough so that we can 
determine that (4) is faster than (6). The 
addition of silica merely made (6) faster than (4). 
In the phenyl chlor acetic acid reaction (6) must 
be faster than (4). 

We can now interpret the function of silver 
salts. The presence of silver ion will keep the 
halide ion concentration low and so prevent 
racemization by reaction (5). This can affect the 
magnitude of the final rotation but not the sign 
of the rotation—that is determined solely by the 
relative rates of (4) and (6). The addition of 
silver ion instead of reducing the surface will 
actually increase it due to the formation of silver 
chloride. Silver ion must therefore increase the 
speed of (4) more than (6). The formation of an 
additicn compound such as 


(6) 


17 Alex McKenzie and Nellie Walker, J. Chem. Soc. 
(London) 109, 1685 (1915). 


A. R. OLSON 


As 
sir: 
C,H;,—C—C 
| OH 
H 


could accomplish this by decreasing the strength 
of the carbon bromine bond and thus increasing 
the speed of the homogeneous reaction. 

In the above discussion we have neglected the 
possibility of ester and ring formation. If an 
ester such as 


Br H 
| 
C,H; — C—COO- 
Br 
| 
H O 


played a prominent part in this reaction, the 
rate law should contain the activity of this 
component to the second power which appears to 
be contrary to experiment. 

The reactions of brom- and chlor-succinic acids 
are somewhat more complex than those of 
phenyl brom acetic acid. Walden!® obtained an 
l-malic acid from /-chlor succinic acid when silver 
oxide was present, and d-malic acid when 
potassium hydroxide was used. Holmberg® and 
Bancroft and Davis!® showed that the rotation of 
the malic acid produced depended upon the 
amount of excess silver oxide. These observations 
can be explained if we assume that in the 
presence of a small amount of silver oxide the 
main effect is to reduce the speed of reaction (5). 
Malic acid and chlor succinic acid which have 
the same rotation must then have opposite 
configurations. In the presence of a large excess 
of silver oxide, not only is (5) slow, but reaction 
(6) is increased. This is in accord with a pre- 
liminary experiment in this laboratory on the 
effect of powdered silica on this reaction. 
Holmberg has also submitted evidence for the 


18 P, Walden, Ber. 32, 1833 (1899). 
19 W. D. Bancroft and H. L. Davis, J. Phys. Chem. 35, 
1251 (1931). 


MECHANISM OF SUBSTITUTION REACTIONS 


formation of a lactone which may hydrolyze in 
two different ways, producing malic acids of 
opposite rotations. Formally, these hydrolysis 
reactions correspond to reactions (4) and (6) 
in that inversion of configuration is produced by 
one substitution on the asymmetric carbon atom, 
and the original configuration is the result of 
two successive substitutions, each of which is 
accompanied by inversion. Various authors have 
considered the possibility of the lack of con- 
figurational change as being due to an even 
number of inversions. 

In the case of the Beckmann rearrangement 
this theory predicts the favorable experimental 
conditions. Thus, if in benzophenone oxime the 
nitrogen oxygen bond is strong enough so that 
primary dissociation does not take place, it 
will be necessary for a substance with an un- 
bonded electron pair, such as chloride ion, to 
approach the nitrogen face opposite to the apex 
to which the hydroxide is bonded, i.e., X in Fig. 3. 
However, because of the free electron pair on the 
nitrogen, such an ion would be subjected to a 
large coulombic repulsion. In a strong acid 
solution hydrogen ion will add to the free pair, 
and reduce this repulsion considerably. This is 


just the condition which Lachman” found was 
necessary for the rearrangement. 

As a further application, consider the thermal 
isomerization of cis dibrom ethylene catalyzed 
by bromine. The usual explanation of this 
reaction is the addition of bromine at the double 
bond followed by a dissociation. Such an 


*° A. Lachman, J. Am. Chem. Soc. 46, 1477 (1924). 


OH 
Fic. 3. 


explanation cannot be correct for tetrabrom 
ethane is rather stable. But if bromine adds to 
the face of the carbon opposite to the apex at 
which a bromine is bound, a reaction corre- 
sponding to a Walden inversion of one of the 
carbons results. The new compound will be 
trans dibrom ethylene. This is the exact analogue 
of reaction (5) which has been experimentally 
demonstrated. 

It is thus seen that this theory of chemical 
reactions, which is based on the energy relations 
involved, and on the quantum mechanical char- 
acter of bonds, does satisfactorily explain the 
results of experiments. In its application to the 
experiments on the Walden inversion it leads to a 
direct correlation between rotation and con- 
figuration. 
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The heat capacity of Ag»SO, has been measured from 
14 to 300°K and the entropy of the salt determined by 
graphical integration. By coupling this value with the 
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entropy of solution and the entropy of silver ion, the 


entropy of sulfate ion in aqueous solution has been 
evaluated. 


HE determination of the entropy of silver 
sulfate is the first! of a series of investiga- 
tions which, along with existing requisite data, 
should serve to establish the entropy of aqueous 
sulfate ion with a good degree of precision. 

The heat capacity of AgeSO, was measured 
from 14 to 300°K by employing the experimental 
method of Latimer and Greensfelder.! Thom- 
sen’s* heat of solution was used and the free 
energy of solution was calculated from solubility 
data. The entropy of Ag* has been recalculated 
by the authors from the most recent data avail- 
able. Twelve independent computations with a 
maximum deviation from the mean of 0.3 e.u. 
give S%Qos.1 of Ag+ = 18.3 e.u. 


MATERIAL 


Silver sulfate of C.P. quality was recrystallized 
twice from concentrated H2SO,. The salt was 
washed and then ignited to 400°C to insure the 
absence of occluded H.SO;. A portion of the 
sample fused to a clear pale yellow liquid. The 
sample showed an Ag content (by precipitation 
of AgCl) of 69.11 percent, whereas theoretical is 
69.19 percent. It is undoubtedly purer than the 
analysis indicates. Specific heat measurements 
were made on a sample weighing 257.444 g. 
(0.82562 mols). 


MEASUREMENTS 


The molal heat capacity values fall on a 
smooth curve which has several points of double 
inflection. These points of inflection undoubtedly 
represent an actual property of the salt since the 
results are reproducible to +0.05 percent. More- 
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over, CaC,O,.H,O! which was measured in the 
same calorimeter and with the same R—T scale 
did not exhibit such a phenomenon. The experi- 
mental results are presented in Table I and 
shown graphically as a function of log 7 in 
Fig. 1. The temperature intervals employed in 
various regions were as follows: below 30°K, 
about 3° per run; 30 to 200°K, about 5° per run; 
above 200°K, from 5 to 8° per run. 


ENTROPY OF SILVER SULFATE 
The entropy of AgsSO, was obtained by graphi- 


- calintegration of theexpression S= 


TABLE I, 


T°K Molal Cp T°K Molal Cy T°K Molal 
14.60 2.12 78.74 17.43 160.89 24.49 
16.18 2.22 82.89 17.96 166.26 24.81 
18.61 3.44 87.18 18.45 171.72 25.18 
21.76 4.64 91.44 18.92 177.95 25.77 
25.32 5.73 95.96 19.38 184.45 26.01 
28.46 6.61 100.89 19.92 191.52 26.41 
32.31 7.76 105.77 20.33 198.96 26.84 
36.77, 9.15 110.29 20.80 206.62 27.34 
40.90 10.42 114.97 21.06 213.81 27.57 
44.78 11.49 119.80 21.52 220.94 27.88 
49.11 12.52 124.38 21.87 229.02 28.23 
53.62 13.28 128.87 22.19 236.88 28.52 
57.41 13.94 133.34 22.43 245.44 28.81 
61.46 14.69 138.00 22.84 255.06 29.35 
65.87 15.17 142.79 23.18 265.41 29.89 
70.03 16.20 147.34 23.61 275.77 30.43 
73.84 16.83 151.74 23.80 285.24 30.80 
75.50 16.98 156.15 24.17 296.20 31.36 


1 Measurements on CaSO4-:2H2O and BaSO, are 
yet unpublished. For complete reference to the general 
theory and previous work, see the authors, J. Am. Chem. 
Soc. 55, 971 (1933). 

2 Thomsen, Thermochemistry, p. 52. Longmans Green 
and Co., 1908. Translated by K. A. Burke. 
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HEAT CAPACITY AND ENTROPY OF SILVER SULFATE 


14 16 16 20 22 24 
Log T 


Fic. 1. The heat capacity of AgsSOs. 


from the curve shown in Fig. 1. The curve 
was extrapolated below 14°K by means of the 
Debye 7* rule, 6=88.01. A summary of the 
entropy calculations is given in Table II. 


TABLE II, 


0.624 e.u. 
47.134 e.u. 


0-14°K 
14-298.1°K 


Extrapolation 
Graphical 


Entropy at 298.1°K of AgsSO;=47.76+0.05 percent. 


HEAT OF SOLUTION 


Thomsen? gives AH of solution at 18°C as 
4480 cal. This value was corrected to infinite 
dilution by means of a value* based on the 
similarity of AgeSO, and K2SO, solutions. This 


*Lange and Streeck, Zeits. f. physik. Chemie A157, 1 
(1931). 

* Randall and Rossini, J. Am. Chem. Soc. 51, 323 (1929). 

*Chloupek and Danes, Collection of Czeckoslovak 
Chemical Communications 4, 124 (1932). 


425 


was then corrected to 25°C by employing the 
partial molal heat capacities of Randall and 
Rossini. We have then: AgsSO,4(s) = AgeSO,(aq) ; 
AFT 098.1 = 4207 cal. 


FREE ENERGY OF SOLUTION 


Chloupek and Danes? find that a saturated 
solution of AgeSO, at 25°C is 0.026888 m. This 
agrees well with the solubility of Harkins and 
others. We have used 0.533" ® for the corre- 
sponding activity coefficient. The free energy 
of solution is given by AF"g95.1= —1363.8 
log [4 (0.533 0.02689)* ] = 6722 cal. 


ENTROPY OF SULFATE ION 


Substituting in the equation AS°=(A//° 
—AF)/T, we have for the reaction, AgeSO,(s) 
=2Ag++S0O,-~— (hypothetical one molal solu- 
tion); AS%o9s.1= —8.4, which yields in turn, by 
application of the third law:—8.4=2X18.3+ 
~—47.8; of SO,- -=2.8 e.u. 

In spite of the uncertainty in correcting 
Thomsen’s? heat of solution to the standard state 
at 25°C, we regard 2.8 as being reliable to 1.0 e.u. 


® Harkins, J. Am. Chem. Soc. 33, 1813 (1911); other 
references in this paper. Hill and Simmons, J. Am. Chem. 
Soc. 31, 821 (1909); Zeits. f. physik. Chemie 67, 594 (1909). 
Rothmund, Zeits. f. physik. Chemie 69, 523 (1910). 
Drucker, Zeits. f. anorg. Chemie 28, 362 (1901). 

7 Lewis and Randall, Thermodynamics and Free Energy of 
Chemical Substances, McGraw-Hill, New York, 1923. 

8 Hass and Jellinek, Zeits. f. physik. Chemie A162, 167 
(1932). 0.533 compares favorably with values obtained for 
similar salts by these workers. 
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In response to requests from our contributors, and after 
a canvass of the opinions of the Editorial Board, the Journal 
of Chemical Physics wishes to announce that a “Letters to 
the Editor’ section will be carried by the Journal. This 
section will accept reports of new work, provided they are 
terse and contain few figures, and especially few halftone cuts. 
The Editorial Board will not hold itself responsible for 


The discovery! that water in which the heavier isotopes 
of Hz and O: have been concentrated by electrolysis has 
a higher boiling point than normal water shows that it 
must be possible to fractionate water by distillation. To 
demonstrate this, 10 liters of water having a specific 
gravity of 1.000053 were distilled at atmospheric pressure 
in a still provided with a 35-foot rectifying column. An 
initial distillate of 200 ml and the final residue of 100 ml 
were compared as to density and found to differ by 64.9 
parts per million, the residue having increased by 53.3 
parts and the distillate having decreased by 13.2 parts. 
Distillation fractionation is thus possible and should find 
practical application in combination with electrolysis 
fractionation. 

The adsorptivity of water by porous materials or 
powders should vary with its isotopic composition and it 
should therefore be possible to fractionate water by 
adsorption. To demonstrate this 300 g of activated 
charcoal was allowed to stand for three weeks immersed 
in 500 g of water having a specific gravity of 1.000053. 
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opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 
not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual page charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 


The supernatant water (176 ml) and the last 123 ml of 
water which could be expelled from the charcoal by 
heating in vacuo at 500°C were compared as to density 
and found to differ by 12.4 parts per million, the super- 
natant water having decreased by 6.5 parts and the ad- 
sorbed water having increased by 6.7 parts. 

The specific gravities reported above were determined 
on carefully purified samples of water and are based upon 
normal water taken as unity. 

Epwarp W. WASHBURN 
EpGarR R. SMITH 
Bureau of Standards, 
Washington, D. C., 
May 15, 1933. 


* Publication approved by the Acting Director, Bureau 
of Standards, Department of Commerce. 

1 Washburn, Smith and Frandsen, J. Chem. Phys. 1, 
288 (1933). 
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